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ABSTRACT 


This  work  deals  with  the  experimental  and  theoretical  investigation 
of  53Cr  and  55Fe.  Excited  states  of  53Cr  were  studied  by  means  of  the 
52Cr (d ,py) 5 3Cr  reaction.  In  the  case  of  55Fe,  the  55Mn(p,ny)  reaction 
was  used.  The  theoretical  calculations  based  on  Nilsson  model  and  the 
unified  model  are  compared  with  experiment  and  other  available  calcula¬ 
tions  . 

The  decay  modes  of  many  of  the  excited  states  of  5  3Cr  below  4.25  MeV 
excitation  energy  have  been  determined  by  means  of  the  52Cr (d ,py) 5 3Cr 
reaction.  Previous  tentative  spin  assignments  of  5/2",  3/2“  and  9/2+ 
to  the  1.006  MeV,  2.324  MeV  and  3.715  MeV  levels  respectively  have  been 
confirmed.  New  spin  assignments  are  the  following:  5/2  or  9/2  to  the 
2.227  MeV  level,  3/2"  to  the  3.190  MeV  level,  5/2  to  the  3.271  MeV 
level  and  1/2  or  3/2  to  the  3.971  MeV  level.  Multipole  mixing  ratios  are 
given  for  many  of  the  observed  transitions. 

The  study  of  the  55Mn(p,ny) 55Fe  reaction  resulted  in  the  determination 
of  decay  modes,  spin  assignments  and  multipole  mixing  ratios  for  levels 
of  55Fe  between  2.5  and  3.1  MeV  excitations.  New  spin  assignments  are 
the  following:  11/2",  (9/2)  for  the  2.542  MeV  level,  9/2  for  the  2.984 
MeV  level  and  (11/2,  9/2)  for  the  3.076  MeV  level.  Measurements  are  in 
agreement  with  the  3/2“  assignment  of  the  3.027  MeV  level. 

The  Nilsson  model  was  used  to  predict  energy  levels  of  53Cr  and  55Fe 
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but  the  fit  to  experimental  results  is  only  fair.  Calculations  based  on 
the  intermediate  coupling  in  the  unified  model  yielded  energy  levels, 
electromagnetic  transition  rates,  branching  ratios,  mixing  ratios,  life¬ 
times  and  spectroscopic  factors. 

These  calculations  resulted  in  an  overall  agreement  with  experimental 
results  for  most  of  the  levels  of  53Cr  and  55Fe  below  ^  2.5  MeV  excitation. 
Comparison  is  also  made  with  other  available  calculations  based  on  the 
unified  model,  the  "Thankappan  and  True"  model  and  the  shell  model. 
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CHAPTER  I 


MOTIVATION 


From  the  measurements  of  electromagnetic  properties  in  nuclei,  it  is 
possible  to  determine  the  modes  of  decay  of  nuclear  levels,  their  spins, 
parities,  electromagnetic  transition  rates,  static  moments,  .  .  .  without 
reference  to  a  particular  theory  of  nuclear  interactions.  These  measure¬ 
ments  constitute  a  good  test  for  the  predictions  of  nuclear  models. 

In  the  last  decade,  there  has  been  great  interest  shown  in  the  study 
of  odd-A,  N  =  29  nuclei  as  they  have  a  simple  v  2p3/2  configuration  and 

are  expected  to  be  easily  described  by  theoretical  calculations.  There 

49  51  53  55 

are  four  nuclei  of  interest:  Ca,  Ti,  Cr  and  Fe.  All  four  nuclei 

20  22  24  26 

have  been  studied  by  (d,p)  stripping  reaction  on  the  nearest  even  nucleus 
(Ka  64,  Ba  64,  Ra  68,  Ma  64  and  references  therein).  These  studies  have 
resulted  in  the  determination  of  energy  levels  up  to  6  to  8  MeV  excitation. 
Comparatively,  there  has  been  little  work  done  to  determine  gamma-ray 
decay  scheme  and  spin  assignments  of  these  levels.  This  explains  why 
theoretical  calculations  were  mainly  concerned  with  reproducing  energy 
levels  and  spectroscopic  factors  rather  than  with  the  electromagnetic 
properties . 

Many  factors  must  be  taken  into  consideration  before  doing  gamma-ray 
work:  target-making  feasibility,  Q-value  of  the  reaction,  problem  of 
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contaminants,  competing  reactions,  gamma-ray  yield  are  among  the  most 
important.  Their  evaluation  will  permit  to  decide  whether  or  not  an 
experiment  is  feasible  and  also  what  type  of  method  is  to  be  used  to 

extract  as  much  information  as  possible. 

4  ft 

°Ca  targets  are  difficult  to  make  and  oxidize  very  easily  even  under 
good  vacuum.  The  presence  of  a  large  amount  of  oxygen,  which  is  not  a 
major  hindrance  for  study  of  stripping  reactions,  virtually  eliminates 
all  possibility  of  the  study  of  the  (d,py)  reaction, if  coincidence  require¬ 
ments  are  imposed,  as  most  of  the  y-yield  originates  from  the  contaminant. 
Self-supporting  targets  of  50Ti  are  easy  to  prepare  and  the  metal  does  not 
oxidize  too  easily.  Unfortunately,  the  gamma  yield  is  too  low  to  permit 
any  p-y  coincidence  study  of  the  (d,py)  reaction.  The  study  of  53Cr  and 
55Fe,  on  the  other  hand,  proved  to  be  quite  feasible  and  is  the  object  of 
the  present  work. 

The  stripping  reaction  52Cr (d ,p) 5 3Cr  has  been  extensively  studied  in 
the  past  few  years  (An  64,  Bo  64,  Ro  65  and  Ra  68).  However,  only  the 
decay  modes  of  levels  below  2.5  MeV  had  been  determined  (Ro  65,  Ra  68) 
prior  to  this  work.  Assignments  of  the  spins  of  the  excited  states  of 
53Cr  had  been  based  on  (n,n'y)  and  (n,y)  experiments  (Va  62,  Ba  65)  and 
on  the  empirical  rules  in  the  analysis  of  pick-up  and  stripping  reactions 
(Ra  68,  Wh  67).  Most  of  these  assignments  were  tentative. 

The  level  scheme  of  33Fe  has  been  determined  from  (d,p)  reactions 
(Fu  63,  Ma  65)  and  little  has  been  done  in  the  study  of  the  gamma-ray 
decay  until  the  recent  work  of  Pilt  (Pi  69)  who  determined  the  spins  and 
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gamma  decay  of  most  of  the  levels  below  2.6  MeV  excitation,  using  the 
55Mn(p,ny) 55Fe  reaction.  The  shell  model  predictions  of  Ohnuma  (Oh  66) 
appeared  to  be  in  very  good  agreement  with  Pilt’s  results.  These  cal¬ 
culations  also  predicted  two  high-spin  states  below  3.0  MeV  excitation 
(11/2  at  2.790  MeV  and  13/2  at  2.975  MeV)  not  reported  by  Pilt.  It 
was  therefore  decided  to  extend  the  study  of  55Fe  up  to  the  3  MeV  region 
in  order  to  investigate  the  possibility  of  the  existence  of  such  high- 
spin  states. 

There  are  two  main  approaches  to  the  theoretical  study  of  odd-A 
N  =  29  nuclei.  First,  shell  model  calculations  in  which  48Ca  is  consid¬ 
ered  as  an  inert  core.  Effective  p-p  and  p-n  interactions  are  used  to 
describe  the  different  configurations  outside  this  core  (Ve  66,  Oh  66). 

The  second  approach  is  based  on  the  intermediate  coupling  model  and 
considers  the  coupling  of  the  odd  nucleon  to  the  neighbouring  even  core, 
which  is  either  in  its  ground  state  or  in  an  excited  state.  If  the 
nature  of  the  excited  states  of  the  core  is  not  supposed  to  be  known,  the 
so-called  "Thankappan  and  True"  approach  is  used  (Th  65) .  Such  calcula- 

r-  o 

tions  were  quite  successful  in  describing  the  low  lying  levels  of  D;)Cu 
(Th  65) .  Recently,  this  type  of  calculation  was  used  to  describe  the 
N  =  29  odd  nuclei  (Ph  70,  La  70).  In  a  more  specific  approach,  one 
considers  the  even  core  as  capable  of  performing  quadrupole  oscillations. 
This  approach  was  used  by  Ramavataram  (Ra  63)  to  describe  N  =  29  odd  nuclei. 

In  the  calculations  of  Ramavataram,  only  levels  with  spins  up  to  7/2 
were  taken  into  consideration.  Moreover,  these  calculations  were  done  at 


' 
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a  time  when  very  little  was  known  experimentally  about  53Cr  and  55Fe. 
Similarly  the  shell  model  calculations  of  Ohnuma  (Oh  66)  and  Vervier  (Ve  66) 
dealt  mainly  with  predictions  of  energy  levels  and  calculations  of  spectro¬ 
scopic  factors.  In  order  to  attempt  a  quantitative  appraisal  of  the  merits 
of  the  respective  models,  more  complete  calculations  were  needed.  Such 
calculations  based  on  the  intermediate  coupling  model  were  undertaken  in 
this  work  while  existing  shell  model  calculations  were  extended  by  Ohnuma 
(Oh  70). 

If  the  odd  nucleon  is  strongly  coupled  to  a  deformed  core,  Nilsson 
model  calculations  are  then  applicable.  Such  calculations  have  been 
fairly  successful  in  describing  57Fe  and  5gFe  (So  67,  So  69).  As  52Cr  and 
5,4Fe  are  not  strongly  deformed  (see  e.g.  Go  62),  Nilsson  model  calculations 
for  5^Cr  and  5^Fe  have  very  little  chance  of  success.  However,  it  was 
thought  that  such  calculations  might  yield  some  fair  results  if  band  mixing 
was  included,  and  they  were  consequently  attempted. 
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CHAPTER  II 


EXPERIMENTAL  INVESTIGATION  OF  53Cr 

The  spins  and  mixing  ratios  for  many  of  the  levels  below  4.3  MeV 
excitation  in  53Cr  were  obtained  using  the  52Cr (d ,py) 53Cr  reaction, 
by  application  of  Method  II  of  Litherland  and  Ferguson  (Li  61) ,  with 
detection  of  the  protons  at  approximately  180°.  Branching  ratios  were 
also  measured  in  this  experiment.  Lifetimes  were  determined  via  the 
53Cr (p , p ’ y) 5 3Cr  reaction  using  the  Doppler  Shift  technique,  but  the 
results  obtained  were  rather  limited. 

Section  1  describes  the  theoretical  aspects  involved  in  this 
study.  Detailed  theoretical  treatment  of  the  study  of  gamma-ray 
angular  distributions  can  be  found  in  articles  by  Litherland  and 
Ferguson  (Li  61),  Poletti  and  Warburton  (Po  65)  and  Rose  and  Brink 
(Ro  67).  Consequently,  the  presentation  of  this  theory  will  be  limited 
to  general  principles.  The  notation  and  sign  conventions  of  Rose  and 
Brink  will  be  used  throughout.  Section  2  describes  the  experimental 
procedure  associated  with  the  study  of  angular  distributions.  A 

description  of  the  Doppler  Shift  Attenuation  Method  can  be  found  in 

•\ 

an  article  by  Litherland  (Li  64) .  This  method  will  be  outlined  in 
section  3.  The  experimental  results  are  given  in  detail  in  section  4. 
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1 .  THEORY 


1.1.  Gamma-ray  Distributions  for  Decay  of  Aligned  States 

Consider  the  reaction  A(a,b)B  where  neither  the  target  nucleus 
A  nor  the  beam  of  particles  a  is  polarized.  It  is  also  assumed  that 

JL  I 

the  residual  nucleus  B  is  an  excited  state  iJ^m-^  with  definite  parity. 
The  state  Ij^m^  is  said  to  be  aligned  if  the  population  parameters 
PCn^)  [probability  that  the  m^h  substate  is  populated]  satisfy  the 
relation  P(m  )  =  P(-m  ).  Because  the  state  |j  m  >  has  definite  parity, 
there  is  symmetry  under  reflection  through  a  plane  normal  to  the  beam 
direction  (Fe  65).  Under  these  conditions,  (unpolarized  target  and 
beam,  intermediate  state  of  definite  parity) ,  there  is  alignment  with 
respect  to  the  beam  axis  as  long  as  the  beam  axis  remains  an  axis  of 
cylindrical  symmetry. 

Consider  the  state  Ij^m^  decaying  by  emission  of  a  gamma  ray  to 
a  state  |j  m2>.  If  q  defines  the  circular  polarization,  (q  =  +1)  and 
if  it  is  the  wave  vector  of  the  gamma  ray,  the  probability  amplitude 
for  such  a  transition  may  be  represented  by: 


Aq  (t) 


m<TT> 

Lm 


J2m2> 


L 

mo 


(R) 


(II. 1) 


, 

’ 
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Where  ^  (R)  is  the  rotation  matrix  describing  the  rotation  R  =  (a$Y) 


mq 


#  ^  \  TT 

which  takes  the  beam  axis  z  into  the  direction  k  5  the  operator  Tt 

Lm 


is  the  electromagnetic  operator  of  parity  tt  [tt  =  0  for  electric 
and  1  for  magnetic  multipolarity] . 

If  the  orientation  of  the  spin  of  the  final  state  is  not  observed, 
then  the  probability  of  transition  from  the  state  |j^nij>  to  any  substate 
of  J2  by  emission  of  a  photon  is 


(II. 2) 


Assuming  that  there  is  cylindrical  symmetry,  the  total  transition 
probability  is  then  equal  to  equation  (II. 2)  weighted  by  the  relative 
probability  of  population  of  each  substate  of  the  initial  state  jj^m^ 
i.e.  it  is  given  by: 


l  p(V  l 


(II. 3) 


This  expression  is  proportional  to  the  number  of  gamma  rays  emitted 
in  the  direction  k.  For  its  evaluation,  equation  (II. 1)  is  substituted 
into  (II. 3).  In  the  resulting  expression,  the  product  of  rotation  matrices 


a,  3  and  y  are  the  Euler  angles  of  the  rotation. 
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can  be  expressed  in  terms  of  expansions  in  Legendre  Polynomials;  the 
operator  matrix  elements  are  reduced  by  application  of  the  Wigner-Eckart 
theorem.  Finally,  the  sums  involving  different  Clebsch-Gordan  coefficients 
are  reduced  by  techniques  of  Racah  Algebra.  The  resulting  angular  distri¬ 
bution  is: 


W(J  +  J  ;kq)  =  £  l  B  (J  )  p  (cose )  (-)q+Ji  J2+L’‘L‘k 

kLL '  k  k 

TT7T  ’ 


x  /2J  +1  (LqL'-q  |  kO)  W(JX  J^L ’  ;kJ2) 


MJI  t‘711  j2>* 


(II. 4) 


in  which: 


Bk  (Jl)  =  l  p  (m! )  (-)J!  (JjmjJj-mJ  kO) 

ml 

There  are  limitations  of  the  multipole  order  L  of  a  given  transition 
->  j  namely:  |  J^-jJ  <_  L  <  J^+J^  (conservation  of  angular  momentum)  and 
L  ^  0  (intrinsic  spin  of  photon  is  1).  Also,  for  a  given  L,  tt  is  restricted 
to  electric  or  magnetic  because  parity  is  conserved  in  this  interaction. 

In  general,  only  the  lowest  orders  contribute  to  a  transition  and  it  is,  in 
most  cases,  sufficient  to  consider  the  two  lowest  multipole  orders  L  and  L'. 
One  then  defines  the  multipole  mixing  ratio  6  as: 


. 


n  • 


. 
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6 


2L’+1 


2L  +1 


L’  >  L 


(II. 5) 


In  this  study,  polarization  is  not  observed.  One  can  then  write  the 
angular  distribution  formula  as: 


R  (LLJ  J  )+  26R,  (IX  'J  J  )+62R,  (LLJ  J  ) 
W(0)  -l  Pk(cos0)  \  +  h  - ^ - — 


(II. 6) 


where : 


R^LL'J  J  )  =  (-)1+J1  J2+L'  L  k  / (2J  +1)  (2L+1)  (2L’+1)  (LlL'-l|kO) 


x  W(J^  J^LL '  ;kJ? ) 


(II. 7) 


Remarks : 

a)  In  equation  (II. 6),  the  B  (J  )  terms  contain  all  the  information 

K-  1 

on  the  alignment  of  the  initial  state  |j1m1>.  The  R^.  terms  depend  specif¬ 
ically  on  the  J1  ->  J2  gamma-ray  cascade.  The  nuclear  properties  on  the 

other  hand ,  appear  exclusively  through  the  mixing  ratio  6, 

b)  In  the  Racah  coefficients  of  equation  (II. 7),  triangular  relations 
limit  k  to  k  <  min  (2l/,  2Ji). 

c)  As  a  consequence  of  alignment 

Bk(Jl)  =  0  If  k  is  odd.  Therefore,  the  angular  distribution  is  limited  to 


r-'Ji.;)  i  ;=>  ;  '  '  "  •  \  . 
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even  order  Legendre  polynomials  and  is  consequently  symmetric  about  90°. 

d)  Equation  (II. 6)  is  established  for  an  ideal  point  detector. 

However,  because  of  the  finite  size  of  the  gamma  detector,  the  distribution 
is  "averaged-out "  over  the  solid  angle  the  detector  subtends.  For  a 
cylindrical  detector  this  attenuation  may  be  corrected  by  introducing 
coefficients,  Q^,  for  each  order  of  Legendre  polynomial.  In  equation  (II. 6), 
each  P^  (cos0)  is  to  be  replaced  by  Q^P^.  (cos0).  For  standard  size  gamma 
detectors  these  coefficients  have  been  tabulated  [see  e.g.  (Wa  59)]. 

1.2.  Alignment  in  the  Case  of  Triple  Correlation 

In  order  to  preserve  alignment  in  the  A(a,b)B  reaction,  the  particles 

are  detected  at  either  0°  or  180°  with  respect  to  the  beam  axis  in  the  case 

of  triple  correlation.  Then,  strong  limitations  occur  in  the  number  of 

magnetic  substates  that  are  populated  (Li  61) .  Because  the  orbital  angular 

momenta  of  the  incoming  and  outgoing  particles  are  perpendicular  to  the 

beam  axis,  only  the  intrinsic  spins  of  the  various  particles  involved 

contribute  to  the  population  parameters.  That  is,  the  largest  projection 

& 

which  can  be  obtained  in  B  is  m  =  J.  +  s  +  s,  where  J.,  s  and  s, 

max  A  a  b  A  a  b 

are  the  intrinsic  spins  of  the  target,  the  incoming  and  outgoing  particles 

respectively.  In  the  52Cr (d,py) 53Cr  reaction,  =  0,  s^  =  1  and  s^  =  V2» 

therefore  m  =  3/?  i.e.  only  P(3/2)  and  P(l/2)  parameters  are  non-zero, 
max 

Because  of  the  finite  size  of  the  detector,  which  in  this  experiment, 
was  an  annular  silicon  surface  barrier  detector  situated  at  180c  with 


j-J 

' 


■tdlaiJiBq  &  j  j;  lox.  a(d«fi)A  ril  ni  Jnsti  i|  :Xfi  avia  tq  03  isbio  n~ 

■ 


11 


respect  to  the  beam  axis,  substates  beyond  m  can  actually  be  populated. 

in  3.x 

Litherland  and  Ferguson  (Li  61)  have  shown  that  if  E,  is  the  half  angle 
subtended  by  the  detector  (£  in  radians),  substates  (m  +  n)  are  populated 

II13X 

2n 

proportionally  to  £  .  For  reasonable  detector-target  distance,  this  effect 

is  generally  negligible  as  £  is  small. 

1.3.  x2  Fitting  of  Angular  Distributions 

As  shown  earlier,  in  (d,p)  experiments  on  even  nuclei  where  the  protons 
are  detected  at  180°,  only  the  magnetic  substates  m  =  3/2  and  m  =  V2  are 
populated.  The  theoretical  angular  distribution  of  the  gamma  rays  is 
expressed  in  terms  of  the  population  parameters,  the  mixing  ratio  6  and  the 
spin  sequence.  A  linear  least-squares  fit  of  the  theoretical  distribution 
to  the  experimental  distribution  is  made  for  a  discrete  set  of  fixed  values 

of  6.  For  each  value  of  6,  the  best  fit  corresponds  to  those  values  where 

2 

X  is  given  by: 

X2  =  -  l  [Y(0  )  -  W(0  )]2  /  E2(6  ) 
n  0  1  1  1 

in  which  E  (0  )  is  the  uncertainty  assigned  to  the  gamma-ray  yield  Y(0^) 
at  angle  0  ,  W(0  )  is  the  theoretical  yield  and  n  is  the  number  of  degrees 
of  freedom,  (number  of  data  points  minus  number  of  free  parameters  in  W(0)). 

Since  6  varies  between  and  +00,  it  is  more  convenient  to  step  through 

2 

the  finite  range  values  of  arc  fan  6.  The  resulting  plots  of  X  versus 
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<5  show  dips  which  correspond  to  possible  solutions  for  at  least  one  set 
of  the  allowed  spins  J1  and  J?.  Sets  of  6,  and  J2  which  correspond 
to  a  x2  lower  than  an  assigned  confidence  limit  are  kept  as  possible 
solutions.  The  confidence  limit  is  a  measure  of  the  probability  that, 
for  the  correct  theoretical  expression,  X  will  exceed  a  certain  value. 
For  example,  with  5  degrees  of  freedom,  the  1%  confidence  limit  is 
3.0.  This  means  that  X  is  expected  to  be  larger  than  3.0  in  1%  of 
the  measurements  if  the  correct  theory  is  used  (in  the  present  case, 
the  correct  spins  and  mixing  ratio).  The  0.1%  confidence  limit 
criterion  is  usually  applied  in  the  determination  of  spins  and  mixing 
ratios,  i.e.  solutions  for  which  X  exceeds  the  value  at  the  0.1%  limit 
are  rejected.  The  0.1%  confidence  limit  criterion  was  adopted  in  this 
work. 


1 . A .  Restriction  on  the  Population  Parameter  Ratio 

The  technique  described  above  takes  no  account  of  the  physical  reality 

of  a  state  except  in  setting  a  population  parameter  to  0  when  a  negative 

value  is  obtained  for  it  in  the  fitting  procedure.  The  question  arises 

then,  whether  all  possible  magnetic  substate  populations  represent  physical 

situations.  This  has  been  discussed  by  Goldfarb  and  Wong  (Go  66)  who  have 

p (3/2) 

calculated  the  ratio  of  population  parameters  £  =  Y(±/2)  in  the  CaSe  °f 
the  2.324  MeV  level  of  53Cr  excited  via  the  52Cr(d,p)  reaction  for  several 
bombarding  energies  between  6.6  and  11  MeV.  They  have  shown  that,  when 
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protons  are  detected  at  180°,  £  varies  between  0.31  and  1.11  for  the 
direct  interaction  mechanism,  whereas  £  remains  almost  constant  and 
equal  to  1.10  for  the  compound  nucleus  mechanism.  In  this  experiment, 
many  weak  non-stripping  states  were  observed  which  are  formed  pre¬ 
dominantly  through  the  compound  nuclear  process.  The  experimental  yield 
curve  showed  no  strong  resonance  structure  and  narrow  resonances  would 
have  been  integrated  over  by  the  thick  target.  Hence  it  was  considered 
that  physically  realistic  situations  would  result  from  limits  on  £  from 
0.3  to  3.0.  In  section  4,  this  restriction  is  used  to  reject  certain 
spin  assignments  for  two  states  in  53Cr. 
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2 .  EXPERIMENTAL  PROCEDURE 


2.1.  Apparatus 

Deuterons  from  the  University  of  Alberta  CN  Van  de  Graaff  accelerator 
were  used  to  bombard  150  yg/cm2  thick  self-supporting  targets  of  chromium 
enriched  to  99.9%  in  52Cr.  Natural  chromium  targets  were  also  used  for 
preliminary  runs  at  4.85  MeV  bombarding  energy.  On  the  basis  of  the 
yield  curve  measurements  between  4.0  and  5.5  MeV,  deuteron  energies  of 
4.11  MeV  and  4.85  MeV  were  chosen  for  the  experiment  since  they  correspond 
to  optimum  excitation  of  the  low  lying  levels  of  53Cr.  The  beam  current 
was  maintained  at  about  200  nA.  A  series  of  three  collimators  lined  with 
platinum  determined  a  beam  spot  of  about  1.5  mm  diameter  on  target. 

The  target  chamber  is  described  in  fig.  1.  The  protons  were  detected 
by  an  annular  silicon  surface  barrier  detector1  located  at  180  with 
respect  to  the  beam  direction.  This  detector,  which  has  a  resistivity  of 
20,000  fi-cm,was  operated  at  a  bias  of  100  volts  and  the  resolution  was 
better  than  25  keV  in  the  region  of  interest.  A  shield  on  the  detector 
mount  (fig.  1)  prevented  elastically  backscattered  deuterons  from  the 
tantalum  beam  stop  from  striking  the  detector  face. 

ft 

The  gamma  rays  were  detected  by  a  7.5  cm  x  7.5  cm  Nal(Tl)  detector 


^  Nuclear  Diode  Inc.  ND  PHG-35-14 
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Figure  1  Schematic  diagram  of  target  chamber  and  detector  assembly 
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placed  10  cm  from  the  target  and  shielded  by  2  mm  thickness  of  lead  to 
eliminate  low  energy  background. 

The  system  isotropy  was  checked  by  measuring  the  angular  distri¬ 
bution  of  the  3.09  MeV  transition  (l/2+  to  1/2“  ground  state)  in  13C 
via  the  13C(p,p')  reaction. 

A  schematic  diagram  of  the  fast-slow  coincidence  of  the  electronic 
system  is  shown  in  fig.  2.  Fast  timing  signals  provided  by  fast  zero 
crossing  discriminators  for  both  the  gamma  and  the  proton  lines  were 
fed  into  a  time-to-amplitude  converter  (TAC) .  The  time  peak  correspond¬ 
ing  to  the  true  coincidence  events  was  approximately  15  nanoseconds 
f.w.h.m.  Two  windows  of  equal  widths  were  set  on  the  time  spectrum, 
one  of  the  time  peaks  corresponding  to  the  "trues"  plus  "randoms"  (T  +  R) , 
the  other  on  an  adjacent  region  corresponding  to  the  "randoms"  (R) .  Both 
gamma  and  proton  lines  were  self-gated  in  order  to  eliminate  the  low 
energy  part  of  the  spectra  and  thus,  reduce  the  counting  rates  and  their 
undesirable  effects.  The  general  coincidence  gate  was  generated  each 
time  proton  and  gamma-ray  pulses  were  in  coincidence  with  either  the 
"trues  +  randoms"  or  the  "randoms"  signal.  The  other  gate  used  for  rout¬ 
ing  was  provided  by  the  "randoms"  signal  only. 

The  linear  signals,  from  the  gamma  and  proton  detectors,  were  fed 
into  two  analog-to-digital  converters  (ADC’s)  .  These  two  ADC's  were 
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Figure  2  Block  diagram  of  coincidence  electronics  used  in  the  correlation 
measurements.  Abbreviations  are  explained  in  the  text.  Most  of 
the  electronic  equipment  was  obtained  from  Ortec  ^  .  Ortec 
model  numbers  are  used. 
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4* 

connected  to  the  on-line  SDS  920  computer,  as  was  a  third  ADC  which 
recorded  the  TAC  spectrum. 

The  on-line  computer  sorted  the  gamma  rays  according  to  charged- 
particle  pulse  height.  This  was  accomplished  by  setting  windows 
around  the  charged-particle  peaks  of  interest  and  routing  coincident 
gamma  rays  into  separate  regions  of  computer  memory,  according  to 
which  window  their  corresponding  charged  particle  pulse  belonged.  Up 
to  16  such  charged  particle  bins  could  be  set  with  each  bin  consisting 
of  256  channels  for  true  plus  random  coincidences,  and  256  channels 
for  randoms  only.  The  sorting  routine  also  provided  1024  channels  for 
the  total  coincident  charged  particles,  1024  channels  for  total  random 
charged  particles,  256  channels  for  the  total  coincident  gamma  ray 
trues  plus  randoms,  and  256  channels  for  the  random  gamma  rays  only. 

A  self-gated  proton  single  spectrum  was  recorded  for  monitoring 
purposes.  The  relative  intensities  of  the  excited  states  to  the 
ground  state  of  53Cr  remained  constant,  and  hence  the  latter  proton 
peak  was  used  to  normalize  the  different  runs. 

The  gamma  rays  were  detected  at  angles  of  0,  30,  45,  60  and  90 

degrees.  For  the  first  series  of  measurements  at  =  4.85  MeV,  only 

the  0,  45  and  90  degree  positions  were  repeated.  For  the  second  series 

at  E,  =  4.11  MeV,  each  angle  was  repeated  3  times, 
d 
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2.2.  Analysis  of  the  Gamma  Ray  Spectra 

A  slight  gain  shift  was  experienced  in  the  gamma  ray  system  during 
the  course  of  this  experiment.  This  was  overcome  by  using  a  computer 
program  that  normalized  all  the  spectra  to  a  standard  gain.  This 
program,  called  LYNE,  (Ca  69)  fits  two  well-defined  peaks  with  a 
gaussian  shape,  determines  accurately  the  peak  positions,  then  renormal¬ 
izes  the  complete  spectrum  to  any  given  gain.  The  different  runs  were 
then  weighted  by  sin  6  and  added  to  give  good  statistics  for  the  deter¬ 
mination  of  branching  ratios. 

The  intensities  of  the  different  gamma  transitions  were  determined 
using  a  gamma  line-shape  analysis  program  GAMRAY  (Sy  68) .  This  program 
generates  the  full  shape:  photopeak,  escape  peaks  and  Compton  distri¬ 
bution  for  each  gamma  ray.  Peaks  and  Compton  distribution  are  smeared 
by  a  near-gaussian  function  and  allowance  is  made  for  absorption  in  the 
crystal  of  Compton  scattered  gamma  rays.  The  parameters  governing  the 
peak  widths,  photopeak-to-escape  peak  ratio  and  photopeak-to-Compton 
ratio  were  reevaluated  using  the  3.564  and  3.48  MeV  gamma  rays  from  the 
26Mg(d,py)27Mg  reaction.  The  absorption  and  efficiency  parameters  were 
checked  against  the  decay  of  2t+Na. 

2.3.  Angular  Distribution  Analysis 

The  gamma  ray  angular  distributions  were  fitted  according  to  the 

O 

theory  described  in  section  1.1  using  a  X  computer  program  (Hu  68). 
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The  distributions  are  fitted  to: 

W(0)  =  +  a2P^(cos0)  +  ai+P^(cos0) 

in  order  to  get  the  best  fit  for  the  Legendre  Polynomial  coefficients  a^. 

These  coefficients  are  also  fitted  using  equation  (II. 6)  [according 
to  the  method  described  in  section  1] .  Angular  distributions  for  several 
transitions  involved  in  the  decay  of  the  initial  state  |  >  can  be  fitted 

simultaneously . 

Finally,  when  limitations  on  the  population  parameters  were  intro¬ 
duced  (as  described  in  section  1.4.),  a  modified  version  of  a  grid  search 
program  (So  68)  was  used.  This  program  is  basically  equivalent  to  the  X 
program  already  described,  and  is  based  on  the  same  formalism  (Ro  67). 
However,  it  was  found  easier  to  modify  for  that  purpose. 


. 
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3.  LIFETIME  MEASUREMENTS  IN  53Cr 

The  53Cr (p ,p ' y) 5 3Cr  reaction  was  used  together  with  Doppler  shift 
measurements  in  order  to  determine  the  lifetimes  of  the  low-lying 
states  of  53Cr.  The  Doppler  shift  attenuation  method  (DSAM)  will  just 
be  outlined  here.  It  has  been  used  extensively  in  the  last  few  years 
and  is  thoroughly  described  in  many  papers  [see  e.g.  (Li  64)  and 
references  therein] . 

The  Doppler  shift  in  the  energy  of  a  gamma  ray  emitted  by  a  nucleus 
moving  with  velocity  v(t)  is,  to  the  first  order  in  v(t)/c: 

E  =  E  (1  +  cos0 )  (II. 8) 

0  c 

where  EQ  is  the  unshifted  gamma-ray  energy  and  0  is  the  angle  of  emission 
of  the  gamma  ray  relative  to  v(t) . 

Since  the  nuclei  are  slowed  down  in  the  target  or  in  the  target 
backing,  the  velocity  is  a  function  of  time,  and  the  Doppler  shift  is 
smaller  than  the  maximum  shift  (v(o)/c  cos0)EQ.  Assuming  exponential 
law,  equation  (II. 8)  becomes: 

E  =  En (1  F(x)  cos0} 

0  c 


where  F(x)  is  called  the  attenuation  factor  and  is  measured  for  a  given 


(6.  ii)  ^  (eaoo^^+i)  0a  =  3 
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stopping  material,  x  is  the  lifetime  of  the  excited  nuclei. 
The  attenuation  factor  F(x)  is  given  by: 


F(t) 


Observed  average  Doppler  Shift _  1 

Full  shift  ~  v(0)x 


v(t)  e 


-t/x 


dx 


One  can  see  that: 

F(x)  =  1  {or  the  observed  average  Doppler  shift  is  equal  to  full  shift} 
when  the  gamma  emission  takes  place  before  the  recoil 
nucleus  is  actually  slowed  down.  This  corresponds  to 
either  a  light  stopping  material  or  a  very  short  life- 
t ime  x . 

F(x)  =  0  {or  the  observed  average  Doppler  shift  is  zero}  when 

the  gamma  emission  takes  place  after  the  recoil  nucleus 
has  been  stopped.  This  corresponds  either  to  a  heavy 
stopping  material  or  long  lifetime. 

The  interpretation  of  the  measured  attenuation  of  the  Doppler 
shift  requires  a  knowledge  of  the  velocity  v(t)  of  the  recoiling 
excited  nuclei  as  a  function  of  time.  In  this  experiment,  the  attenua¬ 
tion  factor  F(x)  was  calculated  using  a  program  called  DSAM  III  based 
on  the  theoretical  calculations  of  Blaugrund  (B1  66) .  The  experimental 
F(x)  was  then  compared  with  the  calculated  F(x). 

The  targets,  approximately  200  pg/cm-2  thick,  consisted  of  chromium 
enriched  to  96.4%  in  53Cr  which  has  been  evaporated  onto  a  tantalum 
backing  0.005  cm  thick.  The  gamma  rays  were  detected  in  a  45  cm3  Ge(Li) 
detector  which  has  a  typical  resolution  of  3.5  KeV  FWHM  at  1.33  MeV. 


■ 
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Gamma  rays  from  137Cs  and  22Na  sources  were  recorded  with  the  53Cr 
spectra  ensuring  an  accurate  calibration  of  53Cr  gamma  rays.  Above 
proton  bombarding  energy  of  2.0  MeV,  all  spectra  were  dominated  by 
gamma  rays  from  the  competing  (d,ny)  reaction.  However,  an  optimum 
yield  for  gamma  rays  from  the  first  two  excited  states  of  53Cr  was 
obtained  at  bombarding  energies  of  3.2  and  3.5  MeV. 
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4.  EXPERIMENTAL  RESULTS 


A  typical  charged  particle  spectrum  taken  in  coincidence  with 
gamma  rays  is  shown  in  fig.  3.  This  indicates  the  relative  intensities 
of  the  excited  states  of  53Cr  and  also  shows  that  two  states  of  29Si 
were  observed  via  the  (d,p)  reaction  on  silicon  contaminant  in  the 
target.  The  energies  of  the  first  two  states  in  53Cr  are  taken  from 
the  Ge(Li)  measurements  reported  in  section  2,  but  the  energies  listed 
by  Rao  et  al.  (Ra  68)  are  used  for  all  the  other  excited  states.  Many 
of  the  excited  states  were  too  weakly  excited  to  obtain  angular  distri¬ 
butions  and  hence  only  their  decay  modes  could  be  established.  Fig.  4 
summarises  all  the  branching  ratios  determined  in  the  present  work. 

For  the  multiplets  near  2.7  and  3.2  MeV,  only  composite  branching  ratios 
are  given  in  Table  1. 

4.1.  The  Levels  Below  2.0  MeV 

The  isotropic  angular  distribution  obtained  for  the  ground  state 
decay  of  the  0.564  MeV  level  is  consistent  with  the  previous  assignment 
of  1/2-  based  on  the  postulated  J-dependent  effects  in  (d,p)  stripping 
(An  64) .  The  energy  of  the  gamma  ray  to  the  ground  state  was 
accurately  determined  to  be  (564.1  +0.2)  keV  from  its  observation  in 
the  53Cr (p ,p  ’y) 53Cr  reaction  at  3.2  and  3.5  MeV  bombarding  energies. 

The  observed  Doppler  shift  obtained  from  measurements  taken  at  0°  and 
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Figure  3  Charged  particle  spectrum  in  coincidence  with  gamma  rays  for  the 
52Cr (d ,py) 5 3Cr  reaction  at  a  bombarding  energy  of  4.11  MeV .  The 
protons  were  detected  near  172°.  The  associated  random 
spectrum  has  been  subtracted.  Proton  peaks  are  identified  by 
the  energies  of  the  states  in  53Cr  to  which  they  correspond. 

The  peak  marked  by  an  asterisk  (*)  has  not  been  identified. 


S4uno3  jo  jaqujnfsj 


Proton  Energy  (Mev) 
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Figure  4  A  summary  of  branching  ratios  and  spin  assignments  of  some 
of  the  excited  states  of  53Cr  below  4.3  MeV.  The  branching 
ratios  have  been  measured  in  the  present  experiment.  The  spin 
assignments  include  our  measurements  as  well  as  all  previously 
available  information. 
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90°  with  respect  to  the  beam  axis  was  less  than  0.2  keV  compared  to 
the  full  Doppler  shift  of  1.5  keV.  This  corresponds  to  a  mean  life¬ 
time  x  >_  0.5  ps  for  the  0.564  MeV  state. 

The  1.006  MeV  level  decays  entirely  to  ground.  It  is  fed  by  an 
=  3  transfer  in  the  (d,p)  reaction  with  a  spectroscopic  strength  of 
S  ~  2.5  (Bo  65).  From  the  angular  distribution  of  the  1006  keV  gamma 
ray,  the  7/2“  assignment  is  just  rejected  on  the  basis  of  the  0.1% 
confidence  criterion,  confirming  the  previous  5/2"  spin  assignment. 

The  extracted  mixing  ratio  for  this  transition  is  -0.27^'^  . 

From  the  5 3Cr (p , p ' y) 5 3Cr  experiment,  the  energy  of  this  transition  was 
measured  to  be  (1006.0  +0.2)  keV.  With  the  gamma  detector  at  0°  and 
90°  with  respect  to  the  beam  axis,  the  Doppler  shift  was  less  than 
0.2  keV  compared  to  a  full  shift  of  3  keV.  This  corresponds  to  a 
lifetime  x  >  1  ps. 

The  1.285  MeV  state  decays  predominantly  to  the  ground  state.  We 
observed  only  a  9%  branch  to  the  1.006  MeV,  considerably  smaller  than 
the  previous  measurements  (Ro  65,  Ra  68).  The  branching  ratio  was 
established  in  the  present  work  by  direct  observation  of  the  278  keV 
gamma  ray,  while  in  previous  measurements  it  was  inferred  from  the 
presence  of  the  1006  keV  gamma  ray.  The  angular  distribution  of  the 
1285  keV  gamma  ray  is  consistent  with  a  pure  E2  transition  from  a  7/2 
state  confirming  the  previous  assignment  from  observation  of  &  =  3 

momentum  transfer  in  both  the  stripping  (d,p)  (Ra  68),  and  the  pick-up 


(p,d)  reactions  (Wh  67). 


3  d 
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The  gamma -ray  spectrum  obtained  for  the  1.537  MeV  state  is  shown 

in  fig.  5.  This  indicates  a  previously  unobserved  branch  of  11%  to 

the  ground  state.  The  angular  distribution  of  the  529  keV  gamma  ray  to 

the  1.006  MeV  state  was  strongly  anisotropic  at  both  bombarding  energies 

and  indicated  a  probable  spin  of  7/2.  This  assignment  is  in  agreement 

with  the  work  of  Whitten  (Wh  67)  who  reported  that  this  state  is  the 

strongest  excited  in  the  pick-up  (p,d)  reaction  with  an  £  =3  transfer, 

n 

indicating  the  state  has  a  predominant  iq  hole  configuration. 

T 

The  1.971  MeV  state  decays  85%  to  the  (3/2“)  ground  state  and  15% 
to  the  1.285  MeV  (7/2-)  state.  This  latter  branch  was  not  observed  in 
previous  work  (Ro  65). 

4.2.  The  Levels  Between  2.0  and  3.6  MeV 

The  2.227  MeV  state  decays  only  to  the  1.537  MeV  (7/2“)  state. 

The  angular  distribution  of  the  690  keV  gamma  ray  shown  in  fig.  6a 
can  be  equally  well  fitted  by  spin  assignments  5/2,  7/2  and  9/2  as 
the  X  plots  in  fig  6b  and  6c  indicate.  The  state  is  weakly  excited 

in  the  (d,p)  reaction  and  exhibits  no  characteristic  stripping  pattern 

2 

(Ra  68).  Hence,  the  X  analysis  was  repeated  with  limits  of  0.3  and 
3.0  placed  on  the  ratio  of  the  population  parameters  as  discussed 
in  section  1.4.  In  the  region  of  best  fit,  this  ratio  was  previously 
either  0  or o*-.,  and  in  the  cases  of  5/2  and  9/2,  the  introduction  of 
the  limits  has  only  a  small  effect  on  the  X2  fit  as  shown  in  fig.  6b. 
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Figure  5 


Gamma  rays  in  coincidence  with  protons  populating  the  1.537 
MeV  level  in  53Cr.  The  lower  spectrum  corresponds  to  the 
simultaneous  random  spectrum. 


Number  of  Counts 


Energy  (MeV) 
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Figure  6 


2 

The  2.227  MeV  state.  The  angular  distribution  (a)  and  X  plots 
(b  and  c)  for  the  690  keV  transition  to  the  1.537  MeV (7/2  )  level 
in  53Cr  are  shown.  The  X2  plots  marked  by  an  asterisk  (*) 
correspond  to  limits  of  0.3  to  3.0  set  on  the  ratio  of  population 
parameters  (see  section  1.4).  The  fits  shown  with  the  angular 
distributions  correspond  to  the  minimum  y2  with  limits  on  the 
population  parameter  ratio. 
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However,  the  new  limits  did  enable  the  spin  assignment  1/2  to  be 
completely  excluded  (fig.  6c). 

The  angular  distribution  of  the  ground  state  transition  from 
the  2.324  MeV  state  was  definitely  anisotropic  with  an  a2  =  0.15  + 

0.05.  This  confirmed  the  previous  tentative  spin  assignment  of 
3/2“  on  the  basis  of  the  J-dependent  effect  in  the  strong  £  =  1 

neutron  stripping  angular  distribution  (An  64) . 

There  is  a  triplet  of  levels  at  ~  2.67  MeV  identified  res¬ 
pectively  as:  2.664  MeV  (£  =  3),  2.676  MeV  (£  =  1)  and  2.715  MeV 

n  n 

(£  =  1) .  The  gamma  decay  in  coincidence  with  the  3  proton  groups 

is  shown  in  fig.  7a.  The  angular  distributions  of  the  2.13  MeV 
and  2.67  MeV  gamma  rays  were  measured.  The  latter  was  isotropic 
but  the  former  enables  a  unique  spin  assignment  of  3/2,  provided 
the  2.13  MeV  gamma  ray  comes  predominantly  from  a  single  level, 
both  1/2  and  5/2  being  rejected  on  the  basis  of  the  0.1%  confidence 
test . 

The  gamma  spectrum  obtained  in  coincidence  with  the  proton  group 
corresponding  to  the  3.190  MeV  is  shown  in  fig.  7b.  The  proton  spectrum 
in  fig.  3  indicates  that  at  most,  25%  of  the  proton  peak  feeding  the 
state  at  3.190  MeV  could  be  attributed  to  the  proton  group  feeding  the 
3.146  MeV  level.  Our  branching  ratio  measurements  are  in  essential 
agreement  with  those  previously  measured  by  Rollefson  e_t  a_l.  (Ro  65) 
except  for  a  weak  (8%)  transition  to  the  1.285  MeV  (7/2-)  level.  The 
angular  distribution  of  the  predominant  decay  mode  of  the  3.190  MeV 
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Figure  7 


)  Gamma  rays  in  coincidence  with  protons  populating  the  triplet 
in  53Cr  at  2.66.4  MeV,  2.676  MeV  and  2.715  MeV. 

)  Gamma  rays  in  coincidence  with  protons  populating  the  3.146 
MeV  and  3.190  MeV  levels  in  53Cr. 


Number  of  Counts  Number  of  Counts 
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Table  I  A  summary  of  Legendre  polynomial  coefficients  fitted  to  the  angular 
distributions  of  transitions  observed  in  the  52Cr (d ,py) 5 3Cr  reaction 


at  a 

bombarding 

energy  of  4. 

11  MeV. 

Initial  Level 
(MeV) 

Final  Level 
(MeV) 

a2 

a 

4 

Multipole  Mixing  Ratios 

0.564 

1" 

2 

0. 

3“ 

2 

-0.12+0.09 

0.07+0.09 

undetermined 

1.006 

5“ 

2 

0. 

3“ 

2 

0.15+0.06 

-0.05+0.08 

1.285 

7“ 

2 

0. 

3" 

2 

0 . 55+0 . 11* 

-0.56+0.18* 

0-  m 

1.537 

7“ 

1.006 

5“ 

0.56+0.10, 

0.19+0.10, 

0  27+0  * 33 

U‘Z/-0.28 

2 

2 

-0.38+0.15* 

0.20+0.195 

2.227 

(!■!) 

1.537 

7~ 

2 

0.58+0.15 

-0.11+0.15 

0.23+°-2E  ;  1.8+°'2  (I) 

-0.15  *  -0.6  \2l 

-6 . 0<6<0 . 2  ~  ) 

—  —  12/ 

2.324 

3“ 

2 

0. 

3“ 

2 

0 . 15+0 . 05 

0.09+0.06 

-5 . 7<_6<0 .09  ;  0.47^6^5.7 

(2.70) 

3" 

2 

0.564 

1“ 

2 

-0.30+0.11 

0.17+0.13 

undetermined 

3.190 

3" 

2 

0.564 

1" 

2 

0.58+0.19 

0.21+0.20 

undetermined 

3.271 

5 

2 

0. 

3“ 

2 

-0.83+0.11 

0.06+0.11 

0.06<6<1.8 

3.625 

3" 

2 

0. 

3~ 

2 

-0.07+0.05 

0.08+0.06 

undetermined 

3.715 

9+ 

2 

1.285 

7“ 

2 

-0.36+0.03 

0.08+0.04 

.  „„+0.09  ,  -+0.8 

°-00-O.O4  5  3  *  7-0 .6 

3.971 

(H) 

0.564 

1" 

2 

0.03+0.08 

-0 . 04+0 . 09 

undetermined 

4.136  ^ 

3+  5+\ 
2  ’2  ) 

2.324 

3~ 

2 

-0.49+0.06 

-0.04+0.07 

-5 . 4<6<0 . 15 ; 0 . 5 5 <6 <4 . 0 

-0 . 05<_6<0 . 7  ;  1 . 4<6<_3 . 8 

4.228  | 

,2  2  / 

1.006 

5" 

2 

0.17+0.04 

-0.03+0.05 

-10  .£6<0. 6 

* 

Results  obtained  at 

Ed  ■ 

■  4.85  MeV. 
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state  to  the  0.564  MeV  (1/2-)  state  will  be  least  affected  by  the  presence 
of  the  3.146  MeV  level.  Hence,  the  angular  distribution  of  the  resulting 
2626  keV  gamma  ray  enables  a  unique  3/2  spin  assignment  to  be  made  to 
the  3.190  MeV  state;  all  other  possible  spin  assignments  were  rejected 
on  the  basis  of  the  0.1%  confidence  limit  criterion.  This  assignment 
is  in  agreement  with  Rao  et_  aM.  (Ra  68)  who  reported  a  possible  5,  =  1 
momentum  transfer  to  this  state  from  the  stripping  (d,p)  reaction. 

The  3.271  MeV  state  decays  20%  to  the  (7/2”)  1.537  MeV  state,  23% 
to  the  (5/2~)  1.006  MeV  state  and  57%  to  the  3/2“  ground  state.  The 
angular  distribution  of  the  3271  keV  gamma  ray  shown  in  fig  8a,  is 
strongly  anisotropic  with  a2  =  -0.83  +  0.11.  The  X  plots  in  fig.  8b 

indicate  that  spin  assignment  3/2  and  5/2  both  give  acceptable  fits  to 

o 

the  data.  The  X  analysis  was  repeated  with  limits  of  0.3  and  3.0  on 
the  ratio  of  population  parameters.  The  new  plots  shown  in  fig.  8b 
enabled  a  definite  spin  assignment  of  5/2  to  be  made.  A  tentative 

=  2  transfer  was  assigned  by  Rao  e£  al .  (Ra  68)  to  this  weak  stripping 
state  inferring  the  3.271  MeV  level  has  positive  parity. 

4.3.  Levels  Above  3.6  MeV 

The  3.625  MeV  state  which  is  strongly  excited  in  the  L’2Cr(d,p) 

reaction  with  an  £  =1  momentum  transfer  decays  only  to  ground  state. 

n 

The  distribution  of  the  3625  keV  gamma  ray  shows  a  very  weak  anisotropy 
a2  =  -0.07  +0.05  which  is  consistent  with  the  previous  tentative  3/2” 
assignment  (An  64,  Ba  65).  Recent  polarization  measurements  (Yu  68) 


35 


2 

Figure  8  The  3.271  MeV  state.  The  angular  distribution  (a)  and  X  plots 

n 

(b)  for  the  3271  keV  ground  state  transition  are  shown.  The  X 
plots  marked  by  an  asterisk  correspond  to  limits  of  0.3  to  3.0 
set  on  the  ratio  of  population  parameters  (see  section  1.4«)»  The 
fits  shown  with  the  angular  distribution  correspond  to  the  mini- 
mum  X  with  limits  on  the  population  parameter  ratio. 
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indicate  that  the  3.625  MeV  level  is  a  1/2“  state.  This  would  mean  that 
our  calculated  error  in  the  a2  coefficient  is  slightly  underestimated. 

The  3.715  MeV  level  is  a  strong  deuteron  stripping  state  with  a 
neutron  momentum  transfer  &  =  4  (An  64,  Ra  68,  Bo  65),  therefore  suggest¬ 

ing  it  can  be  assigned  either  7/2+  or  9/2+.  It  decays  entirely  to  the 
1.285  MeV  (7/2”)  level  with  a  2430  keV  gamma  ray  transition.  The  angular 
distributions  of  the  2430  keV  and  1285  keV  gamma  rays,  shown  in  fig.  9, 
were  analysed  simultaneously.  The  X  plots,  in  fig  9,  indicate  that  the 
7/2  assignment  is  ruled  out  by  the  0.1%  confidence  limit  test  leaving 
9/2+  as  the  spin  and  parity  assignment  of  the  3.715  MeV  level.  Of  the 
two  possible  values  of  the  multipole  mixing  ratio  for  the  9/2"*"  to  7/2 
transition,  the  solution  which  corresponds  to  pure  electric  dipole  is 
more  likely  than  the  solution  corresponding  to  strong  M2-E1  mixing. 

The  3.985  MeV  level  decays  predominantly  to  the  0.564  MeV (1/2”) 
level  with  a  75%  branching  and  has  a  25%  branching  ratio  to  ground. 

o 

The  X  analysis  of  the  isotropic  angular  distribution  of  the  3407  keV 
gamma  ray  enables  the  5/2  spin  assignment  to  be  rejected  leaving  only 
1/2  and  3/2  as  possible  spin  assignments  for  the  3.971  MeV  state.  There 
exists  an  inconsistency  in  the  literature  as  both  the  =  2  momentum 
transfer  in  (d,p)  stripping  reaction  and  the  £  =  1  momentum  transfer 

in  the  (3He,a)  pick-up  reaction  have  been  reported  for  this  level  (Ra  68, 
Bo  65,  DD  69) . 

The  4.136  MeV  level  decays  46%  to  the  2.324  MeV (3/2  )  level,  37% 
to  the  1.006  MeV  (5/2”)  level  and  17%  to  the  (3/2”)  ground  state.  This 
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Figure  9 


The  3.715  MeV  state.  The  angular  distributions  of  the  2430 

o 

keV  and  1285  keV  gamma  rays  in  cascade  are  shown.  The  X 
plots  correspond  to  the  simultaneous  analysis  of  the  angular 
distributions . 
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level  is  excited  strongly  in  the  (d,p)  stripping  reaction  with  an  £  =2 

n 

2 

momentum  transfer.  However,  the  X  analysis  of  the  angular  distributions 
of  the  3128  keV  transition  to  the  1.006  MeV  (5/2“)  level  and  the  1812  keV 
transition  to  the  2.324  MeV  (3/2~)  level  did  not  distinguish  between  spin 
assignments  of  3/2  and  5/2  for  the  4.136  MeV  level  but  did  enable  limits 
to  be  put  on  the  possible  multipole  mixing  ratios  for  both  spin  possibil¬ 
ities  . 

The  4.228  MeV  state  is  the  next  state  which  is  strongly  excited  in 
the  (d,p)  stripping  reaction  with  an  =  2  momentum  transfer.  We  found 
no  evidence  of  a  nearby  excited  state  at  4.204  MeV  which  Rao  et  al . 

(Ra  68)  report  to  be  very  weakly  excited  in  the  stripping  (d,p)  reaction. 
We  measured  a  59%  branching  to  the  1.006  MeV  (5/2“)  state,  30%  branching 
to  the  member  of  the  triplet  at  2.67  MeV  that  decays  predominantly  to 
ground  and  an  11%  branching  to  the  1.285  MeV  (7/2“)  state.  The  analysis 
of  the  3220  keV  transition  to  the  1.006  MeV  level  was  inconclusive  for 
spin  assignment  but  enabled  limits  to  be  put  on  the  possible  multipole 


mixing  ratios. 
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CHAPTER  III 


THE  55Mn(p,ny)55Fe  REACTION 

In  Chapter  II,  a  general  expression  for  angular  distribution  of  gamma 
rays  resulting  from  the  decay  of  aligned  states,  has  been  introduced.  In 
the  same  chapter,  a  specific  case  of  alignment,  namely  the  triple  angular 
correlation,  was  considered.  This  chapter  deals  with  another  way  of  pre¬ 
serving  nuclear  alignment,  namely,  by  ensuring  that  outgoing  particles 
are  predominantly  s-wave.  In  the  following  sections,  the  general  theory 
describing  the  statistical  model  formalism  will  be  outlined.  Specific 
application  of  that  theory  to  the  case  of  the  55Mn(p,ny) 55Fe  will  then  be 
presented . 


'Jia 
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1.  THE  COMPOUND  NUCLEAR  STATISTICAL  THEORY  OF 
GAMMA-RAY  ANGULAR  DISTRIBUTIONS 


1.1.  Introduction 


As  shown  in  Chapter  II,  the  angular  distribution  of  gamma  rays 
resulting  from  the  decay  of  aligned  states  can  be  written  as: 


W(0 )  = 


I  VV  Vcos6) 

even  k  K  1  R 


Rk(LLJiJ2)+26Rk(LL,JiJ2)+62Rk(L,L,JiJ2) 

1  +  52 


(IH.l) 


where  all  the  variables  have  been  previously  introduced. 

This  expression  is  the  product  of  three  terms:  B  (J^)^  which  contains 

all  the  information  on  the  alignment  of  the  initial  state  J^,  the  angular 

part,P  (cos0),  and  the  remaining  term  (hereafter  labelled  A  (LL'J  J  6)) 

K  12 

which  describes  the  interaction  between  the  nuclear  states  and  the  electro¬ 
magnetic  field. 

The  problem  now  is  twofold:  first,  to  determine  the  conditions  for 
alignment  in  a  reaction  such  as  55Mn(p,ny);  second,  to  evaluate  the  term 
in  equation  (III.l)  which  reflects  the  conditions  of  alignment,  namely, 

B^.  In  order  to  achieve  the  latter  goal,  the  population  parameters  (and 
consequently  the  B  ’ s)  will  be  evaluated  after  a  model  describing  the 

K. 


Bk(Jl)  =  I  P(mj )  /2Ji+ 1 


( J^m1 J  ^ -m^ | kO) 


' 
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reaction  mechanism  has  been  introduced. 

1.2.  Conditions  for  Alignment  in  the  A(a,b)B  Experiment 

Assume  the  reaction  A(a,b)B,c  takes  place  in  two  steps:  first, 
formation,  in  the  compound  nucleus,  of  a  state  Y  ,  |  J  m  >,  with  definite 
angular  momentum  and  parity,  then,  decay  by  emission  of  a  nucleon  b 
(fig.  10). 


J2 


B 


J 


3 


Figure  10 

In  the  first  step,  the  beam  direction  is  chosen  as  the  quantization 
axis  (z-axis) .  It  is  moreover  assumed  that  neither  the  beam  of  incident 
particles  nor  the  target  nuclei  are  polarized.  These  conditions  result 

J^m^CRo  67  and  Chapter  II). 


in  alignment  of  the  state 
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In  the  second  step,  a  nucleon  ’b’  is  emitted  but  not  observed.  The 
projection  on  the  z-axis  of  the  orbital  angular  momentum  of  this 
nucleon  is  no  longer  zero  as  it  can  be  emitted  in  any  direction.  One 
then  generally  expects  the  state  B  (|J2m2>)  to  be  de-aligned  as  the 
result  of  the  contributions  of  the  projection  of  on  the  z-axis.  How¬ 
ever,  if  the  reaction  proceeds  near  threshold,  then  s-wave  particles  will 
predominate  and  de-alignment  should  be  minimal. 

As  an  example,  consider  the  case  described  by  fig.  10.  The  incoming 

particle,  a,  has  an  orbital  angular  momentum  and  an  intrinsic  spin  s1 . 

The  outgoing  particle  b ,  has  an  orbital  angular  momentum  &2  and  an 

intrinsic  spin  s2.  J's  are  used  to  define  the  total  angular  momenta  of 

the  different  levels.  Since  the  beam  direction  is  taken  as  the  z-axis, 

the  substates  in  Y  are  populated  up  to  |m1 =  Jq+s^  If  only  s-wave 

particles  are  emitted  in  the  decay  of  Y  then  lm2lmax  =  (^0*  s  l  )"*"  s2  f°r  the 

residual  nuclear  state,  B  .  It  then  follows  that,  in  a  (p,ny)  reaction, 

the  largest  magnetic  substate  to  be  populated  in  the  residual  nuclear 

state  |j_m„>  is  |m~|  =  (JL  +  1) .  In  this  experiment,  the  target  nucleus 

55Mn  has  a  spin  Jn  =  5/2,  therefore  I m0 I  =  7/2. 

r  0  1  z 1  max 

Experimentally,  "just  above  threshold"  may  mean  as  much  as  30  to 
50  keV,  or  more, above  the  threshold  value.  Therefore,  contributions  from 
higher  order  partial  waves  may  be  present  but,  the  effect  on  the  nuclear 
alignment  of  the  residual  state  is  generally  not  very  important. 
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1.3.  Gamma  Angular  Distribution  and  the  Statistical  Model 

Consider  first  the  compound  nucleus  formation.  In  the  channel  spin 
representation,  the  population  parameter  P(m^)  is  proportional  to-.  (Li  61) 

P (m1 )  ~  (s1£1m10|  J^)2  (III. 2) 

If  more  than  one  channel  spin  contributes  to  the  formation  of  J^, 
the  different  contributions  are  added  incoherently: 


PG^) 


i  ( 

si 


si u,  m o 


si  Vi 


(HI. 3) 


If  the  reaction  taken  into  consideration  was  simply  a  resonant 
capture,  with  only  one  partial  wave  contributing,  one  could,  at  this 
stage,  introduce  in  the  coefficient  B^Jj)  the  value  of  P(m1)  defined  in 
equation  (III. 3)  to  get:  (Ro  67) 


B 


k(Jx)  =  I  (~)S1  Jl(2£i+1)  v/2^+1  U^OOlkO)  w(Vl  Vi;ksi)Ts  £  j 

s-^  111 

(III. 4) 


Putting  this  value  of  B^  into  the  expression  W(0),  (equation  III.l), 
would  then  give  the  equation  of  the  distribution  of  gamma-rays  following 
resonant  capture. 


The  problem  here  is  more  complex  as  there  are  other  terms  to  be 
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introduced  which  take  into  account  the  formation  and  decay  of  the  inter¬ 
mediate  state  |  '^2m2>  *  Basically,  one  can  expect  another  term  similar  to 
equation  (III . 4)  which  takes  into  account  the  alignment  of  the  state  Bx. 

As  a  result  of  the  different  couplings  involved,  there  will  be  additional 
Clebsch-Gordan  and  Racah  coefficients.  Following  Sheldon  (Sh  63),  the 
general  expression  giving  W(9)  schematically  reads  as: 


W(0)~  l  transition  parameter  for  the  formation  of  the  compound  nucleus 

I  T 


S1£1J1 


x 


transition  parameter  for  the  formation  of  the  intermediate 


state  J 


's  i  J 

2  2  2 


^  Vr j' 

S2&2  222 


> 


x  A  (LL'J.J.6)  x  P.  (cos0 ) 

y  1  2  k 


(III. 5) 


In  the  second  term,  the  intensity  factor,  {...},  reflects  the  fact 
that,  once  the  compound  nucleus  has  been  formed,  there  are  many  channels 
open.  It  is  a  measure  of  the  probability  of  decay  by  a  particular  channel. 

At  this  stage,  the  only  unknowns  required  to  evaluate  W(0)  are  the 
transition  parameters,!.  A  knowledge  of  the  reaction  mechanism  is  necessary 
to  calculate  these  factors.  The  problem  is  quite  difficult  but,  under 
some  assumptions,  one  can  get  a  good  estimate  of  the  T’s.  The  statistical 
model  has  been  extensively  described  in  the  past  few  years  (e.g.  see  Sh  66, 
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Vo  69,  Pi  69)  and  will  be  outlined  in  the  following  section. 

1.3.1  The  compound  nuclear  statistical  model 

The  model  assumes  that  the  reaction  proceeds  dominantly  through  the 
compound  nucleus  formation  (reactions  of  the  type  (p,ny)  and  (a,ny) 
satisfy  this  condition)  and  that  intermediate  states  have  definite 
parity . 

If  the  reaction  proceeds  via  a  large  number  of  overlapping  levels, 
it  can  be  shown  (Er  63)  that  the  fluctuation  in  the  cross  section  caused 
by  interference  between  these  levels  will  average  out  to  zero.  In 
(p,ny)  and  (a,ny)  reactions  on  medium  weight  nuclei,  the  density  of  levels 
in  the  compound  nucleus  is  large.  Experimentally,  the  use  of  thick  targets 
and  a  broad  beam  energy  resolution  contribute  to  the  excitation  of  a  large 
number  of  overlapping  levels.  This  last  assumption  is  the  basis  of  the 
statistical  model.  As  a  consequence,  it  is  then  possible  to  relate  the 
mechanism  of  the  reaction  to  energy  averaged  coefficients  called  transmis¬ 
sion  functions  or  transmission  coefficients.  These  coefficients  can  be 
calculated  using  codes  based  on  the  formalism  introduced  by  Hauser  and 
Feshbach  (Ha  52)  in  their  theory  describing  the  compound  nucleus  formation 
and  decay.  They  are  directly  related  to  the  T  introduced  earlier.  The 

S  AJ  J 

difference  lies  in  the  respective  momentum  space  representations.  The 
transmission  coefficients  are  usually  calculated  in  the  j-j  representation 
whereas  those  used  earlier  are  given  in  the  Jl-s  representation.  The 
relationship  between  the  two  definitions  is  given  by: 


, 

sd:t 
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t£sJ  -  l  (2j+l)  (2a +1)  W2(£IJs;ja)  T 
a 

where  W  is  a  Racah  coefficient. 

For  the  sake  of  completeness,  it  must  be  pointed  out  that  inter¬ 
ferences  between  various  partial  waves  for  both  incoming  and  outgoing 
waves  should  also  be  considered.  These  interference  effects  are  shown 
to  cancel  out  (Vo  69,  Pi  69). 

In  conclusion,  the  use  of  the  statistical  model  helps  to  describe 
the  mechanism  of  the  reaction.  It  is  assumed  that  the  reaction  takes 
place  in  two  separate  steps:  formation  of  a  compound  nucleus,  then 
decay  through  many  different  channels.  If  all  the  conditions  are  met 
so  that  no  resonance  can  occur,  it  is  possible  to  calculate  the  trans¬ 
mission  coefficients  which  in  turn,  yield  an  estimate  of  the  population 
parameters  of  the  intermediate  state.  The  complete  angular  distribution 
formula  can  then  be  calculated.  It  is  given  by:  (Pi  69) 

W(0)  =  I  (-)J0+J3~VS2(2£1+l)(2j1+l)(2s1+ir1(2J0+l)“1(2J1+l)2' 

kjlJ2 

1 

X  (2J2+1)^"  (AjfcjOOlkO)  w(^1Aij1j1;ks1)  W(JX J1 j l j 1 ;kJQ ) 

Rk(LLJ1J2)+26Rk(LLT  J1J2)+62Rk(L’L’ 1^2) 
x  W(J^ Jj J2J2; kj 2)  x  1  x  2+52 


x  P^  (cos0) 


(III. 6) 
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in  which  t  is  called  the  penetrability  term  and  is  given  by: 


x  = 


T«  i  (Ei>  h  j  (e2> 

«i3l  *232  2 


2JL 2. 

T£^<E> 


The  sum  extends  over  all  open  channels  (including  elastic  scattering) 


■k 


by  which  the  particular  compound  nuclear  state  Y  (J^  can  decay 
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2 .  EXPERIMENTAL  PROCEDURE 


2.1.  Apparatus 

In  this  experiment,  thick  targets  of  natural  Manganese  (100%  55Mn) 
were  bombarded  by  protons.  The  beam  was  collimated  to  a  diameter  of 
1.5  mm.  The  targets  were  prepared  by  mixing  a  small  amount  of  manganese 
powder  into  a  glue  formed  by  dissolving  polyurethane  into  benzene  and 
smearing  the  mixture  onto  a  0.112  mm  thick  tantalum  backing.  The  targets 
were  then  fixed  on  a  holder  constituted  by  a  2  mm  thick  aluminum  plate 
placed  at  30°  with  respect  to  the  beam  direction.  This  plate  also  served 
as  charge  collector. 

Gamma  rays  were  observed  using  a  15  cc  Ge(Li)  detector  which  could 
be  positioned  at  angles  ranging  between  0°  and  90°  to  the  beam  direction. 
The  resolution  of  the  Ge(Li)  detector  was  about  3.5  keV  f.w.h.m.  for 
the  1.332  MeV  peak  of  60Co. 

After  amplif ication,  the  gamma-ray  pulses  were  sent  into  an  analog- 
to-digital  convertor  (ADC) ,  itself  connected  to  the  on-line  SDS  920 
computer  for  live  display.  Self  gating  was  used  to  reduce  the  counting 
rate,  the  cut-off  being  set  at  about  300  keV. 

2.2.  Acquisition  and  Analysis  of  Data 

Following  a  preliminary  study  of  the  gamma-ray  yield  for  bombarding 
energies  ranging  between  3.60  and  5.00  MeV,  the  decay  scheme  established 
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by  Pilt  (Pi  69)  was  slightly  modified  and  extended  up  to  the  3  MeV  region. 
Two  bombarding  energies  were  chosen  for  the  study  of  angular  distribution, 
bearing  in  mind  that  some  de-alignment  is  allowed  when  the  spin  value  of 
the  level  is  large  (>_  7/2)  whereas  for  low  spin  states  it  is  necessary 

to  excite  the  level  as  close  as  possible  to  the  threshold,  yield  permit¬ 
ting. 

Angular  distributions  were  taken  at  five  angles  of  0,  30,  45,  60  and 
90  degrees.  Measurements  were  repeated  at  least  once.  The  beam  current 
was  kept  at  about  100  nA  so  that  the  counting  rate  remained  at  a  reason¬ 
able  value. 

There  was  no  problem  caused  by  contamination  in  this  experiment  as 
the  usual  contaminants,  oxygen,  carbon  and  nitrogen  have  a  very  low  Q- 
value  for  all  reactions  induced  by  protons.  The  gamma  rays  observed  in 
this  reaction,  which  were  not  produced  by  55Fe,  originated  mainly  from 
the  decay  of  levels  produced  through  the  other  open  channels  in  the 
55Mn  +  p  reaction:  (p,p'y),  (p,y)  and  (p,ay).  However,  only  the  first 
two  or  three  levels  of  the  corresponding  residual  nuclei  were  excited. 

Two  gamma  rays  resulting  from  the  27Al(p,p'y)  reaction  were  also  observed. 
The  gamma  rays  at  598  and  693  keV  resulting  from  the  neutron  Inelastic 
scattering  in  Germanium  were  also  present.  A  typical  gamma-ray  spectrum 
is  shown  on  fig.  11. 

The  gamma-ray  decay  scheme  of  55Fe  was  determined  by  observing  new 
transitions  when  bombarding  energy  was  increased  and  by  comparison  with  all 
possible  transitions  from  the  known  levels  of  55Fe  being  excited. 
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Figure  11  A  typical  gamma  ray  spectrum  taken  at  a  bombarding  energy  of 
4.35  MeV  with  the  15  cc  Ge(Li)  detector. 
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Branching  ratios  were  determined  from  the  sum  of  different  runs  weighted 
by  sin0  for  both  sets  of  data  corresponding  to  proton  bombarding  energies 
of  4.20  and  4.35  MeV  respectively.  The  absolute  intensities  were 
obtained  after  correction  for  absorption  and  efficiency.  The  relative 
efficiency  of  the  15  cc  Ge('Li)  detector  was  measured  using  a  56Co  source 
(fig.  12). 

The  strong  410  keV  gamma  ray  corresponding  to  the  decay  of  the  first 
excited  state  (1/2  )  of  55Fe  to  ground  state  was  used  as  internal  monitor 
as  it  has  an  isotropic  distribution. 

The  intensities  of  the  gamma  rays  were  determined  by  summing  over 
the  full  energy  peaks  and  subtracting  the  background  directly  from 
computer  display.  Some  intensities  were  recalculated  using  a  peak 
fitting  program  written  by  Tepel  (Te  66) .  The  results  obtained  from 
both  methods  generally  agreed  within  statistical  error. 

The  experimental  angular  distributions  were  fitted  to  the  expansion: 

W(0)  =  aQ  [1  +  a  2P2  ( cos 0 )  +  a^P^ (cos0 ) ] 

using  a  least  squares  procedure  with  the  Legendre  polynomial  coefficients 
treated  as  free  parameters.  The  experimental  a2  and  a^  coefficients  were 
then  compared  to  the  theoretical  predictions.  The  latter  were  calculated 
for  different  spin  sequences  and  multipolarities  using  the  computer  code 
MANDY,  originally  written  by  Sheldon  (Sh  69).  The  transmission  coefficients 
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Figure  12  Relative  efficiency  curve  for  the  15  cc  Ge(Li)  detector.  The 
experimental  points  correspond  to  measurements  taken  using  a 
56Co  source. 
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needed  for  the  MANDY  calculations  were  obtained  from  the  Hauser-Feshbach 
program  written  by  Davison  (Da  69) ,  using  the  optical  model  parameters  of 
Rosen  (Ro  66)  for  the  proton  channels,  of  Perey  (Pe  63)  for  the  neutron 
channels  and  of  Davison  (Da  69)  for  the  a  channel. 

Spins  were  assigned  on  the  basis  of  the  0.1%  confidence  limit  of 

O 

the  X  fitting.  [This  method  has  been  described  in  Chapter  II.]  The  Rose 
and  Brink  (Ro  67)  sign  convention  for  the  mixing  ratio  6  was  adopted. 
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3.  EXPERIMENTAL  RESULTS 


A  typical  Ge(Li)  spectrum  taken  at  4.35  MeV  bombarding  energy  is 
shown  in  fig.  11.  The  gamma  ray  spectra  were  calibrated  with  22Na, 
80Co  and  208Tl(Th  C")  sources.  The  energy  values  obtained  agree 
essentially  with  the  more  accurate  energy  measurements  of  Robertson 
(Ro  70) ,  which  have  b  een  adopted  in  this  work. 

Figure  13  gives  a  decay  scheme  of  55Fe  with  branching  ratios  and 
spin  assignments  obtained  in  this  work  as  well  as  from  previous  work 
(Pi  69). 

Tables  II  -  IV  summarize  all  the  results  obtained  in  this  exper¬ 
iment:  branching  ratios,  Legendre  coefficients  of  the  angular  distri¬ 

butions,  deduced  spins  and  mixing  ratios. 

3.1.  Determination  of  Branching  Ratios 

The  branching  ratio  values  determined  in  this  work  are  given  in 
Table  II,  where  they  are  compared  with  the  measurements  of  Fischbeck 
et^  aM.  (Fi  66)  and  Pilt  (Pi  69). 

The  1666  keV  has  been  totally  assigned  to  the  decay  of  the  3.076 
MeV  level  reported  by  Sperduto  and  Buechner  (Sp  64)  to  the  1.410  MeV 
(7/2“)  level.  However,  a  gamma  ray  of  this  energy  could  also  corre¬ 
spond  to  the  decay  of  the  2.984  MeV  level  to  the  1.318  MeV  (7/2  ) 


level.  Evidence  of  the  excitation  of  the  2.984  MeV  level  was  seen  at 


iol 

' 
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Table  II 


A  summary  of  branching  ratios  of  the  excited  states  in  55Fe  obtained 
in  the  present  experiment  with  earlier  investigations 


Initial  State 
(MeV) 

Residual  State 
(MeV) 

This  Exp. 

Pilt  (Pi  69) 

Fischbeck  et  al 
(Fi  66) 

2.144 

0 

18  +  2 

17  +  2 

100 

0.410 

3  +  1 

0.930 

43  +  4 

58  +  4 

1.318 

36  +  5 

25  +  3 

2.211 

1.318 

(20) 

(20) 

not  observed 

1.410 

100 

100 

2.301 

0.930 

75  +  5 

80  +  10 

1.318 

15  +  5 

20  +  10 

not  observed 

1.410 

(10) 

(10) 

2.542 

1.318 

100 

not  observed 

not  observed 

2.578 

0 

84  +  2 

(100) 

100 

0.410 

7  +  2  * 

0.930 

6  +  2 

1.318 

3  +  2 

2.871 

0 

88  +  3 

(100) 

100 

0.930 

12  +  3 

2.938 

0 

55  +  5 

100 

55 

1.316 

45  +  4 

45 

2.984 

0.930 

<10 

not  observed 

1.408 

(100) 

3.027 

0. 

65  +  6 

not  observed 

0.412 

35  +  7 

1.408 

<1 

3.076 

1.408 

(100) 

not  observed 

■ 
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bombarding  energies  greater  than  4.25  MeV  with  the  appearance  on  the  gamma- 
ray  spectrum, of  the  1574  keV  gamma  ray  corresponding  to  the  transition  to 
the  1.410  MeV  (7/2")  level.  Unfortunately,  the  bombarding  energy  of 
4.25  MeV  corresponds  to  an  excitation  of  80  keV  above  the  threshold  of  the 
3.076  MeV.  On  this  basis,  is  was  not  possible  to  assign  unambiguously  the 
1666  keV  transition.  Arguments  in  favor  of  assigning  this  transition  to 
the  decay  of  the  3.076  MeV  state  are  the  following.  If  the  1666  keV 
transition  were  totally  assigned  to  the  2.984  MeV  level,  the  measurement 
of  the  intensities  of  the  1666  and  1574  keV  lines  would  indicate  that  the 
decay  is  50%  to  the  1.410  MeV  (7/2")  level  and  50%  to  the  1.318  MeV  (7/2") 
level.  The  two  7/2"  levels  have  configurations  which  are  essentially 

different  (see  Chapter  IV),  one  is  predominantly  a  fy  hole  state  (G1  66), 

2 

the  other  is  the  result  of  the  coupling  of  the  P3  neutron  to  the  5t+Fe  core. 

Y 

Consequently,  one  would  expect  the  2.984  MeV  level  to  decay  predominantly 
to  either  one  of  the  7/2”  states  but  not  to  both  with  equal  strength.  On 
the  other  hand,  at  4.35  MeV  bombarding  energy,  the  analysis  of  the  angular 
distribution  of  the  1574  keV  gamma  ray  yields  an  unambiguous  9/2  assignment 
to  the  2.984  MeV  level,  whereas  the  angular  distribution  of  1666  keV  gamma 
ray  favors  strongly  an  11/2  assignment  (although  9/2  cannot  be  entirely 
ruled  out  on  the  basis  of  the  0.1%  confidence  limit).  Therefore,  the  1666 
keV  gamma-ray  transition  was  entirely  assigned  to  the  3.076  MeV  state.  A 
coincidence  measurement  (e.g.  the  5''+Fe(d,py)  55Fe  reaction)  could  confirm 


this  assignment. 
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Figure  13  The  ^Fe  energy  level  diagram:  The  branching  ratios  for  levels 
below  2.144  MeV  were  determined  by  Pilt  (Pi  69).  For  levels 
>_  2.144  MeV#  they  were  determined  in  this  experiment.  The  £ 
values  are  taken  from  (Fu  63)  and  (Ma  64) . 


3076  - - - (100) 


3027  - 65  +  6—35+7 - <1 -  i 
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3.2.  The  2.211  MeV  Level 


This  level  has  been  assigned  J  =  9/2  by  Pilt  (Pi  69)  using  the 

55Mn(p,ny) 55Fe  reaction,  at  3.50  MeV  bombarding  energy  (210  keV  above 

threshold).  At  4.10  MeV  bombarding  energy  (810  keV  above  threshold), 

the  distribution  of  the  801  keV  gamma  ray  corresponding  to  the  decay 

of  the  2.211  MeV  (9/2")  level  to  the  1.410  MeV  (7/2")  level  is  still 

strongly  anisotropic,  (a2  =  -0.34  +  0.03)  as  shown  in  fig.  14a.  The 
2 

X  plots  shown  in  fig.  14b  indicate  that,  even  at  810  keV  above 
threshold,  the  2.211  MeV  level  can  be  unambiguously  assigned  J  =  9/2. 
The  possible  mixing  ratio  solutions  (0.21  or  2.15  )  are 

consistent  with  those  obtained  by  Pilt  (Pi  69)  (0.19^’^q  or  2.75 


3.3.  The  2.301  MeV  Level 

This  level  has  ^lso  been  assigned  J  =  9/2  by  Pilt  (Pi  69)  from 
the  study,  at  3.50  MeV  proton  bombarding  energy,  of  the  1371  keV  tran¬ 
sition  corresponding  to  the  decay  from  this  level  to  the  0.930  MeV 
(5/2)  level.  The  angular  distribution  of  the  1371  keV  transition  was 
measured  at  4.10  and  4.35  MeV  bombarding  energy  in  this  experiment. 

In  Table  III,  the  Legendre  polynomial  coefficients  obtained  in  these 
measurements  are  given  and  compared  with  Pilt's  measurements.  The 
resulting  spin  assignment  is  also  given. 
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Figure  14  The  2.211  MeV  level:  Angular  distribution  (a)  and  X2  plots  (b) 
for  the  801  keV  gamma-ray  transition  to  the  1.410  MeV  (7/2) 
level . 
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Table  III 


Bombarding  Energy  above  __  Assignment 

Energy  Threshold  a2  a4 

(MeV)  (MeV) 


3.50  (Pi  69) 

125 

0.39 

+ 

0.05 

0.03 

+ 

0.06 

9/2 

4.10 

725 

0.22 

+ 

0.02 

-0.01 

+ 

0.02 

9/2 

4.35 

975 

0.19 

+ 

0.03 

-0.00 

+ 

0.03 

9/2, 7/2 

2 

GO  The  X  plot  corresponding  to  J  =  7/2  has  a  minimum  of  2,4;  this  value 
corresponds  to  about  5%  confidence  limit. 

This  table,  as  well  as  the  results  obtained  for  the  2.211  MeV  level, 
indicate  that  in  this  experiment,  levels  with  spin  9/2  (and  consequently 
with  spin  >_9/2),  remain  strongly  aligned  when  excited  at  a  few  hundred 
keV  above  threshold. 

3.4.  The  2.542  MeV  Level 

The  2.542  MeV  level  was  found  to  decay  entirely  to  the  1.318  MeV 
(7/2-)  state  via  a  1224  keV  transition.  This  level  has  been  reported 
by  Bjerregaard  et  al.  (Bj  64),  Sperduto  and  Buechner  (Sp  64)  from  their 
study  of  the  54Fe (d ,p) 5 bFe  reaction.  Fischbeck  et^  al_.  (Fe  66)  in  the 
study  of  the  (3+  decay  of  55Co  have  found  no  evidence  that  this  level  is 
fed  in  the  (3+  decay.  They  have  concluded  that  this  level  cannot  be 
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Figure  15  The  2.542  MeV  level:  Angular  distributions  (a)  and  X2  plots  (b) 
for  the  1224  keV  gamma-ray  transition  to  the  1.318  MeV  (7/2”) 
level . 
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Legendre  coefficients,  corrected  for  solid  angle  effects,  fitted  to 
angular  distributions  of  gamma  rays  deexciting  levels  of  55Fe 
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assigned  5/2",  7/2”  or  9/2". 

The  angular  distribution  of  the  1224  keV  gamma  ray  at  4.10  MeV 
bombarding  energy  is  shown  in  fig.  15a.  It  is  strongly  anisotropic 
(a2  =  0.30  +  0.03).  The  X2  analysis  (fig.  15b)  yields  two  possible 
spin  assignments:  9/2  and  11/2,  the  former  being  less  likely  than 
the  latter.  The  9/2”  assignment  could  be  rejected  on  the  basis  of  the 
results  of  Fischbeck  e_t  _al.  Positive  parity  for  the  J  =  9/2  assignment 
is  not  favored  as  one  would  rather  expect  the  level  to  decay  via  a  pure 
electric  dipole  transition  in  this  case,  assumption  in  disagreement 
with  the  results  of  the  X  analysis.  Therefore,  the  =  11/2  assign¬ 
ment  is  more  likely  and  corresponds  to  a  pure  E2  transition  (6=0). 

3.5.  The  2.578  MeV  Level 

This  level  decays  primarily  to  the  ground  state  (84%)  .  In  this 
work,  three  weak  transitions,  not  previously  reported  (Fi  66,  Pi  69), 
have  been  observed.  They  correspond  to  decay  to  the  0.410  MeV  (1/2  ), 
0.930  MeV  (5/2")  and  1.318  MeV  (7/2")  levels. 

The  decay  of  the  2.578  MeV  state  is  consistent  with  the  5/2”  spin 
assignment  of  Pilt  (Pi  69).  At  4.10  MeV  bombarding  energy,  the  distri¬ 
bution  of  the  2578  keV  gamma-ray  was  isotropic  within  experimental  error 
No  conclusions  about  the  spin  of  the  level  could  be  drawn  from  this  work 
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3.6.  The  2.818  MeV  Level 


No  evidence  was  found  of  the  existence  of  this  level  reported  in 
the  work  of  Sperduto  and  Buechner  (Sp  64) .  Two  transitions  might  be 
considered  as  originating  from  this  level:  a  transition  of  1408  keV 
corresponding  to  the  decay  to  the  1.410  MeV  (7 / 2 — )  level  and  a  transi¬ 
tion  of  517  keV  corresponding  to  the  decay  to  the  2.211  MeV  (9/2“) 
level.  Because  of  multiple  feeding  of  the  low  lying  states  and  of  the 
uncertainties  in  the  branching  ratios,  these  transitions,  which  would 
be  weak  were  neither  rejected  nor  confirmed. 

3.7.  The  2.871  MeV  Level 

This  level  is  weakly  populated  in  the  55Mn(p,ny)  reaction.  It 

was  found  to  decay  88%  to  the  ground  state  and  18%  to  the  0.930  MeV 

(5/2-)  state.  This  second  branching,  not  reported  by  Pilt  (Pi  69),  was 

measured  by  Fischbeck  et_  al_.  (Fi  66)  to  have  a  strength  of  40%.  The 

angular  distributions  of  the  ground  state  transition  measured  at  4.10 

and  4.35  MeV  proton  bombarding  energy  were  isotropic  within  statistical 

error.  This  limits  the  spin  value  of  the  2.871  MeV  level  to  J  <  7/2. 

There  is  a  level  at  "2.90"  MeV  reported  by  Goodman  et_  al_.  (Go  62)  with 

a  neutron  angular  momentum  transfer  i  =  3  in  the  56Fe(p,d)  pick-up 

n 

reaction.  This  result  is  confirmed  by  Sherr  et  al .  (Sh  65).  If  this 
"2.90"  MeV  state  is  the  2.871  MeV  level,  as  suggested  by  Auble  and 
Rapaport  (Au  70),  the  present  results  would  favor  a  5/2“  spin  assign- 
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ment  as  opposed  to  the  1/2  spin  assignment  suggested  by  Goodman  et  al . 
(Go  62) . 

3.8.  The  2.938  MeV  Level 

This  level  is  also  weakly  excited  in  this  reaction.  It  was 
observed  to  decay  55%  to  the  ground  state  and  45%  to  the  1.318  MeV 
(7/2  )  state,  in  agreement  with  the  results  of  Fischbeck  et_  aM .  (Fi  66). 

The  angular  distribution  of  the  ground  state  transition  was  slightly 
isotropic  (a2  =  -0.10  +  0.05)  at  4.10  MeV  proton  bombarding  energy  and 
nearly  isotropic  (a?  =  0.06  +  0.05)  at  4.35  MeV  bombarding  energy,  limit¬ 
ing  the  spin  of  the  2.938  MeV  to  J  <_  7/2  in  agreement  with  Pilt's  results 
(Pi  69). 

Auble  and  Rapaport  (Au  70)  suggest  the  existence  of  two  levels 

around  2.94  MeV:  one  at  2.940  MeV  which  would  correspond  to  a  stripping 

level  with  a  tentative  neutron  momentum  transfer  l  =1  according  to 

n 

Sperduto  _et_  _al*  (Sp  64),  and  one  level  which  would  be  the  2.948  MeV  level 
reported  by  Fischbeck  et_  ckL.  (Fe  66),  with  spin  value  limited  to  J71  =  5/2“ 
or  7/2“.  No  evidence  of  a  second  level  has  been  found  in  this  study 
which  suggests  that  the  2.938  MeV  level  may  in  fact  be  the  2.948  MeV 
level  observed  by  Fischbeck  e_t  aM.  (Fe  66)  . 
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3.9.  The  2.984  MeV  Level 

The  primary  decay  is  by  a  1574  keV  gamma-ray  transition  to  the  1.410 
MeV  (7/2”)  level.  A  weak  2054  keV  gamma  ray  may  be  due  to  a  transition 
to  the  0.930  MeV  (5 / 2— )  level  with  a  strength  equal  to  less  than  10% 
of  the  total  strength.  A  1666  keV  gamma  ray  which  could  correspond  to 
a  transition  to  the  1.318  MeV  (7/2  )  has  been  assigned  entirely  to  the 
decay  of  3.076  MeV  level  as  explained  in  section  3.1. 

The  angular  distribution  of  the  1574  keV  gamma  ray  (fig.  16a)  is 
strongly  anisotropic  at  4.35  MeV  bombarding  energy  (a9  =  -0.70  +0.11). 

r\ 

The  analysis  (fig.  16b)  yields  a  unique  spin  assignment  J  =  9/2  and 
a  mixing  ratio  of  -0.85^*^  .  It  is  unlikely  that  this  level  has  a 
positive  parity  as  the  experimental  results  indicate  that  the  transi¬ 
tion  would  correspond  to  a  strong  M2-E1  mixing.  This  assignment  J71  = 

9/2”  is  in  disagreement  with  the  b5Co  8-decay  study  of  Fischbeck  et  al. 
(Fi  66)  who  report  no  6  decay  to  this  level  and  consequently  rule  out 
5/2”,  7/2”  or  9/2“  as  possible  spin  assignments.  However,  Haupt  e_t  al . 

(Ha  65)  in  a  similar  work  report  a  1580  keV  gamma-ray  transition  in  their 
Ge  (Li)  spectrum.  Coincidence  measurements  indicate  that  this  transition 
is  fed  by  (0.09  +  0.05)%  of  the  total  8  decay.  As  they  identified  the 
1580  keV  gamma  ray  with  a  1600  keV  seen  in  the  Nal  spectrum,  they  con¬ 
cluded  that  this  transition  was  the  1612  keV  transition  between  the 
2.542  MeV  and  the  0.930  MeV  level.  This  1580  keV  transition  is  likely 
to  be  the  1574  keV  transition  involved  here  and  the  9/2”  assignment  would 
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Figure  16  The  2*984  MeV  level:  Angular  distribution  (a)  and  X2  plots  (b) 
for  the  1574  keV  gamma-ray  transition  to  the  1.410  MeV  (7/2“) 
level . 


A1ISN31NI  3AllV13d 
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Table  V 


A  summary  of  spin  assignments  and  multipole  mixing  ratios  determined 


in  this  study  of 

the  55Mn(p,ny) 55Fe 

Reaction 

Initial  State 

Final  State 

Spin  Sequence 

(keV) 

(keV) 

Assumed 

Mixing  Ratio  Solution 

2211 

1410 

9/2  -  7/2 

n  0 , +0 . 28  .  n  c+0 . 90 

0.21  or  2.15  „ 

-0.70  -0.85 

2301 

930 

9/2  -  5/2 

E2 

2542 

1318 

11/2  -  7/2 

E2 

• 

9/2  -  7/2 

n  cri+0.20  0  --+0.75 

°'60-0.12  or-2-75-1.25 

2984 

1410 

9/2  -  7/2 

0  o  2  •  3  0 
°’85-0.70 

3076 

1410 

11/2  -  7/2 

E2 

9/2  -  7/2 

n  ,-n+0.24  o  -7C-+2.50 

°-60-2.15  °r  -2‘75-2.90 

, 


' 
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then  be  consistent  with  the  |3+-decay  of  55Co. 

3.10.  The  3.027  MeV  Level 

This  level  decays  65%  to  the  ground  state  and  35%  to  the  0.410  MeV 
(1/2”)  state.  A  weak  transition  (<  1%)  to  the  1.410  MeV  (7/2“)  was  also 
observed.  The  angular  distributions  of  the  ground  state  transition  and 
the  2617  keV  transition  to  the  0.410  MeV  were  both  isotropic  at  4.35 
MeV  bombarding  energy.  Both  the  observed  gamma  decay  and  angular  distri¬ 
butions  are  consistent  with  the  J77  =  3/2  assignment  of  this  state,  based 
on  the  radiative  neutron  capture  measurements  of  Earle  and  Bartholomew 

(Ea  66)  and  on  the  i  =1  momentum  transfer  observed  in  stripping 

n 

reaction  (Fu  63,  Ma  64). 

3.11.  The  3.076  MeV  Level 

The  1666  keV  gamma  ray  has  been  entirely  assigned  to  the  decay 
from  the  3.076  MeV  level  to  the  1.410  MeV  (7/2“)  state  (see  section 
3.1.).  The  angular  distribution  of  this  transition  is  shown  in  fig.  17a. 
It  is  strongly  anisotropic  (a2  =  0.38  +  0.06)  and  the  x2  analysis 
(fig.  17b)  yields  two  possible  assignments:  J  =  9/2  or  11/2,  the 
former  being  less  probable  with  a  minimum  x2  of  2.1.  As  no  8+  decay 
has  been  observed  to  this  level  (Ha  65,  Fi  66)>  the  11/2  assignment 


would  be  more  likely. 
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Figure  17 


The  3.076 
plots  (b) 
1.410  MeV 


MeV  level:  Angular  distribution  (a)  and 
for  the  1666  keV  gamma-ray  transition  to 
(7/2  )  state. 
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CHAPTER  IV 


THEORETICAL  INVESTIGATION  OF  53Cr  AND  55Fe 

1.  MOTIVATION 

There  is  evidence  of  considerable  fragmentation  of  the  single  particle 
states  in  the  stripping  and  pick  up  experiments  leading  to  levels  in  53Cr 
and  55Fe.  In  53Cr,  the  first  three  3/2"  states  (the  ground  state,  the 
2.324  MeV  and  the  3.625  MeV  states)  have  been  assigned  almost  the  totality 
of  the  2p3y2  strength,  and  the  1/2”  (0.564  MeV)  state  about  40%  of  the  2p \/2 
strength.  Four  5/2”  states  between  1.0  and  5.0  MeV  excitation  share  about 
50%  of  the  If 5/2  strength.  Most  of  the  lgg/2  strength  has  been  assigned  to 
the  3.715  MeV  state.  The  =  2  transfers  indicate  considerable  fragmenta¬ 
tion  of  the  2d5/2  strength  (Ra  68,  Wh  67). 

Many  inelastic  scattering  experiments  on  53Cr  have  been  reported  and 
all  show  that  the  0.564  (1/2"),  1.006  (5/2"),  1.285  (7/2")  and  1.971 
(>_3/2“)  MeV  states  are  strongly  excited  (Ra  68,  Wh  67  and  Me  66).  The 
angular  distribution  of  inelastic  alpha  scattering  of  the  first  three 
excited  states  is  indicative  of  their  description  in  terms  of  a  particle 
coupled  to  the  excited  52Cr  core  (Me  66).  The  1.285  MeV  level  has  been 
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described  by  a  (2+)  2p 3y  core-particle  weak  coupling  model  (BD  64,  BD  65) 
and  Whitten  (Wh  67)  has  assigned  a  20%  admixture  of  the  (2+)  2p3/2  con¬ 
figuration  to  the  o3Cr  ground  state  configuration.  Measurements  of 
enhanced  E2  transition  probabilities  for  the  low  lying  levels  of  53Cr 
(Me  66,  Ga  69)  indicate  the  existence  of  collective  effects. 

In  55Fe,  less  experimental  data  are  available  but  they  indicate 
that  the  general  features  observed  in  53Cr  are  present  in  this  nucleus 
also.  The  ground  state  (3/2”)  has  been  assigned  about  85%  of  the  2p.^ 
strength  and  the  0.412  MeV  (1/2-)  state  about  45%  of  the  2p1/2  strength 
(Ma  64,  Fu  63).  Four  states  between  1.0  and  4.1  MeV  share  approximately 
85%  of  the  lf5/  strength.  Most  of  the  &n  =  4  (lgg^ )  strength  seems  to 
be  assignable  to  a  level  at  3.8  MeV. 

The  results  obtained  in  this  work  have  shown  the  probable  existence 
of  two  11/2”  states  in  55Fe.  One  of  these  at  3.076  MeV  excitation  energy 
is  apparently  a  member  of  a  band  built  on  the  1.408  MeV  (7/2  )  hole 
state,  with  the  2.211  MeV  (9/2")  state  as  second  member.  On  the  other 
hand, it  seemed  possible  to  describe  many  levels  of  53Cr  and  55Fe  in  terms 
of  another  rotational  band  built  on  a  K  =  1/2”  state.  These  remarks 
suggested  a  possible  application  of  the  Nilsson  model  which  has  been 
applied  to  57Fe  and  59Fe  (So  67,  So  69)  with  some  success.  However,  56Fe 


and  58Fe  are  apparently  strongly  deformed> whereas  54Fe  as  well  as  DZCr 
exhibit  vibrational  properties  (see  e.g.  Go  62,  Fu  64).  This  means  that 
in  both  53Fe  and  55Fe,  the  Nilsson  model  has  little  chance  of  success. 


5  4n 


52/ 


Nilsson  model  calculations  will  be  described  in  section  2. 
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Calculations  involving  the  intermediate  coupling  in  the  unified  model 
have  been  done  for  N  =  29  odd  nuclei  in  1963  by  Ramavataram  (Ra  63) .  At 
the  time,  little  was  known  about  these  nuclei  outside  the  information 
garnered  from  stripping  reaction  studies.  It  was  felt  that  this  type  of 
calculation  could  be  used  again  with  different  parameters,  and  with  the 
inclusion  of  more  levels  to  be  fitted  as  well  as  with  the  calculation  of 
electromagnetic  properties  and  spectroscopic  factors.  For  this  purpose, 
two  programs  were  written  which  are  described  in  the  Appendix.  The 
model  is  also  fully  described  in  the  Appendix  and  again  outlined  in 
section  3.  Calculations  based  on  the  "Thankappan  and  True"  model  and  on 
the  shell  model  have  been  done  for  53Cr  and  55Fe.  Both  models  will  be 
briefly  described  in  sections  4  and  5.  In  section  6,  results  obtained 
from  these  models  will  be  compared  with  the  calculations  based  on  the 
intermediate  coupling  in  the  unified  model.  A  discussion  will  follow 
in  section  7. 
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2.  NILSSON  MODEL  CALCULATIONS 


2.1.  Theory 

The  Nilsson  model  has  been  extensively  described  elsewhere  (Ni  55, 

Da  65)  and  no  attempt  will  be  made  to  derive  the  results  here. 

One  can  describe  odd-A  nuclei  as  constituted  by  an  odd  nucleon 
moving  in  the  average  field  of  the  remaining  nucleons  which  constitute 
the  "core".  In  the  single  particle  shell  model,  the  average  field  is 
supposed  to  be  spherically  symmetric.  In  the  Nilsson  model,  it  is 
assumed  that  the  core  is  permanently  deformed.  Any  oscillations  in 
shape  which  take  place  have  sufficiently  low  frequency  to  permit  the 
particle  to  follow  adiabatically .  In  the  strong  coupling  limit,  the 
motion  of  the  extra  core  nucleon  is  exactly  the  same  in  the  body-fixed 
reference  frame  as  that  of  a  nucleon  moving  in  a  spatially  fixed 
potential  well.  The  average  field  is  no  longer  spherically  symmetric, 
but  has  spheriodal  symmetry  about  the  body-fixed  3-axis. 

The  total  Hamiltonian  of  the  system  contains  two  terms:  the  Ham¬ 
iltonian  of  the  odd  nucleon  and  the  kinetic  energy  of  the  core.  Consider 
first  the  single  particle  motion.  It  is  described  by  the  Nilsson 
Hamiltonian:  (Ni  55) 


(IV. 1) 


.. 
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with  %&  =  -  |  «jf  r2  Y20 


0 

/) 


se 


is  the  Hamiltonian  for  a  spherically  symmetric  harmonic  oscillator 


with  eigenvalues  Nftu)Q;  ^  represents  the  coupling  of  the  particle  to  the 
axis  of  deformation;  C  &.s  is  the  spin-orbit  term;  the  D -  term 
gives  a  correction  to  the  oscillator  potential,  especially  at  large  dis¬ 
tance.  (It  is  sometimes  referred  to  as  the  flattening  parameter  as  it 
serves  to  depress  the  high  angular  momentum  states.) 

Starting  from  this  definition  some  new  parameters  are  introduced 
which  will  express  equation  (IV. 1)  in  a  more  convenient  form  for  computa¬ 
tion.  In  Nilsson's  formalism,  equation  (IV. 1)  becomes: 


—  <^0  =  Khco0  R 


(IV. 2) 


where  R  =  qU-2£.s  -  y£2 

with  the  new  parameters  defined  by: 


< 


1 

2 


C 


y 


and  ri 


6  w0(fi) 

<  T7" 


The  dependence  of  on  the  deformation  6  arises  from  a  requirement 

0 

that  the  nucleus  volume  remains  constant.  At  6  =  0,  u)Q  (0)  =  o)0  is  related 

o  _JL 

to  the  nuclear  mass  through  the  relation  ft  =  41  A  3  MeV . 

The  values  of  the  parameters  <  and  y  are  chosen  such  that  for  zero 
deformation  the  shell  model  levels  be  reproduced.  The  total  angular 
momentum  j"  of  the  particle  is  no  longer  a  good  quantum  number  but  its 
projection  Q,  along  the  symmetry  axis,  (the  3-axis  in  the  body  fixed 
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reference  frame)  is  a  good  quantum  number. 

The  calculation  then  consists  in  the  diagonalization  of  the  matrix 
R  in  a  chosen  representation.  It  is  then  possible  to  obtain  the  eigen¬ 
functions  of  the  particle  Hamiltonian  .  The  basis  vectors  chosen  by 

P 

Nilsson  are  those  in  which  is  diagonal  together  with  £2  ,  £3  (=A)  and 
s3  (=£).  The  eigenvectors  labelled  by  ft  and  another  label  a  are  then 
expanded  in  terms  of  the  basis  vectors: 

av.3> 

£A 

For  a  given  set  of  parameters  y  and  k,  one  can  calculate  the  energy 

eigenvalues  E  of  the  Hamiltonian  'b  as  functions  of  the  deformation 
a  p 

parameters  p  or  6 .  An  example  of  such  a  calculation  for  N  =  3,  y  =  0.25 
and  k  =  0.085  is  given  in  fig.  18,  It  can  be  seen  that  the  shell  model 
single  particle  states  which  are  reproduced  for  zero  deformation  are 
split  into  states  of  different  ft  for  non- zero  deformation.  Each  of 
these  orbits  is  degenerate  in  +  ft.  The  particles  fill  the  orbits  by 
pair.  The  curves  representing  each  state  are  labelled  by  the  value  of 
ft,  the  parity  and  the  triad  [NnaA] . 

In  order  to  determine  the  total  Hamiltonian  of  the  nuclear  system, 
it  is  now  necessary  to  evaluate  the  kinetic  energy  of  the  core.  As 
mentioned  earlier,  vibrations  are  supposed  to  have  low  frequency, 
consequently,  the  vibrational  term  is  neglected  in  the  core  Hamiltonian. 
In  other  words,  it  is  assumed  that  all  the  calculations  correspond  to 
the  vibrational  ground  state.  The  core  Hamiltonian  describes  only  the 
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Figure  18  Nilsson  model  orbits  as  a  function  of  the  deformation  paramete: 
6,  for  y  =  0.25  and  k  =  0,085. 

The  value  of  is  marked  on  the  right  side  of  the  figure  as 
well  as  the  Nilsson  level  asymptotic  quantum  numbers  [Nn3A] . 

The  orbit  numbers  are  those  given  by  Nilsson  (Ni  55)  .  The 

0  -1 

energy  scale  was  obtained  with  hw0  =  41(55) 
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rotation  of  the  core  and  the  total  Hamiltonian  can  then  be  written  as 


+  j2  _ 

T  p  23  L  J 


2(1.1)] 


(IV. 4) 


where  is  the  single-particle  Hamiltonian  and  the  following  term 

describes  the  core  rotation  and  the  coupling  between  rotator  and 
particle . 

The  total  angular  momentum  I  =  R  +  j  with  R  and  j  being  the  core 
and  particle  angular  momentum,  respectively,  (see  fig.  19).  0  is  the 
moment  of  inertia  of  the  rotating  core. 

The  (I.j)  term,  called  either  Coriolis  or  rotation-particle  coupling 
[RPC]  term,  results  from  the  coupling  between  the  intrinsic  motion  of  the 
nucleon  and  the  core  rotation.  In  most  cases,  the  RPC  term  has  a  very 
small  effect  for  low-lying  states  of  an  axially  symmetric  system  and  can 
be  neglected  in  a  first  order  approximation  except  when  K,  the  projection 
of  I  on  the  3-axis,  is  equal  to  1/2.  In  this  case,  the  RPC  term  does 
contribute  directly  to  the  energy  of  states. 

The  total  wave  function  describing  the  state  of  the  nucleus  is  a 
product  of  the  rotation  wave  function  and  the  single  particle  wave 
function.  It  can  be  written  as: 


EIMK>  =  l  C  [ | IMK>  |jfi>  +  (-)1  j  | IM-K>  |j-fl>].  (IV. 5) 

v2  j 


The  eigenvalue  problem  of  the  rotation  Hamiltonian  is  treated  in 
most  textbooks  [see  e.g.  (Pr  62)].  It  can  be  shown  that  the  total  energy 
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Figure  19  Angular  momentum  diagram 

— y 

I  the  total  angular  momentum  is  equal  to  R  +  j  where  R  is  the  angular 
momentum  generated  by  the  collective  motion  of  the  nucleus  and  j  is 
the  angular  momentum  of  the  odd-nucleon,  z  is  the  quantization  axis, 
and  3  is  the  body-fixed  symmetry  axis. 

In  the  ground  state,  R  is  perpendicular  to  3,  (ft  =  K),  i.e.  The 
collective  rotation  takes  place  around  an  axis  perpendicular  to  the 


nuclear  symmetry  axis. 
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eigenvalue  is  given  by: 

2  1 
Ejk  =  Ek  +  Y,  [1(1+ 1)  -  2K2  +  6k  L  a(~)I+T  (I+j)]  (IV. 6) 


06 

where  is  the  energy  of  the  odd  nucleon  in  a  Nilsson  orbit,  and 

a  is  the  "decoupling  parameter"  resulting  from  the  RPC  term,  defined  by: 


Physically,  this  means  that  rotational  bands  are  built  on  each 
single-particle  state.  The  core-momentum,  R,  has  no  component  on  the 
3-axis.  K  is  entirely  due  to  the  particle  angular  momentum  and  is 
constant  throughout  the  band.  All  higher  values  J  >_  K  are  permitted, 
and,  when  K  ^  y,  the  lowest  state  has  I  =  K  and  successive  states  have 
I  =  K  +1,  K  +  2,  .  .  . 

The  vibration-rotation  interaction  can  be  taken  into  account  by 
adding  to  equation  IV. 6  a  term  which  has  the  form:  (Da  65) 

B  [  I  (I  +  1)  +  6  v  1.  a(-)I4T  (I  +  h)2 
R,  2  £ 


The  coefficient  B  is  usually  small  and  can  be  treated  as  an  adjustable 
parameter . 

If  the  particle  is  rather  loosely  bound,  then  the  K-bands  overlap 
strongly  and  the  corrections  due  to  the  RPC  term  become  important.  The 
total  Hamiltonian  then  includes  non-negligible  off-diagonal  matrix 
elements  and  after  diagonalization ,  the  eigenstates  are  a  mixture  of 
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different  K-bands.  Obviously,  K  is  no  longer  a  good  quantum  number  in 
this  case.  Kerman  (Ke  56)  has  considered  the  problem  of  two-band  mix¬ 
ing.  Brockmeier  et  al.  (Br  65)  have  generalized  it  to  many  bands.  The 
band  mixed  states  can  then  be  written  as: 


Iv>  =  l  Cv  (I , K)  |K> 
K 


(IV. 7) 


where  v  is  an  additional  label  for  states  with  same  I.  . 

2.1.1.  Electromagnetic  transitions 

Once  the  total  eigenvalue  problem  has  been  solved,  one  can  calcu¬ 
late  the  electromagnetic  transition  rates  and  static  moments. 

The  procedure  to  follow  is  similar  to  the  procedure  described  in 
the  Appendix  in  the  case  of  the  intermediate  coupling  model.  Brockmeier 
et  al.  (Br  65)  have  calculated  the  transition  rates  formulae  in  the 
band  mixing  case.  For  the  electric  quadrupole  reduced  transition,  they 
obtain: 


vv 


(IV. 8) 


The  matrix  elements  b„  can  be  calculated  from  the  single  particle 

E2 


KK' 


wave  functions  (Br  65)  .  Q 


,  the  matrix  element  of  the  quadrupole 


' 
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operator  is  composed  of  two  terms,  the  intrinsic  quadrupole  moment  Qq 

and  a  single  particle  term,  .  For  nuclei  with  a  large  permanent 

KK 1 

deformation,  Qq  >>Q^,  hence  Q  Q0 .  The  intrinsic  quadrupole 

moment  is  given  as  a  function  of  the  nuclear  deformality,  6,  by:  (Ni  55) 

Q0  —  0 . 8ZRq2  6(1+|  6) 


where  Z  is  the  atomic  number  and  R0 ,  the  nuclear  radius,  is  taken  to  be 
1.2  fm. 

The  relation  between  the  measured  quadrupole  moment,  Q  ,  and  Q0 

s  * 

is  given  by:  (Bo  57) 


s  3K2  -  1(1+1) 
s  (1+1)  (21+3) 


(IV. 9) 


2.2.  Application  of  Nilsson  Model  to  53Cr  and  55Fe 

These  calculations  can  be  divided  into  two  parts.  First,  some 
assumptions  are  made  about  the  possible  rotational  bands  in  the  nucleus 
under  study.  This  allows  the  determination  of  the  rotational  parameter, 

A  =  h2/23,  and,  in  the  case  of  K  =  1/2  bands,  of  the  decoupling 
parameter  a.  Then,  one  tests  these  assumptions  by  including  the  Nilsson 
single  particle  energy  levels  and  by  taking  into  account,  if  necessary, 
the  K-band  mixing.  From  the  results  thus  obtained,  the  Nilsson- 
parameters  are  either  readjusted  or  other  bands  are  considered. 

In  the  program  used  for  these  calculations,  the  different  parameters 
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< ,  y and  6  could  be  varied  freely.  For  a  given  range  of  these  parameters 
a  x2  routine  determined  the  value  of  each  parameter  yielding  the  best 
fit  to  a  given  number  of  experimental  levels.  Band  mixing  was  also 
considered.  As  suggested  by  Brockmeieret  al.  (Br  65),  the  off-diagonal 
RPC  terms  were  allowed  to  vary  freely  in  the  band  mixing  calculations. 

The  fact  that  in  55Fe,  two  sets  of  levels  can  apparently  be 
described  by  rotational  bands  suggested  the  possibility  of  trying 
Nilsson  model  calculations.  The  first  band  is  built  on  a  K  =  1/2 
ground  state.  With  a  decoupling  parameter  a  =  -1.86  and  a  rotational 
parameter  A  =  159  keV  the  ground  state  (3/2”),  the  0.410  MeV  (1/2"), 
the  1.318  MeV  (7/2  )  and  the  2.144  MeV  (5/2  )  level  are  well  reproduced. 
On  the  other  hand,  the  2.301  MeV  (9/2”)  and  the  3.076  MeV  (assumed  to 
be  11/2  )  could  be  members  of  a  second  band  built  on  K  =  7/2".  These 
considerations  imply  that  the  odd  neutrons  should  be  in  the  Nilsson 
orbit  l/2— [321]  in  its  lowest  configuration  and  that  the  deformation 
should  be  positive.  There  are  other  arguments  favoring  positive 
deformation:  negative  deformation  would  result  in  many  hole-state  bands 
built  on  the  1/ 2 [ 330 ] ,  3 / 2 [ 321 ] .  .  .  Nilsson  orbits  in  contradiction 
with  the  fact  that  the  first  strong  hole  state  has  a  spin  J71  =  7/2  , 
impossibility  of  describing  the  levels  of  55Fe  with  a  ground  state  band 
built  on  the  3/2[312]  Nilsson  level,  and  disagreement  with  the  probable 
negative  sign  of  the  quadrupole  moment. 

Starting  from  these  considerations,  it  was  then  necessary  to  deter¬ 
mine  the  parameters  y  and  k.  In  his  paper,  Nilsson  (Ni  55)  fixes  the 


i 
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value  of  y  at  0.35.  This  corresponds  to  a  (f5/_,p3/  )  energy  separation 

'2  2 

which  is  about  half  of  the  (pw  -pq/  )  energy  separation.  [k  determines 

i/ 2  6'2 

the  absolute  magnitude  of  these  energies.]  The  order  of  these  shell 
model  levels  agrees  with  experimental  results  when  the  odd  nucleon  is 
a  proton.  When  the  odd  nucleon  is  a  neutron,  little  is  known  when  the 
f-shell  has  been  filled  and,  in  theoretical  calculations,  the  energy 
separation  of  the  single  particle  states  is  usually  taken  as  a  free 
parameter . 

A  Nilsson  level  diagram  is  shown  in  figure  17.  It  corresponds  to 
55Fe  for  y  =  0.25  and  <  =  0.085.  Those  values  are  close  to  those 
obtained  for  the  best  fits  of  53Cr  and  55Fe. 

The  final  sets  of  parameters  used  in  the  calculation  of  53Cr  and 
55Fe  are  given  in  Table  VI.  The  corresponding  level  diagrams  are 
shown  in  figure  20.  Band  mixing  was  allowed  in  the  calculations. 


Table  VI 

Parameters  corresponding  to  the  best  fits  of  the  energy  levels 

in  53Cr  and  55Fe 


53Cr 


5  5Fe 


y  = 

0.251 

k  =  0.085 

n 

= 

1.35 

K  = 

1/2  band 

[g.s.  band] 

A 

= 

166  keV 

K  = 

3/2  and  K 

=5/2  bands 

A 

= 

166  keV 

K  = 

7/2  band 

A 

= 

77  keV 

y  = 

0.244 

k  =  0.080 

n 

= 

1.44 

K  = 

1/2  band 

[g.s.  b  and ] 

A 

= 

159  keV 

K  = 

3/2  and  K 

=  5/2  bands 

A 

= 

159  keV 

K  = 

7/2  band 

A 

= 

85  keV 

6  =  0.115 
a  =  -2.0 


6  =  0.115 
a  =  -1.86 
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Figure  20  Nilsson  model  calculations:  Theoretical  and  experimental  spectra 
are  shown  for  53Cr  and  55Fe.  The  parameters  used  for  these 
calculations  are  given  in  Table  VI.  All  states  have  negative 
parity . 
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2.2.1.  Discussion 

The  agreement  between  experimental  and  theoretical  level  schemes  is 
not  good  although  some  slight  improvement  was  obtained  when  the  magnitude 
of  the  RPC  term  was  slightly  increased.  The  main  problem  lies  with  the 
5/2  levels  which  cluster  around  2.5 MeV  for  both  33Cr  ancj  55pe  nuciei.  All 
attempts  to  try  to  reproduce  the  first  5/2  state  at  —  1  MeV  resulted  in 
an  overall  poorer  fit. 

The  reduced  transition  probabilities  for  electric  quadrupole  transi¬ 
tions  to  ground  of  the  first  1/2“  and  7/2“  states  of  53Cr  were  calculated 
to  be  152  and  69  e2-fm4  respectively.  These  values  are  to  be  compared 
with  the  experimental  values  of  (200  +  20)  and  (95  +  25)  e2-fm4, 
respectively,  obtained  from  an  average  of  the  results  reported  in  (Me  66) 
and  (Ga  69).  The  electric  quadrupole  moment  of  53Cr  was  calculated  to  be 
-0.12  b to  be  compared  with  the  experimental  value  of  |q|  =  0.03  b  (Fu  69). 

As  explained  earlier,  (section  1),  the  apparent  failure  of  describing 
53Cr  and  55Fe  by  the  Nilsson  model  was  somewhat  expected.  In  52Cr  and 
54Fe,  collective  features  are  observed  but  they  are  better  described  by 
a  vibrational  model  rather  than  by  a  rotational  model  (Go  63) .  Conse¬ 
quently,  the  Nilsson  model  based  on  the  assumption  that  the  core  is 
strongly  deformed  fails  to  reproduce  the  energy  spectra  of  53Cr  and  55Fe. 
Band-mixing  corrections  which  take  into  account  the  fact  that  the  odd 
nucleon  is  not  strongly  bound  to  the  core  are  expected  to  reproduce  cases 
where  the  decoupling  is  partial.  If  the  particle  is  too  loosely  bound 
to  the  core,  the  problem  must  be  completely  reformulated.  The  intermediate 


' 

, 


coupling  approach  seemed  to  be  the  best  alternative.  It  is  described  in 
the  following  and  will  yield  results  much  more  in  agreement  with  the 


physical  situation. 
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3.  THE  INTERMEDIATE  COUPLING  IN  THE  UNIFIED  MODEL 


3.1.  Introduction 


In  recent  years,  considerable  attention  has  been  given  to  the  descrip¬ 
tion  of  odd-A  nuclei  in  which  the  levels  are  expected  to  be  described  by 
a  model  which  considers  a  single  nucleon  moving  in  a  potential  field 
generated  by  the  even  core,  (the  remaining  nucleons),  which  may  be  in  its 
ground  state  or  in  an  excited  state. 

There  are  two  ways  to  approach  the  calculations  in  the  intermediate 
coupling  model.  The  first  approach  was  outlined  by  Bohr  and  Mottelson 
(Bo  52,  Bo  53)  and  later  elaborated  by  Choudhury  (Ch  54)  and  Ford  and 
Levinson  (Fo  55) .  It  assumes  that  the  interaction  between  core  and 
particle  is  known  and  results  from  the  coupling  of  a  single  particle,  which 
has  available  different  orbital  states,  to  a  core  capable  of  performing 
quadrupole  vibrations.  This  model  is  known  as  the  intermediate  coupling  in 
the  unified  model  and  will  simply  be  labelled  as  "intermediate  coupling 
model"  in  the  following.  The  second  approach  was  suggested  by  Lawson  and 
Uretsky  (La  57)  and  de-Shalit  (Sh  61),  and  has  been  developed  by  Thankappan 
and  True  (Th  65).  It  makes  no  assumption  on  the  nature  of  the  excited 
states  of  the  even  core  and  introduces  a  general  type  of  interaction.  This 
approach,  now  referred  to  as  the  "Thankappan  and  True"  (T.T.,  hereafter) 
model  will  be  briefly  outlined  in  section  4. 


. 


::  3 
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3.2.  Outline  of  the  Intermediate  Coupling  Model 

In  this  model,  the  odd-A  nuclear  system  is  described  by  the  coupling 
of  the  odd  nucleon  to  the  remaining  nucleons  constituting  a  core  capable 
of  performing  quadrupole  oscillations.  The  complete  formalism  of  this 
model  is  given  in  the  Appendix.  Only  the  main  features  will  be  recalled 
here . 


The  total  Hamiltonian  of  the  system  consists  of  three  terms:  the 

even  core  Hamiltonian,  the  single  nucleon  Hamiltonian,  >  ,  and  the 

interaction  Hamiltonian,  U  .  .  That  is: 

int 


<D 

9  =  +  'O  + 

c  p 


j  r\i  * 

vJr. 


int 


(IV. 10) 


The  interaction  Hamiltonian,  ,  is  given  by: 


<\Q 

A3 


int 


=  -5^,/t  l  {b9„  +  <“>V  b2,-u}Y2y  (6,(j)) 


\/  5  L  2y 

y 


in  which  £  is  the  dimensionless  coupling  strength,  b  and  b*  are  boson 
creation  and  annihilation  operators,  ^  ^  i-s  t^ie  spherical  harmonic 
operator  and  ftu)  is  the  phonon  energy. 

The  representation  chosen  is  the  representation  for  which  +  ' >  ) 

is  diagonal.  The  resulting  basis  vectors  are  defined  by  | j  NR  IM>  and, 
(neglecting  the  zero  order  energy  point),  satisfy: 


+t  )  |j  NR  IM>  =  (Ntiw  +  Ei)  |  j  NR  IM> 

CD  J 


' 


. 


' 
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where  j  and  R  are  the  angular  momenta  of  the  odd  nucleon  and  the  vibrating 
core,  respectively;  I  =  R  +  j  is  the  total  angular  momentum  of  the  nucleus, 
M,  its  z-component;  N  is  the  number  of  phonons  of  surface  oscillations; 

Ej  is  the  energy  of  the  odd  nucleon  in  a  state  j . 

The  problem  is  then  to  express  the  total  Hamiltonian,  which  after 
diagonalization  will  yield  eigenvalues  and  wave  functions.  Once  the  wave 
functions  have  been  obtained,  it  is  possible  to  calculate  transition  rates, 
static  moments,  etc.  (see  Appendix). 

3.3.  Application  of  the  Intermediate  Coup 1 ing  Model  to  the  Study  of  53Cr 

and  55Fe 

In  the  study  of  5^Cr  and  53Fe,  the  29^h  neutron  was  considered  to  have 
available  the  2p3/2  ,  2p1^2  and  lf5y,  states.  The  single  particle  energies, 
labelled  as  Ej  =  P;/2“P3/2  an<^  e2  =  ^S/z "P3/2*  were  considered  as  free 

parameters;  the  single  particle  spacings  in  ^Ca,  (p ,  -p,  ,  =  2.02  MeV  and 

1  /  2  '2 

f 5  ■  -p3  =  3.95  MeV)  were  taken  as  initial  values.  The  remaining  free 

parameters  are:  Tito,  which  was  taken  as  the  energy  of  the  first  2+  state  in 
the  even  core  and  later  varied  within  reasonable  range,  and  finally  the 
coupling  strength 

It  was  found  that  an  accurate  determination  of  the  different  parameters 
is  not  necessary  to  give  a  good  description  of  the  first  excited  states  of 
53Cr  and  55Fe.  In  Table  VII  are  listed  different  sets  of  parameters  which 
resulted  in  almost  Identical  energy  spectra.  Figure  21  shows  how  the  energy 


. 


(?£ 


■ 


Table  VII 


Sets  of  parameters  used  in  the 
intermediate  coupling  model  calculations 


Nucleus 

Set  £ 

l-Pl/2  -P3/2 
MeV 

£2"f5/2"P3/, 

MeV 

hm 

l 

MeV 

E  (2"*p  exp* 
MeV 

K 

5  3Cr 

A 

2.0 

2.2 

1.5 

1.43 

2.5 

B 

1.8 

2.0 

1.6 

- 

2.25 

55Fe 

C 

2.0 

2.2 

1.45 

1.40 

2.75 

D 

- 

- 

- 

- 

3.00 

E 

2.2 

2.0 

1.40 

- 

2.75 

JL  i 

This  value  is  the  excitation  energy  of  the  first  2~r  level  in 
the  neighbouring  even  nucleus. 
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Figure  21  The  intermediate  coupling  model:  Energy  levels  of  the  total 
Hamiltonian  as  a  function  of  the  coupling  parameter  E, » 


> 
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levels  of  the  total  Hamiltonian  vary  when  the  coupling  parameter  £  increases „ 
(The  remaining  parameters  are  those  of  set  C  in  Table  VII.) 

The  best  fits  were  determined  after  comparison  with  all  the  experimental 
information  available.  In  Table  VII,  they  correspond  to  set  B  for 53Cr  and  to 
set  E  for  55Fe.  Tables  VIII  and  IX  give  the  resulting  wave  functions  for 
low-lying  states  of  53Cr  and  55Fe.  The  complete  theoretical  results  will  be 
presented  and  compared  with  experimental  data  as  well  as  with  other  theoreti¬ 
cal  calculations  in  section  6. 

3.4.  Justification  of  the  Use  of  this  Model 

The  main  assumption  of  this  model  is  based  on  the  existence  of  a 
vibrational  core.  There  seems  to  be  great  controversy  in  the  literature 
on  the  structure  of  the  even  nuclei  in  the  f-shell.  From  a  strict  shell 
model  view  point,  most  of  the  levels  of  a  nucleus  such  as  52Cr  result  from 
the  different  configurations  obtained  from  the  coupling  of  4  f7^  protons. 

With  the  addition  of  the  seniority  quantum  number,  most  of  the  levels  are  in 
fact  reproduced,  but  levels  such  as  the  2.65  MeV  (0+)  in  52Cr  cannot  be 
explained  from  this  picture.  A  second  problem  arises  in  the  calculation 
of  the  electromagnetic  transitions  as  the  shell  model  usually  fails  to 
reproduce  the  enhanced  E2  transitions «  This  problem  however,  is  easily 
solved  by  introducing  the  concept  of  "effective  charge"  (see  e.g.  (Fe  69, 

Ri  69)). 

The  simple  vibrational  model  introduced  by  Bohr  (Bo  52)  predicts  levels 
which  are  equidistant  and  degenerate  for  each  phonon  group.  Recent  calcula- 


Table  VIII 
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Expansion  coefficients  corresponding  to  states  |  E  (keV)  ;  Jtt  >  of  53Cr 


calculated  from  the  intermediate  coupling  model  (Set  B  of  Table  VIII) 


Basic  States 
|NR; j> 

1 0 ;  3  /  2 — > 

| 536 ; l/2“> 

| 1174; 5/2“> 

| 1358 ; 7 / 2"> 

| 1850 ; 5/2“> 

| 2011 ; 3/2“> 

00; 3/2 

0.8582 

0.4077 

12; 3/2 

-0.3884 

0.6325 

0.5743 

0.8393 

0.6273 

0.7549 

20 ; 3/2 

0 . 1115 

-0,3036 

22 ; 3/2 

0.0521 

-0.2453 

0.0593 

-0.1403 

0.4715 

-0.076 

24 ; 3/2 

-0.3926 

-0.3730 

0.1046 

30 ; 3/2 

-0.0172 

0.0277 

32 ; 3/2 

-0.0324 

0  o  0804 

0.0627 

0.0993 

0.0499 

0.1405 

33 ; 3/2 

-0.0005 

-0.0157 

-0.0246 

0.0867 

0.0839 

34 ; 3/2 

0.0946 

-0.1104 

36 ;  3/2 

00; 1/2 

0.6073 

12; 1/2 

-0.2374 

-0.3721 

0.4030 

-0.1009 

20 ; 1/2 

22 ; 1/2 

0.0940 

0.1855 

0.1575 

0,0366 

-0.3037 

24; 1/2 

-0.2875 

30 ; 1/2 

-0.0521 

32; 1/2 

-0.0307 

-0.0740 

0.0742 

0.0516 

33 ; 1/2 

0.0434 

-0.0273 

0.1035 

34; 1/2 

0.0420 

36 ; 1/2 

0.0833 

00 ; 5/2 

-0.4932 

-0.0223 

0.2571 

12 ; 5/2 

0.1308 

0.3284 

0.2205 

-0.1807 

-0.0407 

20; 5/2 

-0.1290 

0.0512 

22 ;  5/2 

-0.0073 

0.1248 

-0.1135 

0.0780 

0.0743 

-0.1046 

24 ; 5/2 

-0.1074 

-0.0442 

0.1081 

0.2430 

0.0716 

30 ; 5/2 

0.0245 

-0.0327 

32 ; 5/2 

0.0142 

0.0609 

0.0368 

-0.0057 

-0.0271 

-0.0548 

33 ; 5/2 

0.0047 

0.0143 

0.0013 

-0.0036 

0.0544 

-0.0761 

34 ; 5/2 

0.0324 

0.0590 

-0.0277 

0.0164 

-0.0760 

36 ; 5/2 

-0.1225 

. 

Table  IX 
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Expansion  coefficients  corresponding  to  states  |e (keV) ;  Jtt>  of  j5Fe 
calculated  from  the  intermediate  coupling  model  (Set  E  of  Table  VII) 

Basic  States 

| NR ; j  >  I  0; 3/2->  I  474 ; l/2“>  1 1013 ; 5/2~>  |ll44;7/2"> 


00 ; 3/2 
12 ; 3/2 
20 ; 3/2 
22 ; 3/2 
24 ; 3/2 
30 ; 3/2 
32  ;  3 / 2 
33 ; 3 / 2 
34; 3/2 
36 ; 3/2 
00 ; 1/2 
12 ; 1/2 
20;l/2 
22 ; 1/2 
24; 1/2 
30; 1/2 
32 ; 1/2 
33 ; 1/2 
34; 1/2 
36 ; 1/2 
00 ; 5/2 
12 ; 5/2 
20 ;  5/2 
22 ; 5/2 
24 ; 5/2 
30 ; 5/2 
32 ;  5/2 
33 ; 5/2 
34 ; 5/2 
36 ; 5/2 


0.8193 

-0.4336 

0.1450 

0.0630 

-0.0241 

-0.0480 

-0.0005 


-0.2417 

0.1135 

-0.0414 


0  o 1430 
-0.0063 

-0.1302 

0.0205 

0.0066 

0.0455 


0.6639 

-0.3029 

0 . 1096 

0.5156 

0.2109 

-0.0679 

0.3343 

-0.1485 

0.0799 

0.0205 


0.5883 

0.0775 

-0.4422 

0.0819 

-0.0211 

0.1172 


-0.3444 

0.1678 


-0.0862 

0.0539 


-0.4261 

0.2147 

-0.1447 

-0.1252 

-0.0422 

0.0290 

0.0450 

0.0012 

0.0734 


0.8101 

-0.1583 

-0.4072 

0.1193 

-0.0257 

0.0483 


-0.2880 


-0.0307 

0.0961 


-0.0222 

0.0844 

0.1178 

-0.0062 
-0.0064 
-0.0309 
-0  o 1411 
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tions  where  anharmonic  terms  are  introduced  remove  the  degeneracy  (A1  69, 

Sh  62) .  Some  calculations  for  j2Cr  have  been  reported  (Wi  62)  and  the 
agreement  is  good  to  the  N  =  2  phonon  states.  However,  the  2.37  MeV  (4^*) 
which  has  seniority  v  =  4  is  not  predicted  in  these  calculations .  In 
5l+Fe,  the  recent  0+  assignment  to  the  2.56  MeV  level  gives  a  good  support 
to  the  vibrational  description  (Kr  69) . 

Another  advantage  of  the  description  in  terms  of  vibrational  model 
is  the  normal  explanation  of  E2  enhancement  through  collective  effects. 

Without  trying  to  enter  the  controversy,  it  can  be  said  that  the 
intermediate  coupling  model  based  on  the  coupling  of  a  nucleon  to  a 
vibrating  core  should  be  quite  applicable  to  the  description  of  53Cr  and 
55Fe.  For  more  than  N  =  2  phonon  states,  the  resulting  multiplets  are 
probably  quite  spread  and  there  is  not  much  experimental  support  to 
the  existence  of  most  of  N  =  3  phonon  states.  However,  as  can  be  seen 
from  Tables  VIII  and  IX,  N  =  3  phonon  state  components  in  the  resulting 
wave  functions  of  calculated  levels  of  53Cr  and  55Fe  are  neglible  for 
levels  with  excitation  below  cr  3  MeV. 

3.4.1.  Additional  remarks 

In  the  study  of  inelastic  scattering  in  the  framework  of  the  collective 
model,  the  reduced  transition  probability  for  electric  quadrupole  radiation, 
B(E2),  can  be  expressed  as:  (Sa  66,  Ba  62) 

2 

B(E2)  =  ~  ZR2  ^  e2 

4tt  5 


(IV. 11) 


' 
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_L 

where  Z  is  the  atomic  number,  R  is  the  nuclear  radius  equal  to  1.2  A  3  fm 
and  32  is  the  deformation  parameter  introduced  by  Bohr  (Bo  52) . 

In  terms  of  the  collective  constants  previously  defined,  32  is  given 
by:  (equation  71  in  (Bo  52)) 


$2 


(IV. 12) 


It  is  interesting  to  compare  the  calculated  values  of  32  using  equation 
(IV. 12)  and  the  parameters  determined  in  this  work,  (Table  VII)  with  the 
experimental  values  obtained  from  inelastic  scattering  experiments  (Ri  69) . 
For  52Cr  the  calculated  value  32  is  equal  to  —  0.16,  the  experimental  32 
ranges  between  0.156  and  0.20.  For  51+Fe  the  calculated  32  is  t^-0.17, 
compared  with  the  experimental  value  of  ~ 0 . 14 .  The  agreement  is  quite 
good  and  suggests  that  the  parameters  determined  in  these  calculations 
fall  in  the  correct  range. 
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4.  THE  THANKAPPAN  AND  TRUE  MODEL 


As  indicated  earlier,  this  model  also  considers  the  coupling  of  the 
odd  nucleon  to  the  neighbouring  even  core,  but  the  approach  is  more  general 
as  no  specific  assumption  is  made  on  the  nature  of  the  interaction, 
de  Shalit  (Sh  61)  assumed  that  under  its  most  general  form,  the  coupling 
between  the  particle  and  the  core  can  be  written  as  a  sum  of  scalar  product 
of  tensors  of  rank  k.  That  is: 


%  =  y  T  (k)  .  (k) 

int  L  c  p 

k  v 


(IV. 13) 


where  T  operates  only  on  the  degrees  of  freedom  of  the  core  and  T 

c  p 

on  the  degrees  of  freedom  of  the  particle. 

On  this  basis,  Thankappan  and  True  write  the  interaction  Hamiltonian 


as : 


=  -  ccj  *1 )  ■ 

c  p 


n 


<vy 


(IV. 14) 


where  J  and  j  are,  respectively,  the  total  angular  momentum  of  operators 
c  p 

for  the  core  and  particle  and  and  are,  respectively,  the  mass  quad- 
rupole  moment  operators  of  the  core  and  particle.  E,  and  r\  are  parameters 
describing  the  strength  of  the  interaction. 

The  resulting  matrix  elements  of  the  total  Hamiltonian  are  functions 
of  £,  q  and  two  reduced  matrix  elements  of  the  operators  Q,  namely: 


■ 
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.J'cllQjIV  and  <J  Tp||QpJ|jp>*  The  particle  reduced  matrix  element  for 

is  evaluated  using  harmonic  oscillator  wave  functions  for  the  radial  integral 

part.  The  reduced  matrix  elements  <J’  II Q  ||j  >  are  considered  as  free  para- 

c  c  c 

meters,  as  the  exact  nature  of  the  core  states  is  not  supposed  to  be  known. 

It  must  be  emphasized,  however,  that  because  of  the  form  of  the  operators 
used  in  the  calculations,  this  approach  is  unlikely  to  be  valid  unless  the 
core  states  are  collective  in  nature  (Tr  65) . 

In  order  to  avoid  dealing  with  an  excessive  number  of  free  parameters, 
it  is  necessary  to  limit  the  number  of  core  levels  taken  into  consideration. 
In  the  calculations  which  will  be  reported  in  this  work  (section  6),  only 
the  core  ground  state  and  first  excited  state  are  taken  into  account.  As 
for  the  intermediate  particle  model  calculations,  the  particle  has  avail¬ 
able  three  single  particle  states:  lpq  /  ,  lp^  ,  and  2ft-/  .  This  results 

0 / 2  >2  <2 

in  calculations  with  five  free  parameters,  at  least,  as  more  free  para¬ 
meters  (e.g.  effective  charge)  are  introduced  in  the  calculations  of  tran¬ 
sition  rates. 


. 


B I ' 
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5.  SHELL  MODEL  CALCULATIONS  FOR  N  =  29  ODD  NUCLEI 


In  the  shell  model  calculations  reported  in  this  work,  the  t+8Ca  is 


considered  as  an  inert  core.  The  Hamiltonian  for  the  ( tt  If  .  )m(v  2p  .  ) 

/ / 2  8/2 

configuration  outside  the  ^Ca  core  can  be  written  as: 


fc  =  fc  +  l$  +%  +  i  €  +  l  % 

C  L  TT-C  V-C  “  do  tt  —  7T  L.  TT-V 

it  pairs  pairs 


pairs 


(IV. 15) 


where  the  indices  c,  tt  and  v  stand,  respectively,  for  core,  proton  and 

neutrons;  m  is  the  number  of  protons  in  the  lf7/  shell  for  the  nucleus 

/  2 

under  consideration;  the  neutron  is  represented  in  its  ground  state 
configuration  2p7  but  higher  configurations  are  also  taken  into  account 
c  is  the  total  interaction  between  the  nucleons  inside  the  core.  The 
second  and  third  term  of  equation  (IV. 15)  represent  the  interactions  of 
the  extra-core  protons  and  neutrons;  the  fourth  and  fifth  terms,  the 
residual  interactions  between  the  extra  core  protons  and  neutron-proton 
pairs . 

The  first  three  terms  contribute  to  a  constant  energy  to  all  the 
states  of  the  (tt  lf7/  )m(v  2p  ,  )  configuration.  The  energy  matrix  can 
then  be  written  as: 


E ( J)  =  const.  +  E  (J  v  ) 6  ,6  ,  +  m  E  (J  v  ,J'v’,j  ;J) 

TT  TT  TT  J  J  VV  TTV  TT  TT  TT  TT  JV 

TT  TT  TT  TT 


iL'-.'JU 
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This  expression  gives  the  energy  of  a  level  of  total  angular  momentum 
J.  is  the  total  angular  momentum  of  the  m  protons;  v^  is  the  seniority 

of  the  proton  group.  is  the  total  angular  momentum  of  the  29t'1  neutron. 

is  the  contribution  of  the  proton  group  and  E  is  the  contribution  of 
a  single  neutron  proton  pair  to  the  matrix  elements. 

In  order  to  evaluate  E^,  it  is  necessary  to  know  the  effective  proton- 
proton  interaction.  The  two-body  matrix  elements  of  this  interaction  can 
either  be  evaluated  from  the  experimental  spectra  of  the  neighbouring  N  = 

28  nuclei  or  calculated  theoretically  [see  e.g.  (Ku  68)]. 

The  neutron-proton  energy  term  is  calculated  by  assuming  some  type 
of  neutron-proton  interaction.  Ohnuma  (Oh  66)  introduces  a  neutron- 
proton  interaction  of  the  form: 


V 

VTT 


v0  f(r  )  { (1-a)  +  a  a.a  } 

u  VTT  V  IT 


where  a  is  an  adjustable  parameter  and  the  radial  function  f  (r  )  is 
gaussian.  Vervier  (Ve  66)  uses  the  same  type  of  interaction  but  his 
radial  function  is  a  6-function. 

These  calculations  do  not  take  into  account  proton  configurations 
other  than  (it  lfy^)™,  that  is,  configurations  arising  from  core 
excitations.  As  a  result  states  in  N  =  29  nuclei,  which  are  predomin¬ 
antly  f  /  hole  states  will  not  be  reproduced. 

7/2 
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6.  COMPARISON  BETWEEN  THEORETICAL  CALCULATIONS  AND  EXPERIMENTAL  DATA 


In  the  following,  the  results  of  the  calculations  based  on  the  inter 
mediate  coupling  model  described  in  section  3  and  in  the  Appendix  are 
compared  with  the  experimental  results  obtained  in  this  work  or  elsewhere 
Other  calculations  will  also  be  presented  for  comparison  purposes.  They 
are,  respectively: 

i)  the  calculations  of  Ramavataram  (Ra  63)  also  based  on  the  inter 
mediate  coupling  model,  but  dealing  with  different  parameters  and 
limited  to  levels  with  spin  J  <_  7/2; 

ii)  the  shell  model  calculations  of  Vervier  (Ve  66)  and  Ohnuma 
(Oh  66,  Oh  70)  described  in  section  5.  Vervier  limited  his  calcu¬ 
lations  to  levels  with  spin  J  <_  7/2;  in  Ohnuma 's  calculations,  J 
is  limited  to  a  maximum  value  of  13/2; 

iii)  the  calculations  of  Philpott  and  True  (Ph  70),  for  53Cr,  and 
of  Larner  (La  70),  for  55Fe,  based  on  the  True  and  Thankappan  (TT) 
model  described  in  section  4. 

In  all  these  calculations,  levels  which  have  a  predominantly  hole  state 
configuration  are  not  predicted.  In  the  following,  calculated  energy 
levels,  branching  ratios,  mixing  ratios,  electromagnetic  properties  and 
spectroscopic  factors  are  compared  with  experimental  data. 
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6.1.  Energy  Spectra 

6.1.1.  53Cr  spectrum 

The  experimental  spectrum  of  53Cr  is  compared  with  intermediate  coupling 
and  shell  model  calculations  in  fig.  22  and  with  the  predictions  of  the  T.T. 
model  in  fig.  24. 

Agreement  is  very  good  for  all  the  levels  below  2.0  MeV  except  for  the 
1.537  MeV  level  which,  because  of  its  predominant  hole-state  configuration, 
is  not  predicted  in  these  calculations.  One  common  feature  is  the  prediction 
of  a  level  with  spin  J  =5/2”  at  about  1.9  MeV.  This  would  indicate  that  the 
1.971  MeV  level  might  be  the  corresponding  level  with  J  =  5/2.  In  their 
calculations  using  the  T.T.  model,  Philpott  and  True  have  tried  different 
fits  based  on  the  spin  value  of  the  1.971  MeV  level.  They  concluded  that  the 
overall  results  favored  a  J  =  5/2  assignment  for  this  level.  Their  results 
presented  here  are  those  which  were  fitted  for  J  =  5/2. 

Above  2  MeV,  Ohnuma ' s  calculations  and  the  new  intermediate  coupling 
model  calculations  apparently  yield  a  better  density  of  states  than  the 
other  calculations.  There  are  indications  that  the  2.169  MeV  state  has  a 
spin  J  =  11/2  (Gu  70).  This  level  could  be  the  level  at  ti2.38  MeV  predicted 
by  Ohnuma' s  shell  model  calculations  and  the  level  at  ^  2.9  MeV  in  this  work. 
By  comparing  with  the  results  of  55Fe  (see  next  section),  the  2.227  MeV  level 
fahich  has  been  assigned  J  =  9/2  or  5/2  in  this  experiment),  could  have  the 
same  configuration  as  the  2.211  MeV  (9/2”)  level  of  55Fe,  as  it  lies  at  about 
the  same  excitation  energy  and  decays  100%  to  the  7/2  hole  state.  Therefore, 
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Figure  22  Energy  levels  and  spins  in  53Cr  as  predicted  by  intermediate 
coupling  model  and  shell  model  calculations  are  compared  with 
experimental  results. 
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Figure  23  Energy  levels  and  spins  of  55Fe  as  predicted  by  intermediate 
coupling  model  and  shell  model  calculations  are  compared  with 
experimental  results. 
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this  level  should  not  be  predicted  by  the  calculations.  There  is  a  9/2” 
state  predicted  around  2.3  MeV  in  the  different  calculations.  Again,  by 
comparing  with  the  results  obtained  for  55Fe,  it  can  be  seen  that  the  2.454 
MeV  level  of  53Cr  might  be  the  predicted  level  with  J  =  9/2.  The  study  of 
branching  ratios  in  section  6.2.  will  give  some  weight  to  this  assumption. 

r  r 

6.1.2.  Fe  spectrum 

The  experimental  spectrum  of  55Fe  is  compared  with  intermediate  coupling 
and  shell  model  calculations  in  fig.  23  and  with  the  predictions  of  the  T.T. 
model  in  fig.  24. 

The  same  general  remarks  apply  as  in  the  case  of  53Cr.  Both  shell  model 
calculations  of  Ohnuma  and  the  new  intermediate  coupling  model  calculations 
favor  a  spin  assignment  J  =  11/2  for  the  2.542  MeV  state. 

In  terms  of  density  of  states,  the  shell  model  calculations  of  Ohnuma 
and  the  present  intermediate  coupling  calculations  are  the  most  successful. 
The  former  calculations,  however,  seem  to  give  a  better  overall  fit  to  the 
experimental  spectrum  than  the  latter.  Ohnuma ' s  and  Larner's  calculations 
predict  a  J  =  1/2  level  at  ct  2  MeV.  This  level  might  correspond  to  the 
1.919  MeV  state  as  the  neutron  capture  study  of  Earle  and  Bartholomew  (Ea  66) 
favors  a  spin  3/2  for  the  2.052  MeV.  On  the  other  hand,  the  intermediate 
coupling  model  calculations  would  rather  favor  spin  assignments  J  =  3/2 


for  the  1.919  MeV  and  2.052  MeV  levels. 
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Figure  24  The  T.T.  model:  Calculations  of  energy  levels  in  53Cr  (Ph  70) 
and  55Fe  (La  70)  are  compared  with  experimental  results 
obtained  in  this  work. 
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6.2.  Branching  Ratios 

6.2.1.  5 3Cr 

Some  branching  ratios  measured  in  this  work  on  53Cr  are  compared  with 
theoretical  predictions  based  on  the  shell  model  (Oh  70)  and  on  the  inter¬ 
mediate  coupling  model  calculations  in  Table  X. 

The  agreement  is  generally  good  and  it  can  be  seen  that  the  calculated 
branching  ratios  are  consistent  with  a  spin  assignment  J  =  5/2  to  the  1.971 
MeV  state.  It  can  also  be  seen  that  the  assumption  that  the  2.454  MeV  level 
might  be  assigned  J  =  9/2  is  consistent  with  the  measured  branching  ratio. 

6.2.2.  55Fe 

In  Table  XI,  the  branching  ratios  for  some  levels  of  55Fe  measured  in 
this  work  or  reported  by  Pilt  (Pi  69)  are  compared  with  Ohnuma's  shell 
model  calculations  and  with  the  present  intermediate  coupling  model  calcu¬ 
lations  . 

The  overall  agreement  is  qualitative.  In  particular,  the  present 
calculations  would  rather  favor  a  J  =  1/2  assignment  rather  than  J  =  3/2 
for  the  2.052  MeV  level,  in  contradiction  with  the  measurements  of  Earle 
and  Bartholomew  (Ea  66).  This  contradiction  may  just  be  fortuitous,  as  it 
was  noticed  that  the  wave  functions  of  the  two  J  =  3/2  levels  at  —2.0  MeV 
excitation  energy  were  extremely  sensitive  to  slight  changes  of  K.  This 
is  caused  by  exchange  of  the  main  configurations  of  these  two  J  =  3/2 
levels  when  £  —2.75  (see  fig.  21). 


. 
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Table  X 

Calculated  and  measured  branching  ratios  for  some  levels  of  53Cr 
Calculations  are  based  on  the  intermediate  coupling  model  (this  work) 


and  the  shell  model  (Oh  70) . 


Initial  State 

Final  State 

Experiment 
(this  work) 

Theory 

(this  work) 

Shell  Model 
(Oh  70) 

1.006 

5/2" 

0 

100 

100 

100 

1.285 

7/2" 

0 

91  +  6 

55 

40 

1.006 

9  +  4 

45 

60 

1.971 

(5/2") 

0 

85 

87 

63 

1.006 

6 

14 

1.285 

15 

7 

23 

2.324 

3/2" 

0 

100 

50 

0.564 

25 

1.006 

23 

2.454 

0 

6+ 

15+ 

0.564 

<1 

1.006 

60 

26  87++ 

64 

1.285 

40 

68  13 

21 

3.625 

1/2"' 

0 

100 

94 

0.564 

3 

1.006 

3 

Assuming  J  =  5/2  for  the  2.454  MeV  level 
Assuming  J  =  9/2  for  the  same  level 


, 

Table  XI 
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Experimental  and  theoretical  branching  ratios  for  levels  of  55Fe  up  to 
2.578  MeV.  The  theoretical  calculations  are  based  on  the  intermediate 
coupling  model  (this  work)  and  on  the  shell  model  (Oh  70) 


Experiment  Theory 


Initial  State 

Residual  State 

(See  Chapter  III) 

This 

work 

Shell  Model 
(Oh  70) 

0.930 

5/2" 

0. 

98  +  1 

,w 

LOO 

83 

0.410 

2  +  i 

17 

1.318 

7/2“ 

0. 

96  +  1 

65 

33 

0.930 

4  +  1 

35 

67 

1.919 

(1/2, 3/2) 

0. 

68  +  3 

57f 

59++ 

40  + 

0.410 

32  +  3 

43 

22 

60 

0.930 

- 

19 

2.052 

(1/2, 3/2) 

0. 

23  +  2 

56* * 

37** 

0.410 

77  +  2 

44 

72 

59 

0.930 

- 

- 

27 

4 

2.144 

5/2" 

0. 

18  +  2 

41 

22 

0.410 

3  +  1 

4 

- 

0.930 

43  +  4 

53 

24 

1.318 

36  +  5 

2 

54 

2.301 

9/2" 

0.930 

80 

86 

-100 

1.318 

20 

14 

2.470 

3/2" 

0. 

100 

2.542 

(11/2") 

2.301 

— 

29 

1.318 

100 

71 

2.578 

5/2' 

0. 

84  +  2 

1 

0.410 

7  +  2 

5 

0.930 

6  +  2 

28 

1.318 

3  +  2 

66 

assuming  1/2"  for  1.919  MeV  level 

^  assuming  3/2  for  1.919  MeV  level 

*  assuming  1/2"  for  2.052  MeV  level 
__ 

assuming  3/2  for  2.052  MeV  level 
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6.3.  Electromagnetic  Properties 

6.3.1.  Reduced  transition  probabilities 

The  calculated  reduced  transition  probabilities  and  lifetimes  are 
shown  in  Table  XII  and  XIII  for  53Cr  and  55Fe,  respectively.  The  B(E2) 's 
were  calculated  for  neutron  charge  en  =  0  and  en  =  1.  It  can  be  seen  that 
the  B(E2)'s  are  not  too  sensitive  to  the  effective  charge  parameter.  This 
is  to  be  expected  as  most  of  the  contribution  to  the  electric  quadrupole 
transition  comes  from  the  core. 

In  Table  XIV,  the  B(E2)’s  for  transition  to  ground  of  the  lowest 
states  of  53Cr,  are  compared  with  the  experimental  data  of  Meriwether 
(Me  66)  and  Galparin  (Ga  69).  From  the  weak  coupling  model,  one  would 
expect  the  B(E2)  transitions  in  the  odd-nucleus  to  be  equal  to  the  B(E2; 

2+  -*■  0)  for  the  decay  to  ground  of  the  first  2+  state  of  the  neighbouring 
even  nucleus.  In  the  case  of  52Cr,  the  B(E2;  2+  -*  0)  is  in  the  order  of 
50  to  140  e2-fmL<,  i.e.  6  to  16  w.u.  (Ri  69,  Me  66);  in  54Fe,  the  measured 
B(E2;  2  +  0)  is  of  the  order  of  cr  100  e2-fm4,  i.e.  —  11  w.u.  (Fu  64,  Si  65). 
In  the  simple  picture  of  the  intermediate  |2  ,P3/2>  coupling,  the  ratio  of 
B(E2)  for  decay  to  ground  of  levels  of  the  |2+,P 3/  >  quadruplet,  to  the 
B(E2;  2+  +  0)  in  the  even  nucleus  is  expected  to  increase  from  unity  for 
those  states  with  spin  1/2,  5/2  and  7/2  and  to  decrease  from  unity  for  the 
3/2  state  as  the  coupling  increases.  Results  in  Table  XIV  show  that  this 
simple  picture  does  not  hold  for  53Cr.  The  main  reason  is  probably  because 
of  the  mixing  of  states  in  terms  of  single  particle  configurations  and  of 
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Table  XII 

Electromagnetic  transition  strengths  and  lifetimes  in  53Cr  calculated  from 
the  intermediate  coupling  model.  The  lifetimes,  calculated  for  en  =  0, 


are  compared  with  experimental  results  when  available » 


Initial 

Residual 

B(E2) 

(e2-fm4) 

B  (Ml) 

Lifetime  x 

State 

State 

en  =  0 

en  -e 

(eti/2Mc) 2 

calc . 

experiment 

0.564 

1/2 

0. 

3/2 

101 

179 

1.97 

0.16  ps 

>  0.50  psa) b) c) 

1.006 

5/2 

0. 

3/2 

54.1 

74.8 

0.009 

4.3  ps 

i  a) 

>  1  ps 

1.285 

7/2 

0. 

3/2 

112 

104 

forbidden 

1.1  ps 

1.971 

(5/2) 

0. 

3/2 

36.8 

29 

0.031 

94  f s 

1.285 

7/2 

44.0 

81.8 

0.011 

' (4+2)  fsb^ 

^ (6.2+2)  fsC; 

2.324 

3/2 

0. 

0.564 

3/2 

1/2 

50.3 

0.002 

32.9 

0.024 

0.073 

0.085 

24  fs 

1.006 

5/2 

19.9 

35.0 

0.236 

[2 . 190] * 

_  /  _ "k  * 

3/2 

0. 

3/2 

0.002 

0.0001 

0.0002 

2.454 

(9/2) 

1.006 

5/2 

159 

215 

forbidden 

0.66  ps 

1.285 

7/2 

17.1 

25.9 

0.0047 

[2.90]* 

11/2 

1.285 

7/2 

193 

237 

forbidden 

0.25  ps 

2.454 

9/2 

0.182 

0.179 

1.13 

Only  negative  parity  states  have  been  considered 
calculated  energy 

third  3/2  state  in  the  calculations 


Ref .  a)  This  work 

b)  Bo  63 

c)  Ba  60 


Table  XIII 
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Electromagnetic  transition  strengths  and  lifetimes  in  55Fe  calculated  from 


the  intermediate  coupling  model.  The  lifetimes,  calculated  for  en  =  0, 
are  compared  with  the  recent  measurements  of  Robertson  (Ro  70)^ 


Initial 

Residual 

B(E2) 

(e2-fm4) 

B  (Ml) 

Lifetime  t 

State"*"  * 

State 

en  =  0 

en  =e 

(eh/  2MC)2 

Theory 

Experiment 

0.410 

1/2 

0. 

3/2 

130 

210 

1.89 

0.43  ps 

0.930 

5/2 

0. 

3/2 

67.9 

90.4 

0.014 

3.8  ps 

1.318 

7/2 

0. 

3/2 

129 

161 

forbidden 

0.56  ps 

0.930 

5/2 

0,287 

0.40 

1.10 

1.919 

(1/2) 2 

0. 

3/2 

17,9 

1.46 

0.029 

A 

110  fs 

0.410 

1/2 

forbidden 

0.051 

)  <25  «)  fs 

1.919 

(3/2) 2 

0. 

3/2 

44.5 

27.4 

0.096 

44  fs 

0.410 

1/2 

0.0078 

0.0040 

0.085 

2.052 

(1/2) 2 

0. 

3/2 

17.9 

1.46 

0.029 

110  fs 

0.410 

1/2 

forbidden 

0.051 

> (31  +  8)  fs 

2.052 

(3/2)3 

0. 

3/2 

0.020 

0.043 

0.0008 

76  fs 

0.410 

1/2 

82.6 

87.5 

0.105 

2.144 

5/2 

0. 

3/2 

34.1 

27.4 

0.039 

0.410 

1/2 

46.9 

59,8 

forbidden 

47  fs 

55  +16  fs 

0.930 

5/2 

21.3 

35.5 

0.351 

^  -12  S 

1.318 

7/2 

57.1 

96.5 

0.040 

2.301 

9/2 

0.930 

1.318 

5/2 

7/2 

180 

21.4 

245 

30 . 7 

forbidden 

0.0089 

0.8  ps 

„  „  +0.7 
°‘9  -0.23PS 

2.470 

3/2 

19  +  10  fs 

2.542 

(IV  2) 

1.318 

7/2 

210 

260 

forbidden 

1  ps 

>  0.66  ps 

2.301 

9/2 

0.537 

0,633 

1.20 

2.578 

5/23 

0. 

3/2 

1.82 

35.3 

,0002 

0.410 

1/2 

27.3 

17.1 

forbidden 

31  fs 

(67  +  9)  fs 

0.930 

5/2 

6.37 

6.71 

0.116 

1.318 

7/2 

0.987 

1.25 

0.438 

f  All 

states 

have  negative  parity 

++  The 

lower 

index  i 

represents 

the  i*-k 

state  of 

the  same  spin 

as  obtained 

from*  theory . 

It 

nas  been  omitted 

for  the  first  states. 

' 
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Table  XIV 

B(E2)  values  in  units  of  for  decay  to  ground  of  lower 

states  of  53Cr.  Experimental  values  are  based  on  (a, a') 


experiment  (Me  68)  and  coulomb  excitation  (Ga  69) 


Experiment 

Theory 

Energy 

(MeV) 

(Me  66) 

(Ga  69) 

This 

e  =0 
n 

work 

e  =6 

T.T.  model 
(Ph  70)  b> 

Shell  model 
(Oh  70) °) 

0.564 

178 

220 

101 

197 

146 

40 

1.006 

35 

14 

54 

75 

35 

22 

1.285 

69 

120 

112 

139 

63 

117 

1.971 

33 

29 

62.1 

47 

a)  Intermediate  coupling  model  calculation 

b)  eR  =  1.64  e 

c)  ep  =  2e ;  en  =  e 
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Table  XV 


Calculated  and  measured  mixing  ratios  in  5jCr  and  55Fe 
Calculations  were  obtained  from  the  intermediate  coupling  model 


Nucleus 

Transition 

Experiment 

Theory 

53Cr 

1.006  +  0 

-0.27+°’^  a) 
-o.31 

0.67 

2.324  -*  0 

-5.7  <  6  <  0.09  a) 
0.47  £  6  1  5.7 

-0.51 

5  5Fe 

0.930  ->  0 

0.36  +  0.11  b) c) 

0.53 

2.144  ->  0 

4*  00  c.  ^ 

lo°  -1.3 

0.52 

2.470  -*  0 

-1.0  +1-°  c> 

—  nc 

0.12 

a)  this  work 

b)  Ba  60 

c)  Pi  69 


. 
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more  than  one  phonon  component.  The  intermediate  coupling  model  calculations 
show  that  qualitative  agreement  is  obtained.  In  the  T.T.  calculations,  good 
agreement  is  expected  as  there  are  two  free  parameters  to  fit  three  pieces 
of  data  (Ph  70) .  In  the  shell  model  calculations  of  Ohnuma,  the  agreement  is 
poor  although  an  effective  charge  of  lewas  taken  for  neutron  and  proton. 

This  tends  to  reinforce  the  confidence  in  the  collective  description.  All 
calculations  predict  a  non-negligible  B(E2)  for  the  1.971  MeV  state.  Although 
this  state  is  well  excited  in  inelastic  scattering  experiments  (see  e.g. 

(Me  66)  and  Wh  67))  no  B(E2)  measurement  has  been  reported. 

6.3.2.  Mixing  ratios 

Very  few  mixing  ratios  have  been  accurately  measured  and  comparison 
between  theory  and  experiment  can  only  be  semi-quantitative.  Results  are 
shown  in  Table  XII,  where  the  experimental  mixing  ratios  of  some  transitions 
in  53Cr  and  55Fe  are  compared  with  the  predictions  of  the  intermediate 
coupling  model.  The  same  phase  convention  is  used.  The  overall  agreement 
is  good  except  for  the  decay  to  ground  of  the  1.006  MeV  (5/2-)  state  of  53Cr, 

for  which  there  is  apparently  a  contradiction  in  sign  between  experimental 

\ 

and  theoretical  results.  It  must  be  noted  that  there  is  agreement  between 
theory  and  experiment  for  the  mixing  ratio  of  the  corresponding  transition 
in  the  0.930  MeV  (5/2”)  state  of  55Fe. 

6.3.3.  Lifetimes 

Tables  XII  and  XIII  also  show  the  calculated  measured  lifetimes.  The 
overall  agreement  is  good,  mainly  for  the  55Fe  results.  However,  for  some 
levels  like  the  2.324  MeV  (3/2")  in  53Cr  and  the  1.919  MeV  (1/2,  3/2), 


- 

' 

tC’. 
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2.052  MeV  (1/2,  3/2)  and  the  2.578  MeV  (5/2)  in  55Fe, disagreement  between 
experimental  and  theoretical  branching  ratios  indicates  that  the  correspond¬ 
ing  theoretical  lifetimes  should  be  considered  with  caution.  For  the  other 
levels,  it  can  be  concluded  that  the  calculated  wave  functions  give  a  good 
description  of  their  configuration. 

6.3.4.  Static  moments 

Static  moments  for  53Cr  and  55Fe  calculated  from  different  models  are 
shown  in  Table  XVI.  Only  for  53Cr  are  experimental  values  available. 


Table  XVI 

Calculated  and  experimental  static  moments  in  53Cr  and  55Fe 


Electric 

quadrupole  moment  Q 

(in  barns) 

Nucleus 

Exp . 

(Fu  69) 

This  work 

T.T.  model 
(Ph  70) 

Shell  model 
(Oh  70) 

Single 

neutron 

state^ 

53Cr 

+  0.03 

-  0.09 

-  0.15 

-  0.17f 

-  0.04 

55Fe 

-0.125 

-  0.16f 

-  0.04 

i 

calculated  for 

tf 

calculated  from 

p  n 

e  ■  =  e 

eff  eff 

formula  (V 

»  1 

.1.)  (Bo  53) 

Magnetic 

dipole  moment  y  (in  n.m.) 

Nucleus 

Exp . 

(A1  53) 

This  work 

Shell  model 
unquenched 

(Oh  70) 
quenched 

Schmidt 

limit 

53Cr 

-  0.474 

-  1.83 

-  1.12 

-  0.47 

-  1.91 

55Fe 

-  1.62 

-  1.03 

-  0.40 

-  1.91 

■ 


■ 

. 
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The  first  published  measurement  of  the  electric  quadrupole  moment  of 

r  q 

Cr  (Te  61)  corresponded  to  a  value  Q  =  -0.03b.  Since  then,  there  have  been 
many  contradictory  results,  and  at  present  the  value  adopted  by  Fulmer  and 
Cohen  (Fu  69)  is  |q|  =  0.03b  .  It  is  interesting  to  note  that  all  the 
calculations  indicate  that  the  sign  should  be  negative.  All  theoretical 
calculations  yield  a  value  too  large  in  magnitude  when  compared  with  the 
experimental  result,  except  for  the  single  particle  calculation.  The  only 
explanation  could  be  either  that  the  ground  state  wave  function  obtained 
from  the  different  calculations  does  not  give  a  good  description  of  the 
ground  state,  or  that  for  some  unknown  reason,  some  terms  should  cancel 
out.  The  situation  is  identical  when  one  compares  the  experimental  value 
of  the  magnetic  dipole  moment  of  ^3Cr  to  the  theoretical  estimates.  However, 
in  this  case  the  single  particle  estimate  (i.e.  the  Schmidt  limit)  is  about 
four  times  larger  than  the  experimental  value.  It  appears  that  the  best 
value  is  obtained  when  "quenching"  effects  are  taken  into  account.  The 
theories  which  try  to  explain  why  the  anomalous  magnetic  moment  of  the  odd 
nucleon  is  somewhat  "quenched"  in  the  presence  of  other  nucleons  usually 
involve  consideration  of  the  mesonic  field  (see  e.g.  (Pr  62)  p.  74,  and 
references  therein) .  This  might  explain  why  all  the  models  fail  to  describe 
the  magnetic  dipole  moment. 

6.4.  Spectroscopic  Factors 

In  Tables  XVII  and  XVIII,  experimental  spectroscopic  factors  obtained 
in  (d,p)  stripping  reactions  are  compared  with  values  calculated  from 


vswoi  .es  J&ntiaas  j  o3  l1  ^  io  Jr;9ui  ,<  sloqlfc  alJW- >  ^d3  j 
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different  models,  including  the  present  intermediate  coupling  model  calcu¬ 
lations.  From  these  results; it  appears  that  most  of  the  calculations  yield 
a  fair  agreement  with  the  experimental  results.  A  point  of  interest  is  to 
compare  the  two  intermediate  coupling  calculations, bearing  in  mind  that  one 
of  Ramavataram’ s  main  goals  was  to  check  his, against  the  measured  spectro¬ 
scopic  factors.  For  the  low  lying  states  the  results  are  almost  equivalent 
except  for  the  single  particle  strength  of  the  first  5/2  level.  Ramavataram ' s 
better  results  for  this  level  suggest  that  the  (f5 ^  — p 3 ^ 2 )  single  particle 
energy  might  be  smaller  than  the  value  used  in  the  present  calculations .  On 
the  other  hand,  these  calculations  seem  to  describe  better  than  Ramavataram’ s 
the  configuration  of  the  2.324  MeV.  In  the  case  of  the  T.T.  model,  the 
results  are  not  so  good  and  very  likely  imply  that,  because  only  two  core 
states  are  taken  into  consideration,  the  odd-A  nuclei  wave  functions  are  too 
limited . 


, 


Sample  spectroscopic  factors,  (2J+1)S,  for  the  52Cr (d ,p) 5 3Cr  reaction 
The  experimental  values  were  obtained  by  averaging  the  results  given  in  (An  64,  Bo  65, 

SS  64  and  Ra  68)  and  are  compared  with  theoretical  values  calculated  from  different  models. 
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Table  XVIII 

Sample  spectroscopic  factor  (2J+1)S  for  the  54Fe(d,p) 55Fe  reaction 
The  experimental  values  were  obtained  by  averaging  the  results  given  in  (Ma  64)  and 
(Fu  63)  and  are  compared  with  theoretical  values  calculated  from  different  models. 


121 


'X) 

vO 

1 - 1 

0)  CL) 
T3  > 
O 

E 


i — I  vO 

a)  vo 
A 

CO 

O 


<U 

03 

O 

E 


o 

h-J 


H 

H 


i— I 
0) 

O^, 

Ew 

bo  co 
(2  co 
•H 

i — I 

P.  P3 
d 
o 
a 

Q) 

•h  ^ 

d  Ci 
•H  O 
03  > 

a> 

E  co 
M  *H 
0)  42 
+J  H 
(3 


H> 

a 

a) 

E 

•H 

a) 

P 

X 

w 


t= 

•“5 


bO  > 
s-i  CU 
a>  g 
(3  ^ 
W 


00 

o 

00 

1 - 1 

vO 

O 

o 

CM 

CM 

oo 

o 

UO 

CO 

vO 

CO 

o 

o 

CO 

o 

o 

o 

O 

o 

o 

o 

CM 

> 

CM 

00 

CM 

o 

CO 

UO 

00 

O 

O'. 

O'! 

CM 

uo 

Co 

CO 

CM 

• 

® 

• 

• 

• 

• 

• 

• 

CO 

o 

CO 

o 

o 

O 

o 

o 

CM 

1 - 1 

w 

Kj- 

CM 

M3 

vO 

CO 

CM 

rH 

vO 

U0 

vO 

CO 

• 

• 

• 

• 

• 

• 

CO 

1 - 1 

CO 

O 

o 

1 — 1 

CM 

CM 

CO 

CO 

w 

o 

1 - 1 

CO 

Co 

uo 

uo 

oO 

vO 

CO 

1 - 1 

CM 

« 

• 

• 

• 

• 

• 

CM 

o 

CM 

o 

O 

o 

✓“—s 

CM 

CO 

s~' 

00 

CO 

CO 

CM 

CM 

vO 

CD 

lO 

O 

CM 

00 

CM 

00 

CM 

o 

1 - 1 

O 

o 

O 

1 — 1 

O 

CO 

uo 

o 

CM 

vO 

CM 

CM 

c- 

• 

1 — 1 

1 - 1 

rH 

1 - 1 

O 

O 

• 

o 

• 

• 

• 

• 

® 

• 

• 

o. 

o, 

o, 

O, 

o, 

o, 

o, 

o, 

+ 1 

+  1 

+  1 

+  1 

+  1 

+ 

+ 1 

+ 

+ 

UO 

uo 

uo 

\D 

rH 

CM 

o 

UO 

ON 

00 

1 - 1 

CM 

UO 

CM 

o 

CM 

o 

CM 

o 

O 

o 

O 

o 

o 

1 

1 

1 

1 

/'“'V 

1 

1 

1 

1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

rH 

uo 

CO 

UO 

CO 

U0 

CO 

y—s 

tH 

i — 1 

CO 

rH 

rH 

CO 

1 - 1 

CO 

'w' 

o 

O 

Co 

CM 

o 

00 

rH 

CO 

rH 

U0 

O' 

CM 

Co 

CO 

o 

i — i 

•<r 

U0 

o 

• 

o 

o 

O 

1 - 1 

CM 

CM 

CM 

CM 

CO 

co 

bO 


*)  Normalized  at  2.50  for 


■ 

' 

. 

. 


122 


7.  CONCLUDING  REMARKS 

Considering  the  simplicity  of  the  intermediate  coupling  model,  it  appears 
that  theoretical  predictions  are  in  good  agreement  with  experiment.  The 
model  reproduces  well  most  of  the  properties  of  the  levels  below  2.5  MeV 
excitation,  except  the  static  moments,  which  all  models  discussed  here  fail 
to  reproduce.  In  most  cases,  the  contributions  from  coupling  to  3-phonon 
states  are  negligible  but  contributions  from  2-phonon  states  must  be  taken 
into  account.  As  a  result,  the  overall  agreement  obtained  from  these  cal¬ 
culations  is  better  than  those  based  on  the  T.T.  model  in  which  only  the 
coupling  to  ground  and  first  excited  state  of  the  core  is  considered.  It 
seems  that  if  one  wants  to  keep  the  advantages  of  the  more  general  approach 
suggested  by  Thankappan  and  True,  one  should  include  contributions  from 
higher  excited  states  in  the  core,  while  trying  at  the  same  time  to  limit 
the  number  of  free  parameters  by  assuming  that  "some"  collective  features 
are  the  same  for  different  states  of  the  core.  This  should  be  a  reasonable 
assumption  in  most  cases. 

There  are  many  common  points  between  intermediate  coupling  and  shell 
model  calculations.  This  suggests  that  the  form  of  interaction  taken  in 
the  unified  model  is  valid.  Each  model,  however,  ignores  some  levels  of 
the  excited  core:  in  shell  model  calculations,  the  0+  2-phonon  vibrational 
state  is  not  taken  into  account,  in  the  intermediate  coupling  model,  a  level 
like  the  2.37  MeV  (4+)  state  in  52Cr  is  ignored  as  it  is  not  explained  by 
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the  vibrational  model.  It  is  not  possible  for  the  time  being  to  evaluate 
the  importance  of  these  neglected  terms.  Both  models  also  fail  to  reproduce 
the  hole-states  such  as  the  second  7 / 2 —  states  in  53Cr  and  55Fe.  It  must  be 
pointed  out,  however,  that  these  model  calculations  presented  here  can  be 
qualified  as  "simple".  By  taking  more  particles  into  account  and  using  the 
"particle-hole"  formalism,  undoubtedly  one  should  get  better  results,  but 
then  the  problem  lies  in  prohibitive  matrix  sizes  and  computation  difficult¬ 
ies  . 

For  the  time  being  the  important  question  is  to  try  and  estimate  how 
the  overall  results  would  be  modified  when  hole  states  are  included.  From 
the  study  of  pick-up  reaction,  which  shows  that  the  cross  section  for  say, 
the  7/2  hole  state  is  much  larger  than  the  pick-up  cross  section  for  the 
other  7/2  state,  and  from  the  measured  branching  ratios  which  show  that 
within  experimental  error  a  level  like  the  9/2”  state  does  decay  to  either 
7/2  states  but  not  to  both,  it  seems  to  be  possible  to  assume  that  mixing 
between  such  states  is  probably  small. 

If  one  accepts  the  validity  of  the  type  of  interaction  used  in  the 
intermediate  coupling  model  calculations,  it  appears  that  there  are  possi¬ 
bilities  of  improving  the  calculations.  Rather  than  considering  fmo  as  a 
free  parameter,  it  would  be  better  to  couple  the  odd-nucleon  to  the  known 
vibrational  levels  of  the  core:  for  instance,  in  53Cr  the  1.434  MeV  (2+) 
would  be  the  N  =  1  phonon  state,  the  2.65  MeV  (0+)  would  be  the  0+  N  =  2 
phonon  state,  etc.  This  assumes  that  the  general  theory  is  still  valid 
and  that  the  odd-particle  does  not  appreciably  affect  the  core  properties. 


, 
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The  overall  results  might  not  be  greatly  modified  but,  with  the  elimination 
of  one  free  parameter,  it  should  then  be  possible  to  get  a  better  estimate 
of  the  single-particle  energy  spacings .  Another  possible  improvement  should 
result  from  the  use  of  effective  charge  and  effective  gyromagnetic  parameters 
in  the  calculation  of  electromagnetic  transitions.  These  parameters  could 
be  determined  by  fitting  either  B(E2)'s  or  the  measured  static  moments  but 
this  requires  that  the  experimental  values  be  accurately  known  to  make  such 
a  procedure  worthwhile. 

No  mention  has  been  made  of  the  positive  parity  states.  In  53Cr  and 
53Fe,  there  is  apparently  no  positive  parity  state  below  3.5  MeV  excitation. 
This  is  above  the  excitation  which  has  been  considered  as  the  limit  of 
validity  of  the  present  model.  Therefore,  positive  parity  states  can  be 
ignored  in  these  calculations. 


' 


, 


CHAPTER  V 


CONCLUSIONS 


In  the  study  of  53Cr,  some  improvement  has  been  brought  to  Litherland 
and  Ferguson  Method  II,  by  limiting  the  ratio  of  population  parameters.  A 
better  knowledge  of  the  reaction  mechanism  might  result  in  more  restrictive 
conditions  of  this  ratio.  This  would  then  limit  the  number  of  possible  spin 
assignments  and  might  also  yield  more  accurate  mixing  ratios.  Experimentally, 
this  method  requires  rather  sophisticated  electronic  set-ups.  On  the  other 
hand,  if  the  bombarding  energy  is  well  chosen,  many  levels  can  be  studied 
at  the  same  time  (the  only  limiting  factor  being  the  number  of  "bins" 
available) .  As  the  gamma  ray  spectra  are  recorded  separately  by  sorting 
according  to  the  energy  of  the  different  proton  groups,  branching  ratios  can 
be  accurately  determined  for  most  of  the  states  except  in  the  case  of  closely 
spaced  multiplets.  High-spin  levels  which  are  non-stripping  states  are 
weakly  excited  in  the  (d,py)  reaction  and  another  type  of  reaction  must  be 
used  to  investigate  such  levels. 

In  the  study  of  55Fe  via  the  55Mn(p,ny) 55Fe  reaction,  the  experimental 
method  is  quite  appealing  by  its  simplicity.  This  reaction  also  proves  to 
be  quite  adequate  for  the  investigation  of  high-spin  states.  On  the  other 
hand,  levels  with  spin  J  <_  5/2  are  not  easily  assigned  with  this  particular 
reaction.  As  the  population  parameters  are  calculated  from  the  statistical 
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model,  there  is  one  less  free  parameter  in  the  x2-fitting  calculations  and 
one  then  should  expect  a  more  accurate  determination  of  the  multipole  mixing 
ratios.  With  a  suitable  choice  of  bombarding  energies,  many  levels  can  be 
investigated  at  the  same  time.  This  is  important,  as  a  careless  application 
of  this  method  could  be  quite  machine-time  consuming.  In  any  case,  it 
appears  that  the  number  of  levels  which  can  be  investigated  by  this  method 
is  limited.  When  too  many  levels  are  being  excited,  the  resulting  gamma 
spectrum  becomes  quite  complex  and  some  transitions  are  masked  by  other 
gamma  rays.  The  determination  of  decay  modes  also  becomes  difficult  because 
of  multiple  feeding  of  the  lower  states. 

In  55Fe,  it  was  hoped  to  find  the  13/2  level  predicted  by  the  calcula¬ 
tions  of  Ohnuma  at  2.975  MeV.  No  evidence  was  found  of  the  existence  of 
such  a  level  although  one  could  speculate  about  the  unobserved  (?)  2.818  MeV 
level  being  such  a  state. 

Theoretically,  the  failure  of  Nilsson  model  calculations  confirmed  that 
52Cr  and  5uFe  cannot  be  regarded  as  strongly  deformed  nuclei.  Thankappan 
and  True- type  calculations  are  quite  appealing  by  their  general  approach. 
However,  they  appear  to  have  a  rather  limited  range  of  application  and  could 
only  be  improved  by  including  more  core  excited  states. 

For  the  present  state  of  knowledge,  intermediate  coupling  and  shell 
model  calculations  have  proved  to  be  successful  in  describing  most  of  the 
experimental  features  reported.  As  mentioned  earlier,  both  types  of  calcu¬ 
lations  could  be  refined  and  this  may  be  necessary  when  higher  excited  states 
are  considered.  The  failure  to  reproduce  the  electric  quadrupole  moment  of 
53Cr  remains  unexplained.  Measurements  of  the  electric  quadrupole  moment  of 
the  other  N  =  29  nuclei  may  give  some  clue  to  the  problem. 
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Abstract 

Two  programs  based  on  the  intermediate  coupling  in  the  unified 
model  are  described.  The  first  program,  CORPAR  I,  calculates  energy 
levels  in  the  framework  of  this  model  for  a  single  particle  having 
up  to  three  orbitals  available,  coupled  to  a  vibrating  core  where 
up  to  three  phonon  state  vibrations  are  considered.  The  second 
program, CORPAR  1 1,  cal cul ates  electromagnetic  transitions,  mixing 
ratios,  lifetimes  and  nuclear  static  moments  using  the  wave  functions 
calculated  by  CORPAR  I. 
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CHAPTER  I 
CORPAR  I 


Title:  CORPAR  I 

Type:  Main 

Source  Language:  Fortran  I V (G ) 

Origin:  T.P.G.  Carola,  Nuclear  Research  Centre 

Abstract:  This  program  calculates  the  total  Hamiltonian  HT  in  the  frame 

of  the  intermediate  coupling  model.  After  di  agonal  izati  on , 
the  energy  eigenvalues  and  wave  functions  are  given  for  each 
value  of  the  final  spin  state  I  and  the  coupling  parameter  £. 


1.  Description  and  theory 

1.1  Introduction 

The  intermediate  coupling  in  the  unified  model  was  first  outlined  by 
Bohr  (Bo  52)  and  Bohr  and  Mottelson  (Bo  53).  Later,  it  was  elaborated  by 
Choudhury  (Ch  54)  and  Ford  and  Levinson  (Fo  55).  This  model  has  since  been 
used  (see  e.g.  G1  60,  Ro  63,  Ch  67,  Ru  68  and  Pa  70)  to  describe  the  spectra 
and  other  relevant  physical  quantities  such  as  reduced  transition  proba¬ 
bilities,  lifetimes,  spectroscopic  factors  ...  in  the  case  of  odd  nuclei 
whose  even  mass  neighbouring  nuclei  exhibit  a  vibrational  spectrum. 

In  this  model,  the  odd  mass  nucleus  is  treated  as  a  coupled  system 
comprising:  i)  a  doubly  even  core  capable  of  performing  collective 
quadrupole  oscillations  of  frequency  go;  ii)  an  extra  core  nucleon  which 
has  several  single  particle  levels  available. 


- 
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In  the  program  described  here,  the  single  particle  has  up  to  3 
orbitals  available.  Core  excitation  up  to  3  phonons  are  considered. 

The  energies  of  the  resulting  states  of  the  coupled  system  depend 
upon  the  following  parameters: 

Hw  phonon  energy 

e 1  energy  difference  between  the  first  and  ground  state 

orbi tal 

e2  energy  difference  between  the  second  and  ground  state 

orbi tal 

€  coupling  parameter. 

Note : 

As  described  extensively  in  (Bo  53), the  limiting  cases  of  this  model 
correspond  either  to  weak  coupling  (for  nuclei  near  double  closed 
shells)  or  strong  coupling  (case  of  deformed  nuclei)  and  are 
usually  treated  differently. 

1.2  Formalism 

The  core-parti cle  system  is  described  by  the  coupling  of  a  single 
particle  of  angular  momentum  j  with  the  quadrupole  oscillations  of  the 
nuclear  surface  of  the  core. 

The  total  Hamiltonian  is  the  sum  of  three  terms: 


Hr 


H  (a„  )  +  H  (x)  +  H,  (a  ,x) 
s'  2y'  spv  int  2y 


(1) 


where:  a2)J  are  the  collective  coordinates  as  defined  by  Bohr  and  Mottel- 
son  (Bo  53) 

x  are  the  single  particle  coordinates. 

The  surface  Hamiltonian  Hs  is  defined  by: 


V-  :  ••i-mo'lt 


I  f 


4 


M“2y)  =  I  {r  B2  I  “2u  1 2  +  2-C2l%l2}  <2) 

in  terms  of  the  mass  parameter  E>2  and  the  nuclear  deformabi  li  ty  C2. 

It  is  equivalent  to  a  system  of  harmonic  oscillators  of  frequency 


The  particle  Hamiltonian  is  a  constant.  It  corresponds  to  the 
motion  of  the  odd  nucleon  in  the  effective  average  potential  of  the  spheri¬ 
cal  core.  The  correspond!' ng  eigenvalues  are  then  simply  the  single 
particle  energy  levels. 

Finally,  the  Hamiltonian  H  represents  the  coupling  of  the  single 
particle  to  the  vibrating  surface.  It  can  be  written  as: 

Hint  =  -  k  I  “2pV0-*)  <4> 

r* 

where  k  is  the  coupling  constant  introduced  by  Bohr  and  Mottelson  (Bo  53); 

Y  is  the  spherical  harmonic;  0  and  <j>  are  the  polar  angles  of  the  particle. 
It  is  more  convenient  to  write  the  surface  variable  in  terms  of 

creation  and  annihilation  operators  b*  and  b  .  The  expression  for  H 

y  l-i  int 

then  becomes: 

Hint<“2y>  -  -  k  Z  /if  <by  -  (->y  b*y>  Y  2y<9’*>  <5> 


* 

It  is  actually  the  radial  average  k  =  <k(r)>  and  can  be  considered 
as  a  constant. 


' 
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The  representation  of  the  system  is  in  the  angular  momentum  space.  More 
specifically,  the  eigenvectors  of  the  uncoupled  system  |jm>  and  |N  R  m'>, 
representing  the  single  particle  and  the  core  respectively,  are  taken  as 
the  basis  states,  j  is  the  angular  momentum  of  the  single  particle,  m 
its  projection,  N  is  the  number  of  phonons  of  angular  momentum  R. 

The  basis  vectors  of  the  coupled  system  are  then  given  by: 

|j  NR  IM>  =  l  ( j  R  mm' | IM) | j  m> |NR  m'>  (6) 

m,m' 

where  (j  R  mm'|IM)  is  the  Clebsch-Gordan  coefficient  resulting  from  the 

-> 

angular  momentum  coupling.  The  total  angular  momentum  I  =  j  +  R.  From 
this  definition,  it  follows  that:  (neglecting  the  zero-point  energy) 

(H  +  H  )|j  NR  IM>  =  (Nfia)  +  E-s )  j  j  NR  IM>  (7) 

s  p  -> 

with  Ej  correspondi ng  to  the  energy  of  the  single  particle  with  total 
angular  momentum  j.  Using  eq.  (6),  one  can  now  evaluate  the  matrix 
elements  of  Hi  t: 

<j'  N 1 R '  IM|Hint|j  NR  IM>  =  -  k  < j  1  N'R!  IM|J(b  +  (-)V  ) 

’  u 

x  Y  . (e  ,$)  |  j  NR  IM>  (8) 

z  y 

For  the  time  being,  we  will  write  -k 
of  equation  (8)  as  <Hint> . 


as  K  and  the  left  hand  term 


JWOfl]  IS 1 


■n  ■ 

•  ■  i 
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<Hint>  =  K  {<y  N'R'  IM-I  I  by  Y2y|j  NR  IM> 

y 

+  <j‘  N'R'  IM|E(-)V|  y,  |j  NR  IM> }  (9) 

y  y  y 

Consider  the  first  term  of  eq.  (9).  Using  eq.  (6),  one  gets 

<  II  by  Y2yl  =•  =  I  < J  1  v | Y 2  | jm> . <N 1 R ' v | b  | NR  n‘>(j'  R'  vv'|IM) 

y  mm' 

vy.  x  ( j  R  rrm '  ]  IM)  (10) 

-k 

One  can  now  apply  the  Wigner-Eckart  theorem,  bearing  in  mind  the  differ¬ 
ence  between  an  operator  acting  on  the  ket  (like  Y  )  and  the  bra 

Z  pi 

(like  b  )  vectors . 

<  II  by  V2y|  >  =  I  (j  2  m  y  |  j  '  v )  <  j  '  |  |  Y2  |  |  j  >  ( R 1  2v'y|Rni')<N,R,||b||NR> 

y  mm' 

y 

W  x  (j1  R'  vv1 | IM) ( j  R  mm' [ IM)  (11) 

What  follows  is  purely  Racah  algebra  (see  e,g.  Ro  57  -  pages  38  and 
110).  In  order  to  get  a  Racah  coefficient  out  of  the  sum  of  our  Clebsch- 
Gordan  coefficients,  some  rearranging  is  necessary: 

"  II  b  Y 2  |  >  =  l  I  <j,||V2||j><N'R'||b||NR>  (-)R'+2‘R(j  2  m  y|j'v) 

y  mm'  y 

x  (2  R'  y  v'  |  Rm' )  ( j  1  R'  w'  |  IM)  ( j  R  mm*  j  IM) 

02) 

* 

The  Wigner-Eckart  theorem  states  that: 

< j '  m ' | tlm I J  m>  =  ( j  L  m  M | j '  m ' )  < j '  |  | T2 1  I J> 
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This  can  then  be  more  simply  expressed  as: 

'  I  I  bu  Y2yl  >  =  (-)R'_R  <J"  I  |v  |  IjxN'R'  I  |b|  | NR>  /(2j'  +  l)(2R+l) 

M 


x  W ( j  2  I  R';j'R)  (13) 


<  I  I  by  Y2|  >  =  (-)*'  R+j'+R  I  y^TT  <j'  I  I Y2 1  |j>  yjR+r  <N ' R '  |  jb|  |NR> 

y 

x  W  ( j  j  '  RR';2I)  (14) 

The  treatment  of  the  second  term  of  eq.  (9),  is  equivalent  to  what 
has  just  been  done.  The  <N'R' | |b| |NR>  terms  contribute  only  when  N'  <  N, 
whereas  the  <N'R' | |b*| |NR>  terms  contribute  only  when  N 1  >  N .  As  the 
total  matrix  is  symmetric,  it  suffices  to  calculate  only  elements  on 
one  side  of  the  diagonal.  Limiting  to  the  case  N'<  N,  one  gets: 

<  I Hint  |  >  =  (-)R'+r'1  K  <j  1  |  |  Y2  |  |  j> 

x  /2R+T  <N ' R 1  I  I b I i NR>  W(jj'  RR 1 ; 21) 

for  N 1 <N  (15) 

One  can  now  introduce  the  dimensionless  coupling  parameter  defined  by 

5  =  k  /5/2ir.-ho).C  (16) 


(15)  then  becomes: 
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'  lHintl  >  =  (-)R,+3''I+1  yi  -no.  c  /2j'  +  l  <j  '  |  |  v2  |  |  j> 

x  /2R+1  <N ' R 1 | | b | | NR>  W  ( j  j  '  RR 1 ; 21 ) 

for  N'<N  (17) 

Some  authors  write  the  phase  factor  as  (-)R  +I.  Both  expressions 

obviously  yield  the  same  result. 

The  reduced  matrix  elements  for  Y2  and  b  are  the  only  terms  left 
to  be  explicitely  defined.  The  former  can  be  found  in  many  texts  (see 
e.g.  Ro  59,  page  4)  and  we  have  in  the  program  used  one  of  the  simplest 
forms,  namely: 


v'ZT+T 


J>  =  (-) 


3  2 


1  “l1/2 

(2j+l )  (2tj  '  +  1 ) 

4tt 


( j  j •  I -I 

VJJ  2  2 


20)  if  £  +  £ '  is  even 


=  0 


if  £  +  £'  is  odd 
(18) 


Detailed  calculations  of  the  reduced  matrix  elements  of  b  can  be 
found  in  (Ch  54)  and  (Fo  55).  Their  values  are  tabulated  in  Choudhury's 
paper  (Ch  54)  and  have  been  used  in  this  program. 


2.  Description  of  the  program 

The  program  CORPAR  I  calculates  the  matrix  elements  of  the  total 
Hamiltonian  HT  for  given  values  of  I,  £,  Ej  and  tiu,  diagonalizes  the 
matrix  so  obtained  and  yields  energy  eigenvalues.  The  eigenvalues  are 


written  as: 
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E  IM> 


I 

j  f  NR 


A  ( E )  |  j  NR  IM> 

J  ,IMK 


(19) 


The  expansion  coefficients  A.  (E)  are  printed  out.  The  basis 

3  fN , R 

states  |j  NR  IM>  are  written  as  |j  N  R>  in  the  program  output. 

Fig.  1  shows  a  schematic  flow  diagram  of  the  program.  This 
diagram  introduces  the  main  subroutines.  The  remaining  subroutines 
are  used  to  calculate  Racah  coefficients,  Clebsch-Gordan  coefficients, 

Y2  reduced  matrix  elements,  etc... 

2.1  Input  Output  Data 

The  single  particle  states  are  defined  by  their  total  angular 
momentum  j  which  can  take  a  maximum  of  three  values  corresponding  to 
the  maximum  of  three  orbitals  available.  At  this  point,  it  must  be 
again  emphasized  that  the  orbital  angular  momenta,  £,  MUST  have  the 
same  parity.  The  separation  energy  between  the  single  particle  states 
are  defined  by  the  variable  EPS.  The  phonon  energy  is  called  EFONON 
(in  MeV)  and  the  coupling  parameter  XI(£)  is  defined  by  XMIN,  XMAX  and 
XSTEP  for  its  minimum,  maximum  and  step  value  respectively. 

The  coupling  constant  k  has  been  fixed  to  40  MeV,  as  suggested  by 
Bohr  and  Mottelson.  This  value  is  commonly  used  by  most  authors  (Ch  54, 
Pa  70).  For  given  £  andHoo,  the  value  of  the  nuclear  deformabi  1  i ty  C2 
(CTWO  in  the  program)  is  calculated.  Finally,  for  comparison  purposes, 
two  other  coupling  variables  are  calculated:  X0  is  the  x  used  by  Choud- 
hury  (Ch  64)  and  HARVEY  Y  VALUE  is  the  coupling  parameter  introduced  by 
Harvey  (Ha  63) . 
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Figure  1 

Program  Flowchart 


In  order  to  simplify  the  typing  and  page-setting,  the  use  of  flow¬ 
chart  symbols  has  been  greatly  reduced. 


FLOWCHART 


COMMENTS 


INIT(l)  initializes  everything 


INPUT  reads  all  parameters 
necessary  for  a  set  of  calcula¬ 
tions 


MIXJ  defines  all  possible  combina¬ 
tions  of  j  and  R  to  get  I 


INIT(2)  initializes  the  energy 
matrix 


1 
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DIAGME  calculates  diagonal 
matrix  elements 

NOD I AG  calculates  off-diagonal 
matrix  elements 

HERDAG  is  the  diagonalization 
subroutine 
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2.2  Input  Format 


Card  # 

Variable  and  Format 

Comments 

1 

ITAG ,KPUNCH  (211) 

ITAG  =  1  Start  a  calculation 

i  1  END 

KPUNCH  -  1  All  results  for  the 
following  calculation 
to  be  punched  out 
(see  CORPAR  II) 

2 

TITLE 

80  alphameric  characters 

3 

I (K)  [K=l  ,6]  (6F10.5) 

Read  spin  values  of  the  final 
states 

4 

J,J1  ,J2  ,EPS1  ,EPS2 ,EF0N0N  (6F10.5) 

Defined  in  the  text 

5 

XMIN ,XMAX ,XSTEP  (3F10.5) 

For  one  value  of  £  only,  make 
XMAX=XSTEP=0 . 

Example : 

Set  of  calculations  for  N  =  29  nuclei 


Card  # 

Vari  able 

Comments 

1 

ITAG  =  1 

2 

The  following  calculations  for  N =  29  odd  nuclei 

3 

0.5  1.5  2.5  3.5  4.5  5.5 

I  ranging  from  1/2  to 
11/2 

4 

1.5  0.5  2.5  2.02  2.2  1.5 

0 

£i  =  P i~p 3  -  2.02  MeV 

e2  =  f 5" P.3  =  2-2  MeV 

2  2 

fiu>  =  1.5  MeV 

5 

2.5  4.0  0.5 

£  ranging  from  2.0  to 
4.0  by  steps  of  0.5 

. 


■ 

’.01  / 

. 

13 


2.3  Corresponding  Output 

THE  FOLLOWING  CALCULATIONS  FOR  N=29  ODD  NUCLEI 

FIRST  ORBITAL  J=1.5  SECOND  ORBITAL  J=0,5  THIRD  ORBITAL  J=2.5 

FINAL  SPIN  STATES:  0.5  1.5  2.5  3.5  4.5  5.5 

EPSIL0N1  =  2 o 020  MeV  EPSIL0N2=2. 200  MeV  PHONON  ENERGY=1.500  MeV 

XMIN=2 . 500  XMAX=4. 000  XSTEP=0.0 

(on  your  next  page) 

******NUMBER  OF  BASIS  STATES  IS  10 

BLOCK  DATA  INDEX  KM  1  THIS  IS  FOR  SPIN  0.500 

COUPLING  PARAMETER:  XI=2.500  CTW0=1 35 , 8123  MeV  X0=0 .7217  HARVEY  Y  VALUE=0.4431 
EIGENVALUES  IN  MeV 

8.3230  7.3489  6.8185  5.5221  5.0847  3.9581  3.2850  2.5516  1.9603  -0.1727 

EIGENVECTORS  BUILT  ON  THE  |N,R,J>  STATES 

1 0 ,0 , l/2>  1 1 ,2 ,3/2>  1 1 ,2 ,5/2>  |2,0,l/2>  |2,2,3/2>  |2,2,5/2>  |3,0,l/2>  |3,2,3/2>... 

-0.095  etc.... 

• 

e 

o 

o 

• 

0.5352  0.6673  0.3295  0.1978  -0.2871  etc.... 

Remarks : 

i)  The  index  KM  is  to  be  used  with  the  program  CORPAR  II. 

ii)  One  can  see  in  this  case  that  the  lowest  eigenvalue  for  spin  I  =  1/2 
has  an  energy  -0.1727  MeV,  Its  wave  function  can  be  written  as 

|-0.1 727(l/2)>  =  0.5352 | 0,0,l/2>  +  0.6673| 1 ,2 ,3/2>  +  0.3295 | 1 ,2, 5/2>  +  ... 

iii)  Some  calculations  of  the  same  type  sometimes  include  a  least  squares 
fitting  procedure  for  the  energy  levels.  It  was  felt  (after  using  this 
method  in  conjunction  with  Nilsson  Model  calculations)  that  this  kind  of 
procedure  is  somehow  uncertain  or  results  in  lengthy  calculations  caused 
by  slow  convergence.  The  former  case  occurs  when  too  few  levels  are 


. 


' 
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being  fitted.  They  are  then  heavily  weighted  and  the  overall  result 
is  usually  poor.  On  the  other  hand,  when  too  many  levels  are  being 
fitted,  it  appears  that  the  convergence  is  extremely  slow  when  more 
than  one  free  parameter  is  being  used. 

2.4  Suggestions  for  the  Use  of  this  Program 

To  simplify  the  search,  it  is  advised  to  normalize  all  the  energy 
terms  with  respect  to  the  phonon  energy.  For  different  values  of  EPS1 
and  EPS2  (in  fiu)  units),  one  should  plot  the  results  for  a  reasonable 
span  of  £.  Once  the  ranges  of  5,  EPS1  and  EPS2  are  chosen,  the  calcula¬ 
tions  can  be  achieved  with  absolute  energy  units  in  a  limited  number  of 
runs . 


2 . 5  Machine  Resources 

e.g.  6  spin  (from  1/2  ->  11/2)  and  £  stepped  from  0  to  5  by  =  0.5 
this  corresponds  to  a  total  of  66  di agonal i zations 

TIME:  <  2mn 
LINES:  <  4K 
SPACE:  100K 
CLASS:  A 

If  5  is  not  stepped,  there  are  only  6  di agonal i zations  i.e.  for  a  6  spin 
value  case,  one  could  run  under  class  Q. 


- 
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CHAPTER  II 


CORPAR  II 


Title:  CORPAR  II 
Type:  Main 

Source  Language:  Fortran  IV (G ) 

Origin:  T.P.G.  Carola,  Nuclear  Research  Centre 

Abstract:  Using  the  wave  functions  calculated  by  "CORPAR  I",  this 


program  calculates  electromagnetic  transitions,  lifetimes, 
mixing  ratios  and  static  moments. 


1.  Description  and  theory 

In  the  following,  the  reduced  transition  probabilities  will  be 
explicitly  derived  for  magnetic  dipole  and  electric  quadrupole  transi¬ 
tion.  The  derivation  of  the  static  nuclear  moment  formulae  will  be 
omitted  as  it  is  basically  equivalent  to  the  B(M1)  and  B(E2)  calculations. 

1.1  Electromagnetic  Transitions 

The  transition  probability  T(x)  for  radiation  of  a  photon  of  multipole 
order  X  and  frequency  ooA  is  given  by  (B1  52) 


B( x )  is  the  reduced  transition  probability  defined  by: 


' 

P  9 
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B(X)  =  2irr  l  j<I  lM'  |M(X ,y )  |  IM>  | 2  (2) 

MM'  y 

where  M(x,y)  is  a  multi  pole  operator  of  order  X;  |  IM>  is  the  initial 
state  and  1 1 1 M ' >  is  the  final  state. 

The  values  of  the  multi  pole  operators  used  here  are  those  introduced 
by  Bohr  and  Mottelson  (Bo  53)  for  the  system  resulting  of  the  coupling  of 
a  single  particle  to  the  quadrupole  oscillations  of  the  core  surface. 

The  multi  pole  operators  are: 

i)  Electric  quadrupole  operators:^ 


MP(2,y)  = 


<eP  +  eeff>  r2  ^(e.  +  )  +  R*  a*  (3) 


where  ep  is  equal  to  e  for  protons  and  o  for  neutrons,  e  ff  is  an 
effective  charge  parameter,  r  the  radial  coordinate,  RQ  the  nuclear  radius. 
Z,  e  have  their  usual  meaning.  All  the  remaining  parameters  have  been 
defined  in  Chapter  I.  The  effective  charge  parameter,  eeff,  has  been 

added  here  in  place  of  the  -  \  term  introduced  by  Bohr  and  Mottelson. 

A2 

This  term,  which  corresponds  to  "recoil  effects",  is  sometimes  introduced 
in  calculations  but  it  has  little  theoretical  support.  As  a  matter  of 
fact,  some  authors  have  shown  that  this  term  should  cancel  out  in  complete 
calculations  (Fa  59). 


f  * 

In  the  literature  one  usually  finds  Y_  instead  of  Y_ 

2  y  2y 


However,  for 

phase  consistency,  it  is  necessary  to  use  the  complex  conjugate  of  the 

spherical  harmonics.  We  recall  that:  Y*  (0,<f>)  =  (-)^  Y  _  (0,cf>). 

2y  2 ,  y 


I 


. 


. 
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ii)  Magnetic  dipole  operators: 


Mmo,p) 


e~h 

2Mc 


£ 

z 


S 

z 


(4) 


g  and  g  are  the  orbital  and  spin  g-factors  respectively.  The  core 

X/  S 

•  efi 

gyromagnetic  ratio  gR  is  taken  as  Z/A,and  i s  the  nuclear  magneton. 

1.2  Calculation  of  the  B(E2) 

We  write  the  basis  vectors  of  the  coupled  system  in  the  angular 
momentum  space: 


j  NR  IM>  =  l  (j  R  mm' |IM)  |j  m> | NR  m'>  (5) 

m,m' 


Introducing  the  operator  (3)  into  equation  (2)  and  making  use  of  (5), 
it  can  be  seen  that  B(E2)  will  have  the  form: 


B  ( E  2 ) 


I  l 

NR  j'  N'R' 


2 


(aA  +  3B  +  yC) 


with  a,  3  and  y  constants  to  be  defined  later. 
A,  B  and  C  are  given  by: 

A  =  < j 1  N'R'  I'M' I Y*  I j  NR  IM> 

1  2\i 1 

B  e  < j  1  N'R'  I'M'  | b  | j  NR  IM> 

C  =  <  j  1  N'R'  I'M'  I  (-)yb*  I  j  NR  IM> 

1  y 1 


(6) 


(7) 


. 


fit. 
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1.2.1  Calculation  of  A 

A  E  <j'  N'R'  I'M' | Y*  |j  NR  IM>  =  <j '  N ' R '  I ' M ' | ( - ) y  Y_  _  |j  NR  IM> 

2U  2  f  —  y 

Using  (5) 

A  =  l  (j1  R'  vv' | I ' M ' ) ( j  R  mm' |IM)  < j '  v|(-)y  Y.  _  |j  m>  <N'R'  v 1 |NR  m‘> 

vv'  2'  y 

mm' 

Using  the  orthonormality  of  the  basis  vectors  | NR  m'>  and  applying  the 
Wi gner-Eckart  theorem,  one  gets: 

A  =  l  SNN'  {rr'  V.'lj'  R'  vv '  1 1 1 M 1 ) ( j  R  mm1 1 IH) 

VV  ' 
mm' 

x  (-)y(j  2m-  u|j'v)  <j'  |  | Y2  |  |j> 

All  the  following  steps  are  self-explanatory.  (Basically  one  rearranges 
the  various  Clebsch-Gordan  coefficients  in  order  to  get  a  Racah  coeffi¬ 
cient).  (see  e.g.  Ro  57) 

A  =  l  ( j '  R'  v  m ' | I ' M ' ) ( j  R  mm ' |IM)(-)y(j  2  m  u | j '  v) 
vmm' 

x  <j'l|Y2||j>  &m, 


=  I  (-)3  ’_V  /Mir1  j'  M'  -  v|R  m')(-)R+j"I(R  j  m'm|IM) 

vmm'  “ 

(-)u(-)j'n/^~  O' 


I 


X 


j  V  -  m 1 2  -  p)  <j ■ 1 1 Y2 | | j>  6nn, 


. 


. 
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(21 '+1 ) (2j '+1 )  v  /_\R-j ,+j-i+j+j '+j-2 


A  H-) 

vmm' 


x  (1 1  j  ‘  M '  -v|R  m')(R  j  m'm|IM)(j'j  -  v  m|  2u)  <j  '  |  |  Y2 1  |  j>  6^,  6^, 


/(2I'  +  l)(2j'+l)  (-)R_I'j  W{I  ’  j  '  Ij;R2)  (I  ■  2  M'  y|IM)  <j'||Y,,||j>  <5NN,  «RR, 


(8) 


1.2.2  Calculation  of  B 


B  =  < j '  N '  R '  I '  M' |b  | j  NR  IM> 


The  method  is  the  same  as  in  the  calculation  of  A. 


I  ( j 1  R1  vv ' 1 1 1 M' ) ( j  R  mm' | IM)  <j  1  v  | jm>  <R' v ' |b  | R  m ' > 
vv'mm1  M 


l  6...  6  (j1  R'  vv' | I ' M 1 ) ( j  R  mm' | IM) (R '  2  v'  y|R  m'>  <N 1 R ' | | b | | NR> 

vv'mm'  33  Vm 


l  (j1  R1  m  v1 | I ' M ' ) ( j  R  mm' | IM) (R 


v  'mm' 


l  (-)R,+v,Jf£rr  (I'  8'  v'lJ 


v '  mm' 


X  (-) 


R'  -V 


V 


2R+1 


(R'R  v*  -m'  1 2 


2  v'  y | R  m ' )  <N ' R ' | |b| |NR> 

-  m) ( j  R  mm1 | IM) 

-  y )  <N'R'  | | b |  | NR>  6 


6  . 

33 


I 


i  y-  ('••  • «  r 1 1:  \<  n  (  m 
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B 


(21 '+1 ) (2R+1 ) 

”5(2j+l) 

x  (j  R 


j  ^_^R'+v'+j+R-I+R'-v' 


v  'mm' 


-m  - m'  1 1  -M)(R'R  v'  -m'  |2 


(I1  R'  -  M'  vr  |j  -  m) 
-y)  <N 1 R 1 | | b | |NR>  6 


I 


=  \/(2I'+l)(2R+l)  (-)2R'+j+R  1  W(I'  R’  IR;j2) (1 1 2  -M'  -y|I  -M) 

x  <N ' R ' I  I b I |NR>  fi, 

3D 


=  ^/ ( 2 1 '+1 ) (2R+1 )  (_)2r,+3+R  i+i'+2  I  W(j(I  R 1 R •  2J ) ( 1 1  2  M'  y|IM) 

x  <N ' R 1 | | b | | NR>  6  (9) 


1.2.3  Calculation  of  C 


C  =  <j'  N 1 R 1  IM |  (-)W  b*  |j  NR  IM> 

y 

=  ( - ) y  l  ( j 1  R '  w ' 1 1 1 M 1 ) ( j  R  mm ' | IM )  < j 1  v | j  m>  <R '  v '  |  b* | R  m ' > 

mxn'vv'  y 

=  (-)u  l  <5 . .  ,  6  ( j  '  R'  vv1  1 1  'M'  )(j  R  mm1  |  IM)  (R  2  m'  -y|R'  v ' )  <N 1  R'  |  | b* |  |NR> 

imm'vv'  Vm 

The  calculation  is  the  same  as  for  A  and  B.  The  final  result  is: 

C  =  y/(2I '+1 ) (2R'+1 )  (_)r,+2r+d-i+i,+2-i  W(I.i  R'R*2j)(I 1  2  M'  y | IM) 

x  <N'R' | |b*|NR>  6  (10) 
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1.2.4 

As  seen  in  eq.  (6),  the  B(E2)  will  be  proportional  to  a  sum  of  the 
square  of  (A+B+C). 

(A+B+C)  =  /(2I  '+1 )  (I'  2  M'  M-M 1  |  IM)  [^ST+T  W(jj  '  1 1 1  ;2R)  (- )R_I_j  <jl||V2||j> 

aNN'  6rr'  +  (42R+TT  (-)?K<+j+R"2I+I'+^/<N'Rl|ib|iNR> 

+  v^F+T  (-)R'+X+j-2I+r+X<N'R'  |  |  b*  |  |NR>}  S  W(R'R  1 1 1  ;2j )  ] 

(11) 

In  order  to  simplify  the  phase  factors,  the  COMPLETE  expression  can 
be  multiplied  by  any  phase  factor,  bearing  in  mind  this  expression  will 
ultimately  be  squared.  If  the  complete  expression  is  multiplied  by  (-)  , 

one  then  gets:  (Note  that  in  the  phase  factor  -R  or  +R  is  equivalent  as  R  is 
i  nteger) . 

(A+B+C)  =  /(2T  +  1)  (I'  2  M'  M-M 1  |  IM) 

x  [(-)I,+j'R  W(jj'  II';2R)  ^JriT<j'||Y2||j>  6nn,  6rri 

+  (-)I_j  !  (-)RI  /2R'+1  <N'R' I |b*| |NR> 

+  (-)R  /2R+T  < N ' R " | | b | | NR> j  W(R' R  I1 1 ;2j )  6 . . ,  ] 

(12) 

To  get  the  B(E2),  one  must  now  introduce  the  complete  quadrupole 
operator  (the  constants  have  been  left  out)  and  include  the  complete 


. 
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eigenvalues  functions  defined  by: 

|E  IM>  =  l  A  (E)|j  NR  IM>  (13)  (see  equation  (19)  in 
j,N,R,  :'N'R  part  I) 


The  notation  used  by  Choudhury  and  O'Dwyer  (Ch  67)  has  been  followed 
closely.  Keeping  with  these  authors,  the  following  terms  are  now  introduced 


<ep+eeff)  <r2> 


i  (e  +e  ) 

5  p  eff  0 


(14) 


3  “T\  OJ  7  q2 

47  ,/  2C  Ze  Ro 


(15) 


Using  eq.  (12)  and  the  expansion  (13),  together  with  the  definition  of 
B(E2),  one  finally  gets: 


B(E2)  =  (21 1 +1 ) 


i+i  ' 


.1  (E)  xW-DI,+J-R^<J 

jNR  J  J 

j'N'R' 


2  I  IJ 


W(jj  ■  1 1 '  ;2R)  Sm,  +  n2(-)I_j  [(-)*'  /2R'  +  1  <N'  R'  |  |b*|  |NR= 


+  (-)R  /2R+T  <N'R'  I  |b  I  |NR>]  W(R'R  1 1 1  ;2j)  6  j3 ,  (16) 


where  the  relation  l  1(1'  2  M'  M-M'|IM)|2  =  21+1  has  been  used.  The 

MM'  ' 

6  ,  term  has  been  introduced  for  the  sake  of  completeness,  (a  is  the 

orbital  angular  momentum  of  the  particle  in  the  state  |jm>). 


. 
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1.3  Calculation  of  B(M1) 

Following  the  same  method  as  described  in  the  case  of  the  calculation 
of  B(E2),  it  can  be  seen  that  three  terms  have  to  be  calculated  which  can 
be  labelled  as : 


A  E  <j' 

N1 

1  R ' 

I 

'  M ' 

1  & 

Z 

1  j 

NR 

IM> 

B  E  <j' 

N1 

1  R ' 

I 

'M ' 

lSz 

1  j 

NR 

IM> 

o 

III 

A 

C—i. 

N1 

1  R ' 

I 

1 M 1  | 

R 

z 

1  j 

NR 

IM> 

(17) 


1.3.1  Calculation  of  A 

One  starts  by  using  the  expansion  of  | j  NR  IM>  (eq.  5),  then  introduces 
an  expansion  in  the  (£ ,s)  space  in  order  to  get  the  eigenvalue  of  the  operator 
£  .  The  other  steps  deal  mainly  with  re-expressi ng  the  different  Clebsch- 

Z 

Gordan  coefficients  in  order  to  regroup  them  into  Racah  coefficients.  The 
different  steps  read  as  follows: 


A  =  l  (j 1  R1  vv 1 | I'M1 ) (j  R  mm' | IM)  <j'v|£z|jm>  <N 1 R 1  v 1 | NR  m ' > 

vv ' 

mm' 


I  ^NN ' 

VV  ' 
mm' 


RR' 


6  R' 


v  m 


vv 


I'M' ) ( j  R  mm'  j IM)  A' 


where  A'  =  < j 1  v|£z|jm> 

In  order  to  evaluate  A',  let  us  express  |jm>  in  terms  of  the  new  subspace 


. 
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jm>  =  l  (£  s  mm  |jm)  |£>|s> 

36  S 


m  m 
£  s 


A'  =  l  (z  s  m  m  | jm) (a  1  s'  m'  m'|j'v)  <£'  m'|£  |£  m  >  6 

36  S  36  S  36  Z  36  SS 


m.m 
£  s 

I  ' 

m„  m 
£  s 


The  last  term  is  easily  evaluated  using  the  Wi gner-Eckart  theorem 
(see  e.g.  Bo  69,  page  82). 


A'  =  l  (£  s  m£ms|jm)(£'  s'  m'0  m^|  j  1  v' )  (£l  m0  y  |  £ 1  m\)  <£'||£7!||£>  6e, 


mm 
£  s 

m' m' 

£  s 


£  s 


£  1  £  1  1  z 1  1  s' s 


=  J  (£  s  m£ms|jm)(£'  s  mpmc  |  j  1  v ' )  (£l  m?  y  |  £ '  mp )  /£(£+! )  l 


m.m 
£  s 


£  s 


£  1  £ ' 


£  £ 


£+s-j  £-m 

I  (-)  (s  £  msmjjm)(-) 


£ 


mm 
£  s 


mu 

2 


1/2 


(j '  £  m'-m  Is  m  ) 

£  s 


x  (-) 


£-m  +£+£-1 

X/ 


2£+l 


1/2 


(££  -  m£mjly)  /£(£+! ) 


r  £-j-l/2-2m 

1  (-)  (s  £  m  m  | jm) ( j '  £'  m'-m  Is  m  ) 

S  36  X/  o 

1/2 


mm 

£  s 


x  (££-m£m£ | 1 y ) 


(2,j  '  +  1 )  (2£+l ) 

3x2 


/£(£+!  ) 


A.  =  (_)£-j+l/2  /zyTT  /£(£+!  )(2£+l)  W(j  '  £  j  £  ;l/2  l)(j'lm  y|jm)  (18) 
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Now 

A  ■  I  ( 

vmm' 

X 

where  { 


Fi  nal  1  y ,  , 

A  =  (-) 

x 

1.3.2 
fact  that 

B  =  (-)* 


-)£_j+1/2  yzyrr  Au+i)(2£+i)  w(j'j  u;  1/2  1) 

{(j'  1  rn1  y  |  jm)  ( j  1  R'  vv'|rM')(j  R  mm '  |  IM)}  6nn  ,  6^,  (19) 


}  =  Z  ( J  '  R  ’  vv  '  |  IM ' )  ( j  R  mm1  |  IM)  ( j  '  1  vy|jm)  6  6 

.  I  JNJN  K.K 

vmm 


=  l  \  mrr  (-)j'  V(1'  J"  M'-v|R  m')(-)R+j  i(R  j  m'm  I  IM) 


vmm' 


2R+1 


x 


(_)  j '-v+j '+j-i  ty/y+T 


(j'j  -  V  m|l  p)  6hn,  Srk, 


=  (-) 


R-I+2 j+ j ' 


W( j  1  j  I  ’  I ;  1 R ) (I '  1  m'  y  |  IM) 


\  becomes 


:  i+Jt-l/2+j+j  '  /(2I.  +  1)  /(2j  '  +  1 )  ( 2j  i-l )  /£(£+1)(2£+l) 

W(  j  '  £  j  £;l/2  1)  W(j'j  I'I;lR)(r  1  M'  u  |  IM)  (20) 


The  calculation  of  the  B  term  is  equivalent,  but  one  uses  the 
s  =  1/2;  therefore  one  gets  : 

i+£-1/2  /2J-I+Y  /(2j'+l)(2j+l)  /3/2 

W(j 1  1/2  j  1/2;1  £ )  W(j 1 j  I' I ;IR) (1 1  1  H'  w 1 IM)  (21) 


X 


n 
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1.3.3  In  order  to  evaluate  C,  one  will  use  the  relation: 


<N ' R ' M ' | R | NRM>  =  /R(R+1)  (R  M  1  ylR'M') 


C  =  l  (j  R  mm' | IM) ( j '  R'  vv' |I'M')(R  1  m'  y | R'  v')  /R(R+1)  6^,  6rrI  6 


vv ' 

mm' 


■  I  (-) 

v ' 

mm' 


R  +V  ^Ij+l 1  0 '  R'  _M'  vlj  -  m  )(j  R  mm' | IM) 


x  (-)R+m'  /2£itl  (R'  R  v'  -  m 1 1  M)  /rTR+TJ  «nn,  «rr, 


_  v  I  sj+R-I.  2I'+1  .  / 2R+1  ,tl  „  . .  > 

-  H-)  V2jTrv~(I  R  - M  m  Ij  -  m> 


mm' 


x  (j  R  -  m  -  m'  |  I-M)  ( R  R  v'  m'|l  -  y)  v/R(R+l) 


=  (-)j+R_I  yrr+r  ymr  w(tr  ir;ji ) (i ■  i  -  m(  -  u | i-m)  /r(r+i) 

=  (-)I,+:I"R  STT+T  /R(R+1 )  (2R+1 )  W(I'R  IR;jl)(I'  1  M'  y|lM)  (22) 


In  order  to  get  B(M1),  one  now  proceeds  as  in  the  calculation  of  B(E2) 

(see  eq„  (13)  and  what  follows).  To  be  able  to  compare  the  complete  result 
with  the  result  of  Choudhury  and  O'Dwyer  (Ch  67),  it  is  necessary  to  multiply 
throughout  (before  squaring)  by  the  phase  (-)1  +I+2J„  It  is  then  possible 


This  difference  in  phase  -  which  does  NOT  affect  the  final  result  - 
apparently  comes  from  another  definition  of  the  basis  vectors,  namely: 

|j  NR  IM>  =  l  (Rj  m'mjlM)  |jm>|NRm'> 

Comparing  with  the  definition  used  here,  one  sees  that  the  Clebsch-Gordan 

coefficients  differ  by  a  phase:  R+i-I 

(Rj  m'mjlM)  =  (-)  J  (j  R  mm' | IM) 


. 


' 


. 
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to  square  and  sum  the  square  of  the  remaining  Clebsch-Gordan  coefficient. 
The  final  result  reads  as: 


B  ( Ml ) 
r*i ' 


_3_ 

4tt 


(uN)  (2I'+1) 


I 

j  j  'NR 


A.  (E) 

jNRI v  ' 


Aj ' NRI ' 


x  l(-)I,+R+i'1/2[(2j+l)(2j'+l)]>/2 

*  *  I  l/2 

X  [ga(-)3+:i '  [*.(A+1  )(2A+l )]  W(j'j  M;1  1/2) 
+  gs  v^72  W(jj  1  1/2  1/2;U)]  W(j'j  I '  I  ;1R)  «w, 

+  gR(-)I+j+R  tR(R+l )  (2R+1 ) 1 1/2  W(RR  I '  I  ;1  j )  «i(1  6..J 

eti 

where  uN  is  the  nuclear  magneton 


(23) 


1.4  Static  Moments 

The  method  used  to  calculate  the  static  moments  is  the  same  method  we 
have  used  in  all  the  previous  calculations.  The  results  will  now  merely 
be  stated: 


1.4.1  Electric  quadrupole  moment 

The  electric  quadrupole  moment  in  a  given  state  of  spin  I  and  energy 
E  is  defined  as:  (Bo  53,  Ch  54) 


<E  INI 


/  16tt  / 

v— ( 


e  +e  „)  r2  Y*  (e,*)  +  —  Z  R2  a?„|E  IH> 
P  eff  20 


'0  a20 


Q(E;I) 


M=I 


. 
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where  the  terms  of  this  equation  have  been  defined  earlier. 
Using  the  wave  functions  defined  by  eq .  (5),  one  gets: 


-  f  Ro  I  a,„t(e)  a,,„101t(e) 


'i  +  l  )  (21  +  3)  _.^R  j NRI 

j  ’N’R1 


j  NRI 


y¥(ep+eeff)(-)I+j'R%'  6RR'  ^T<j'l|Y2||j>  W(jj'  II;2R) 


+  z 


Hu) 


C(-)I_j  6...  6££,[(-)R'  >/2Rt+T  <N'R'  I  |b*|  |NR> 


+  (-)R  V/2R+T  <N'R'  |  |b|  |NR>]  W  ( RR '  1 1 ;  2  j ) 


(25) 


1.4.2  Magnetic  dipole  moment 


Similarly  the  magnetic  dipole  moment  y  is  defined  by 


y(E;I)  =  <E  IM I uN(ga  iz  +  gs  sz  +  gR  Rz)|E  IM>X= 


or 


y(EI) 


I=M 


=  y 


I (21+1 ) 


N\/  I  +  l  AjNRI^  Aj'N'R'I^  [9£^ 

j  'N'R' 


I+j-R 


6nn'  6rr'  5jj'  6 AA ' 


x  /R ( R+ 1 )  ( 2 R+ 1 )  W(RR  II  ;lj)  +  g£ (- )R’£+1/2+j+j  ' -i  ^  6rr| 


x  /(2j+l)(2j'+l)  /£(£+!  )(2£+l)  W  ( j  j  '  II ;  1 R)  W(££  jj';l  ±) 


+  9S(-) 


£+l/2+R-I 


6NN'  ^RR'  <5££' 


X  /372  /(2j+l)(2j'+lT  W(jj  1  II  ;1R)  W(jj 1  1/2  l/2;U)j 


(26) 
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2.  Description  of  the  program 

The  program  CORPAR  II  is  basically  simple.  It  reads  in  the  wave  functions 
output  by  CORPAR  I  and  two  main  subroutines  called  EMTRAN  and  STATIC  are 
used  to  calculate  the  B's  and  nuclear  moments  respectively.  The  remaining 
of  the  program  is  constituted  by  INPUT-OUTPUT  subroutines  and  various 
functions  to  be  used  with  EMTRAN  or  STATIC. 

2 . 1  Input-Output  Data 

In  order  to  facilitate  possible  expansions  of  this  program,  each  set 
of  data  to  be  read  is  preceded  by  a  card  identified  by  an  index  ITAG 
(Format  II ) . 

The  first  set  of  data  -  before  reading  wave  functions  -  is  preceded 
by  ITAG=1  and  contains:  Z,  atomic  number;  A,  atomic  mass;  AMI,  AM2,  AM3 
angular  momentum  of  the  odd  particle  in  the  respective  orbits  available; 

AN1 ,  AN2,  AN3  the  corresponding  orbital  angular  momenta  of  these  orbits 
(Format  8F10. 5) . 

Normally  these  values  are  input  once  and  for  all  (unless  calculations 
are  done  for  different  nuclei!).  Then  data  are  input  at  will,  being 
identified  by  the  control  index  ITAG. 

Other  cases: 

ITAG  =  0 

Each  BLOCK  DATA  output  by  CORPAR  I  (i.e.  for  a  complete  set  of 
wave  functions  for  given  £,  e1,e2,  “Ru  and  for  different  final  spin  I) 
must  be  preceeded  and  followed  by  a  blank  card. 


V 
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ITAG  =  3 


Read  EFF  (F10.5) 


If  EFF=0 


charge  is  e  for  proton,  o  for  neutron 


If  EFF^O  and  f  10  charge  is  (1+EFF)  e  for  proton,  EFF  x  e  for  neutron 


proton,  -  —  e  for  neutron 


If  EFF=10 


ITAG  =  4 


Read  Ti tie  (80  alphameric  characters) 

Read  DIFF,  KI ,  KF,  INDI,  INDF ,  NQ  [F10.5,  5(12, 8x)]. 

DIFF  is  the  energy  of  the  gamma  ray  for  the  transition  considered. 

The  KI  and  KF  coefficients  define  a  block  data  [i.e.  a  set  of  wave  functions 
. . for  a  given  final  spin  I].  These  coefficients  are  printed  by  CORPAR  I, 
The  INDI  and  INDF  coefficients  define  the  level  considered  with  each  block 
KI  and  KF  respectively  (see  example).  If  NQ  =  2,  only  static  moments  are 
cal  culated. 

ITAG  =  5 

Read  COEFK  and  CTWO  (2F10.5) 

If  these  values  are  not  entered  COEFK  =  k  is  set  at  40  MeV ; 

CTWO  =  C2  is  calculated  from  this  value. 

ITAG  =  6 


End  of  calculation. 


' 
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2.2  Example 

This  is  a  sample  input  corresponding  to  calculations  of  e-m  transi¬ 
tions  in  55Fe.  The  data  punched  out  by  CORPAR  I  (defined  by  DATA 
hereafter)  corresponded  to  £  =  2.75,  tuo  =  1.40  MeV,  e1  =  2.02  MeV, 
e2  =  2.2  MeV,  1/2  £  I  £  11/2. 


INPUT: 

ITAG  =  1 

26.  55.  1.5  0.5 

ITAG  =  0 

A 

DATA 

ITAG  =  0 
ITAG  =  3 
0.5 

ITAG  =  4 

55  FE  CALCULATES  STATIC 
0.  2  2  16 

[i.e.  KI  =  KF=2 ,  here  spin  1=3/2; 

ITAG  =  4 

55  FE  FIRST  1/2  TO  G.S. 

0.412  1  2 

[i.e.  KI  =  1 :  s pi n  1  =  1/2 ; 
calculated  here  the  1/2 

ITAG  =  6 


2.5  1.  1.  3. 


16  2 

I NDI= I NDF= 1 6 ,  here  ground  state; 
NQ=2:  calculate  static] 

16 

:  lowest  eigenvalue;  to  be 
3/2 (g  .s - )  0.412  MeV  transition] 


10 

I N  D I  =  1 0 


END 
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Corresponding  OUTPUT 

1st  Page: 

Atomic  number:  26.  Mass  number:  55. 

THE  G  FACTORS  USED  HERE  ARE: 

GL=0 .0  GS=-3.8263  GR=0.4727 

THE  FOLLOWING  CALCULATIONS  ARE  DONE  WITH  A  CHARGE  EQUAL  TO:  0.500  IN  UNITS  OF  E 
CALCULATIONS  FOR  XI=2.750  HBAR-0MEGA=1 .400  MeV 

END  OF  DATA  READING 

THE  PARTICLE  HAS  A  SPIN  J=3/2  IN  THE  GROUND  STATE.  THE  CORRESPONDING  L  IS  1 
THE  NEXT  ORBITS  AVAILABLE  ARE:  J2=l/2  WITH  L=1  AND  J3=5/2  WITH  L=3 

Next  Page : 

55FE  CALCULATES  STATIC 

THE  LEVEL  WITH  SPIN  1.500  HAS  A  QEI=-0.143  EOO  BARNS 
MAGNETIC  DIPOLE  MOMENT  -0.162  E01  N.M. 

55FE  FIRST  1/2  TO  G.S. 

KI,  KF,  INDI,  INDF  AND  NQ  ARE  RESPECTIVELY  1  2  10  16  0 

DECAY  FROM  THE  1=1/2  STATE  TO  THE  1=3/2  STATE  E (GAMMA)=0 . 41 2  MeV 

B ( Ml )  =  0.189E  01  (N.M.)SQ  T(M1)  =  0.233E  13  1/SEC  TAU1  =  0.430  E-12  SEC 

B(E2)  =  0.167E  03  E2-FM4  T(E2)  =  0.244E  10  1/SEC  TAU2  =  0.411  E-09  SEC 

DELTA  =  -  1.85  DEGREES 

DELSQ  =  0.105  E-02 

TAU  =  0.429  E-12  GAMMA  =  0.153  E-02  EV 


' 
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All  the  variable  names  have  their  usual  meaning.  The  B's  and  T's  have 
already  been  defined  (see  eq.  1  and  2).  TAU1  and  TAU2  are  the  corres¬ 
ponding  partial  lifetimes;  DELTA  is  the  mixing  ratio  in  degrees  (i.e. 
arctan  6);  DELSQ  is  62;  TAU  is  the  lifetime  (in  sec.)  for  the  level 
decaying  through  this  transition;  GAMMA  is  the  corresponding  width  in 
eV. 


2.3  Machine  Resources 

CORPAR  II  is  a  fast  program.  For  most  of  the  calculations,  i.e.  for 
20  to  30  transitions,  the  following  requirements  will  be  quite  sufficient 

TIME:  £  2mn 
LINES:  £  2K 
SPACE:  100K 
CLASS:  A 
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